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MEETING OF 1959 MARCH 13 
Professor R. O. Redman, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


Donald van-Burren Dun, 17 Victoria Gardens, London W.11 (proposed by 
W. B. Caunter); and 

Anthony Fynn, S.J., Riverview College Observatory, Riverview, Sydney, 
Australia. 


‘The election by the Council of the following Junior Members was duly con- 
firmed: 

‘Thornsteinn Saemundsson, University of London Observatory, Mill Hill 
Park, London, N.W.7 (proposed by C. W. Allen); 

Nigel Oscar Weiss, University Dept. of Geodesy and Geophysics, Cambridge 
(proposed by N. M. Hill); and 

Bernard Douglas Yallop, Mathematics Dept., Imperial College, London, 
S.W.7 (proposed by G. J. Whitrow). 


Seventy-seven presents were announced as having been received since the 
last meeting, including :— 

G. Alter: Two renaissance astronomers: David Gans — Joseph Delmedigo 
(presented by the author); 

\. Becvar: Atlas eclipticalis 1950.0 (presented by the Czech Academy of 
Sciences); and 

M. Waldmeier: Leitfaden der astronomischen Orts-und Zeitbestimmung, 2nd 
edition (presented by the author). 


MEETING OF 1959 APRIL 10 
Professor R. O. Redman, President, in the Chair 


‘The President announced that the Council had elected the following 
Associates of the Society: 

Nicholas Ulrich Mayall, Lick Observatory, Mount Hamilton, California, 
U.S.A.; 

William Wilson Morgan, Yerkes Observatory, Williams Bay, Wisconsin, 
U.S.A.; and 

Alexi Boriseivitch Severny, Crimean Astrophysical Observatory of Academy 
of Sciences U.S.S.R., Pochtove, Crimea, U.S.S.R. 
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The election by the Council of the following Fellows was duly confirmed :— 


Fazal Ahmed, Royal Observatory, Edinburgh (proposed by H. A. Briick); 

Edward Doylerush, Tyn Leydiant, Roewen, Conway, Caerns. (proposed by 
J. Davis); 

Jan A. Hégbom, 25 Victoria Street, Cambridge (proposed by M. Ryle); 

J. D. North, 35 North Way, Headington, Oxford (proposed by C. W. Allen); 

Dennis John Pearson, 40 Beechwood Road, Uplands, Swansea (proposed by 
B. E. Featherstone); 

John Roland Shakeshaft, Cavendish Laboratory, Cambridge (proposed by 
M. Ryle); 

Harry Sykes, Central Electricity Board, Butterworth, Malaya (proposed by 
A. F. O’D. Alexander); 

Arthur George Tearle, 26 Wilsden Avenue, Luton, Beds. (proposed by 
H. C. King); and 

Alan Burgess, Physics Department, University College, London, W.C.1 
(proposed by M. J. Seaton). 


The election by the Council of the following Junior Member was duly 
confirmed :— 


John Edward Basil Ponsonby, 26 Upper High Street, Thame, Oxon (proposed 
by J. Heywood). 


Eighty-five presents were announced as having been received since the last 
meeting, including :— 


Royal Astronomical Society: Franklin Adams Charts (presented by Mrs. 
Wilfred Merton); and 


G. J. Whitrow: The structure and evolution of the universe (presented by 
the author). 


MEETING OF 1959 MAY 8 
Professor R. O. Redman, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


Simon Fairfax Humphrey John Archer, Norman Lockyer Observatory, 
Sidmouth, Devon (proposed by D. R. Barber); 

Serge Colombo, Institut Henri-Poincaré, Paris, France (proposed by 
V. C. A. Ferraro); 

John Aldworth Gould, 24 Salisbury Street, Swindon (proposed by A. P. 
Lenham); 

Donald Harvey Tarling, 35 Thornhill Road, Sparkhill, Birmingham 11 
(proposed by J. M. Bruckshaw); and 


Wilfred Wallbank, 6 Waterside Cottages, Everton, Chorley, Lancs. (proposed 
by W, Granger), 
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The election by the Council of the following Junior Members was duly 
confirmed :— 


Petros Serghiou Florides, Royal Holloway College, Englefield Green, Surrey 
(proposed by W. H. McCrea); and 


Gillian Peach, Royal Holloway College, Englefield Green, Surrey (proposed 
by M. R. C. McDowell). 


Eighty-two presents were announced as having been received since the last 
meeting, including :— 
International Astronomical Union: Second conference on co-ordination of 
galactic research (presented by I.A.U.); and 


R. D. Davies and H. P. Palmer: Radio studies of the universe (presented by 
Routledge & Kegan Paul Ltd.). 





A NEW DETERMINATION OF THE CENTRE TO LIMB CHANGE 
IN SOLAR WAVE-LENGTHS 


M. G. Adam 


(Communicated by the Director, University Observatory, Oxford) 


(Received 1959 May 4) 


Summary 


The centre to limb increase in wave-length, which is observed for 
medium strength solar lines, occurs almost wholly over the last 10 per cent 
of the way across the disk. So far, measurements have not been made in this 
region with nearly high enough disk resolution. The Oxford 35 m telescope, 
which gives a solar image of radius 165 mm, is used to make detailed observa- 
tions of limb effect for 261 disk points in latitudes near to the solar equator. 
The observed points all lie close to the limb, at radial distances frorn the 
centre between 85 and 100 per cent of the disk radius. The individual 
observations show considerable fluctuations about a mean curve, which are 
attributed to the local velocity fields which are known to exist on the Sun. 
The systematic differences which appear in the measures from the east 
limb and the west limb, for radial distances greater than 98 per cent of the disk 
radius, are shown to be the result of light scattering. In the observational 
determination of the limb-effect curve it is most important that distortions, 
due to velocity fields and scattered light, should be taken into account. 

The new limb-effect curve is in fair agreement with the old determinations, 
over the range in common. At the extreme limb there is a shift with respect 
to the centre wave-lengths which, expressed as a Doppler velocity, amounts 
to 0:53 kmsec™'. This indicates at the extreme limb a red shift of 0-84 km 
sec~' compared with the Einstein shift of 0-636 km sec". 





Introduction.—F ollowing Evershed’s pioneer work on the change in wave- 
lengths across the solar disk (8), detailed observations have been made at Potsdam 
(9) and at Oxford (1), and these have formed the basis of theoretical discussions of 
the subject. Although the two series of observations are in good agreement, there 
are still peculiarities and deficiencies in the observational relation between wave- 
length and disk position, which make further measurements a matter of some 
importance for the problem of solar wave-lengths. On the visible disk the 
Fraunhofer line wave-lengths do not differ from their centre values by as much as 
0-001 A until 85-90 per cent of the way towards the limb, and after this they 
increase rapidly to become some 0-007 A longer than the centre wave-lengths. 
Fig. 1 shows this centre to limb behaviour, for medium-strength lines, as observed 
at Oxford in 1948. The ‘limb effect’ A(@*), by which we mean now wave-length 
observedat disk position &* minus wave-length observedatthe disk centre, isexpressed 
as a Doppler velocity in km sec! and is plotted against cosine 6*, where 6* is the 
angle between the line of sight and the solar radius at the point of observation. 
This representation of limb effect is strictly relevant to the radial streaming 
interpretation, and is more suitable for Freundlich’s path length interpretation 
than the usual plot of A(@*) with sin 9*, The vertical lines through the observed 
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points in Fig. 1 show the estimated uncertainty in A(@*), and the horizontal lines 
show the range in cosine 6* covered by each observed point. 

The first necessity in improving the observations is obviously to work with 
much higher resolution on the solar disk, and thus to reduce the range in cosine 
6* covered by each observed point. The Potsdam observations were made with a 
resolution probably comparable with that shown in Fig. 1, and the number of 
observed points isthe same. In making more precise observations it is important 
to take into account the recent work on motions in the solar photosphere (12), (12), 
(15), (17); this work shows the existence of long-lived small scale and large scale 
velocity fields on the Sun, which will superpose Doppler shifts on the limb effect 
observed at any point. Such velocity-field effects may be reduced by making 
observations at many close points on the disk, and by combining observations 
from several different regions on the Sun. 
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Centre Cosine @* Limb 
Fic. 1.—Solar limb effect as observed in 1948 
Wave-length region 6100 A. 14 lines: mean Rowland number 6. 


For the theoretical interpretation of limb effect, it is observations near to the 
solar limb, and particularly at 6*=7/2, which provide the most sensitive test. 
We see that in Fig. 1 the extreme point covers the range in cos 6* 0-25 — 0-00, and it 
is here especially that high disk resolution is desirable. But this 25 per cent range 
in cos * with which we are particularly concerned corresponds to a range of only 
3 per cent in the distance from limb to centre, and detailed observations can only 
be made in this region with a very large solar image. Moreover it is the region 
within which the disk intensity drops very sharply from 60 per cent to zero, and 
this means that the observations are very susceptible to any effects introduced by 
the scattering of light, instrumentally or by the Earth’s atmosphere. 

The Oxford 35 m telescope, with its large solar image, and the associated high- 
dispersion grating spectroscope form excellent equipment for a redetermination of 
the limb effect curve. The paper gives an account of the first results for solar 
latitudes near to the equator, and cosine 6* values between 0-5 and oo. The 
measures are differential and give only the change in wave-length from centre to 
limb. The new observational material and its measurement and reduction to give 
the observed limb effect are described in Section 1. These results show that the 
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scattering of light must be taken into account in deriving the limb effect for cosine 
6* < 0-2, and this problem is investigated in Section 2. The paper concludes with 
a short comparison of the final limb-effect curve with theoretical predictions. 


1. The observational material, measurement and reduction 


1.1. Solar spectra.—For the detailed point-to-point observations on the solar 
disk which we now wish to make, it is essential to adopt the method of measuring 
the solar lines against telluric lines in the same wave-length region. Only in this 
way can we avoid spurious displacements arising from any differences in collimator 
illumination. The new grating spectroscope gives sufficient dispersion to work in 
the red at 6300 A with atmospheric oxygen lines as standards. This isa favourable 
region in which to measure medium-strength solar lines, and the constant strength 
of the oxygen lines gives them a considerable advantage over the water vapour 
lines for our purpose. 

The 35 m telescope (16) gives an image 330 mm in diameter, corresponding to 
a scale 5”°9 per mm; this was used in conjunction with the 11 m spectroscope in 
May and June 1957 to obtain solar spectra at 6300 A for the determination of limb 
effect. ‘The image was guided so that the spectroscope slit was always radial, and 
each plate normally carried two spectra, one in the upper hemisphere of the solar 
image and one in the lower hemisphere. The two spectra were for points equi- 


distant from the disk centre and lying closely along the same diameter. Spectra 
from the disk centre were obtained from time to time during the observing period, 
and also some spectra, described in Section 2, for estimating the amount of 
scattered light. 70-80 plates were obtained covering the full range in cosine 6* 
across the disk. For the plates which included the actual limb the rocking 


plate guiding device was used to secure as sharp a limb as possible. Ilford 
half-tone panchromatic plates (R50) were used for the limb effect and centre 
spectra and these were developed in Parkhurst two-solution developer; Ilford 
R 40 panchromatic plates and Press Contrast developer were used for the scattered 
light plates. 

Because of the rotation of the solar image the full series of plates covers a wide 
range of latitudes in the two hemispheres. For a first investigation measurements 
have been confined to six heavily exposed limb spectra which were especially 
favourable for measurements close to the extreme limb. These were obtained 
fairly close to the solar equator, and include the limb and about 20 mm on to the 
disk; the range in cosine 6* is o-o to o-5. The details for these six spectra are 
given in Table I: ten centre spectra obtained on seven days have also been 
measured, and three of the scattered light spectra 


TasBLe I 
Details of measured limb plates 
Plate No. Date U.T. 1957 Exp. Region Solar latitude in degrees 
18a May 22°5906 40 Secs. East limb — 11°53 to — 12°36 
18b 22°5946 40 secs. West limb + 11°09 to + 13°73 
19a 22'6036 40 Secs. West limb + 13°85 to + 16°92 
19b 22°6081 40 secs. East limb — 16°99 to — 18°62 
20a 22°6165 40 secs. East limb — 19°66 to — 21°70 
20b 22°6211 40 secs. West limb + 19°23 to +23°21 


The Bausch and Lomb grating was used in the spectroscope in the first order, 
giving a dispersion of 1-4mm A~! at 6300A, or in terms of Doppler velocities 
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33°98 km sec! per mm shift on the plate. Following preliminary measurements 
in the region it was decided to work with seven lines, four oxygen lines to serve as 
standards, and three solar lines. The selected lines are shown in Table II: the 
oxygen lines alternate with the solar lines and span the whole wave-length range: 
An analysis of the measurements on the east and west limbs, showed that these 
oxygen lines were apparently free from blends with unsuspected solar lines. The 
I.A.U. 1928 tables contain wave-lengths for six of the seven lines: the wave-length 
given for the oxygen line 6302-0 A is a mean value from Babcock’s measures (4) 
and the Allegeheny Observatory—Bureau of Standards wave-lengths (6). The 
Rowland numbers of the lines are given in column three, and for the iron lines we 
have also shown the multiplet number and lower state excitation potential as 
given in Miss Moore’s Multiplet Table (1945). 
i 
Taste Il ' 
Measured lines 
1.A.U. wave-length Element Row. Mult. Excit. pot. (low) Transition 
inA No: No: in volts 
6295°960 Atm. O, 


6297°799 62 2°21 a®*P-y*D* 
6299°228 


6301°508 816 3°64 z§Pe-D 
[6302000] 


6302°499 816 3°67 2*P*-fD 
6302°764 


1.2. Micrometer measurements and the determination of line shifts.—The spectra 
have been measured with a Hilger photomeasuring micrometer which is fitted with 
a transverse motion on the plate stage (15); this motion greatly facilitates the 
measurements at many close points in the same spectrum. All the measurements 
of centre, limb, and scattered light spectra have been made with a diaphragm in the 
focal plane of the eye-piece, which limits to 0-145 mm the extent of the spectrum 
which is visible, perpendicular to the dispersion. 

The measurements in the limb spectra were designed to give a uniform spread 
in cosine 6* values. This means that the measured points were very close together 
near to the extreme limb and at gradually increasing intervals further in on the disk. 
The limb plates were marked, perpendicular to the dispersion, with fine scratches 
at mm intervals from the inner edge of the spectrum. Two sets of marks were 
made, one to the red of the wave-length region at 6337 A and the other to the violet 
at 6278A. The limb formed the other edge of the spectrum, and its position was 
estimated visually with respect to the mm marks. The results from the two ends 
showed an uncertainty of the order of 0-04 mm in its position. In a normally 
exposed spectrum Miss Hart (11) has found that visual estimates agree well with 
photometric determinations. But, as Dr Treanor has pointed out, there may be 
systematic displacements with heavy exposures. The R 50 plates, however, have 
a very steep characteristic curve, which would reduce any such effect and the fair 
agreement of the estimates from the two ends of the plate, where the photographic 
density is noticeably different, indicate that it is small. Measurements 
to the dispersion were made as soon as the photographic density at the limb edge of 
the spectrum permitted it; this gave measurements for four or five points in the 
scattered sky light beyond the estimated limb. The first 10 measured points 
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were at 0-1 mm intervals, the next 10 at o-2 mm intervals, then 10 at o-5 mm inter- 
vals and the remaining 13 or 14 at 1mm intervals. This arrangement gives a 
fairly uniform spread in cosine 6* values; we now have 24 measured points in the 
region covered by the first observed point in Fig. 1. For the six limb spectra 
measurements were made in all at 261 points on the disk. 

Disk positions were calculated in terms of p the radial distance of the measured 
point from the disk centre, and py the disk radius. We obtain p,-p from the 
estimated limb position and the particular point of measurement, both known with 
respect to the mm scratches. The image radius p, at the time of observation was 
obtained from a series of direct measures on the solar image which gave 
Py = 167°1 mm in the spectroscope slit plane, for an angular radius of 16’. The 
spectroscope magnification can be assumed to be unity. 

In deducing line shifts from the micrometer measurements we have followed 
the procedure as described by Professor Plaskett (15). The measurements in 
the centre spectra were used to set up ‘standard’ micrometer readings for the 
seven lines. The solar lines were first of all corrected for the observer’s motion 
with respect to the solar centre, and the mean reading from the ten spectra was 
then taken as the standard setting for any line. For the solar lines this reduces the 
effect of velocity fields on the Sun; the standard error of the mean reading for a 
solar line was equivalent to a velocity of 0-04 km sec~! (average value for the three 
lines). This estimate includes measurement errors in the solar and oxygen lines 
as well as any real variations due to velocity fields, though estimates of the measure- 
ment errors showed that the scatter was chiefly due to field effects. 

For the limb spectra the micrometer measurements of the oxygen lines in each 
spectrum were used to set up a correction curve for the spectrum, to reduce the 
solar measurements to the dispersion of the standard settings. ‘The correction 
curve was assumed to be linear over the very short wave-length range with 
which we are concerned. The corrected shifts for the solar lines were then 
coverted into Doppler velocities with factors 34-00, 33:99 and 33°97 kmsec~! per 
mm shift for the three lines in the order of Table II. The three velocity shifts 
were then used to form a mean velocity shift for the disk point. In order to point 
out the main features of the observational problem, we show in Fig. 2 the mean 
velocities so obtained for the east limb spectrum and the west limb spectrum, from 
Plate 18 in Table I. The velocities shown include solar rotation and limb effect 
and also effects due to the observer’s motion; these last change only very 
slightly over the observed range in p/p,. On the east limb we see that solar 
rotation and limb effect combine to give a nearly constant observed velocity, 
whereas on the west limb both effects give a red shift, and the observed velocity 
increases steeply as we approach the limb. An outstanding feature of the 
observations is the fluctuation, which increases with distance from the limb, about 
the general form of the two curves. In the light of the work on photospheric 
motions which we have quoted earlier, it appears that we are again observing the 
phenomenon of velocity fields on the solar surface. An analysis of the 
departures of the individual line velocities from the mean for any one disk 
point, gives 0-08 km sec~' for the standard error of a single line velocity. From an 
examination of these residuals we conclude that there are no systematic differences 
in the centre-limb behaviour of individual lines within + 0-03 kmsec~!. Assuming 
that measurement errors in the solar and oxygen lines are comparable in size, we 
then estimate fluctuations due to this cause at 0-06 km sec~, (standard error of a 
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single disk-point velocity). We see that the observed fluctuations are definitely 


greater than this and must indicate real variations from point to point on the solar 
disk. 
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"1G. 2.—Velocities observed on the east and west solar limbs (Plate 18). 
Open circles*—west limb. 
Filled circles—east limb. 

1.3. Correction for Earth’s motion and solar rotation.—In the correction for 
Earth’s motion and the calculations for solar rotation we have followed very 
closely the scheme set out by Professor Plaskett (14). In almost all these calcu- 
lations it was possible to work at representative values of p/p, across the slit, for 
the epoch of exposure, and then use interpolated values for the 44 measured points. 
The corrections for Earth’s motion were found to vary linearly across the slit, 
and all the disk-point velocities obtained as described in the previous section were 
first of all reduced to velocities on the Sun (v,,,), by applying this correction. 

We assume that v,,, at any point X on the disk is due wholly to solar rotation 
and limb effect A(6*), The component of solar rotation along the line of sight 
will be R, cos y, where R, is the linear velocity of solar rotation at the latitude B 
of X, and y, is the angle between the line of sight and the direction of solar rotation 
at X. Then from Section 3-3 in (14) we have 


cos y, sin (L — L,) cos By sin 8*cosec 6 . (1) 
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in which (L—L,) is the difference between the solar longitude of X, and L, the 
longitude of the centre of the projected disk, and B, is the latitude of the disk 
centre, @ is the angle between the line of sight to the disk centre and the solar radius 
tothe point X. If we then use sin 6* ~ p/p, and sin (L — L,) =sin @ sin (P—¢) sec B 
in equation (1) we obtain 

cos B cos y, = sin (P—¢) cos B,. p/py (2) 


¢ and P are the position angles of X and of the N point of the disk respectively. 
¢ we obtain for the epoch of observation from the known orientation of the solar 
image. A series of observations of the E-W line at different times of day gave 
$ =231'2 +h for the upper limb, and less 180° for the lower limb, in which & is the 
hour angle of the Sun in degrees. If we could assume solid body rotation for 
the Sun, R, would be R, cos B, where Ry is the linear velocity at the equator, and 
the rotation correction we require would then be given by the product of R, and 
the right hand side of equation (2), and would vary only with p/po, across the slit. 

Although the departures from solid body rotation are small for our latitude 
range of 12°—23°, they are not negligible and must be taken into account. Instead 
of R, = R, cos B we have, with a Faye type formula, R, = R, cos B (1 — (b/R,) sin*B) 
in which 6 is an observational constant. For b/R,, on the basis of determinations 
by Adams and Lasby (3), and the revision of their values by Freundlich, Brunn and 
Briick (9), we have taken a value 0-229. The line-of-sight component of solar 
rotation v, is then given by 

v,= R,, cos y, = R, cos B (1 — 0:229 sin*B) cos y, 

or in terms of equation (2) 


v,= R, [(1 — 0-229 sin*B) sin (P—¢) cos By. p/py| =f. Ro. (3) 


For any one slit position the factor f varies only with p/p, and the quantity 
0-229 sin?B. Across any one slit position 0-229 sin?B varies by at most o-o1 and its 
maximum value for the highest latitude used is 0-035. For each slit position a 
curve was constructed from which the value of 0-229 sin?B could be read off for 
each p/pp. 


We have then 261 observational equations of the type 
Vovs, =f. Ro + A (0*) 


in which the unknowns are R, and A(@*). * A(@*) is a function of p/pp, the form of 
which we are investigating, and in order to obtain R, we must combine equations 
at the same value of p/p). In order to reduce velocity-field effects and to use the 
directly observed velocities we have made solutions for R, combining observations 
over short ranges in p/pp. We know that A(@*) varies only very slowly for our 
inner disk points, and further out, where A(@*) changes more rapidly, we must 
restrict ourselves to very short ranges in p/pp. We then obtain R, from each range 
and a mean value for A(@*) over the range. Table III shows the result of least 
squares solutions of the observational equations for R, and A(@*); the range in 
p/p) we have given in terms of cosine 0*, and we have also shown the number of 
observational equations used from each limb. In the lower half of the table we 
give the results obtained very close to the limb: the number of observational 
equations is reduced because of the small range in cos #*. As we shall see in the 
next section it appears that these results, obtained so near to the limb, are affected 
by scattered light. ‘The purpose of our present calculations is to obtain a value for 
R, which we can then use to find A(@*) with high disk resolution: we see that the 
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Taare III 
Least squares solutions for R, and A(0*) 
Range in cos 6* No. of obs: vel: R, km sec A(@*) km sec 

E Ww 
0°47-0°42 15 15 2°002 "156 
0°41-0'°35 15 15 1'986 0°186 
0°33-0°28 15 15 1°992 0248 
0'28-0°22 15 15 2°007 0°354 

1°997 


0°08-0'06 8 7 1°946 0°384 
0°06-0'04 4 3 1'941 0-398 
0°04~-0°02 3 3 1-918 0°456 
0°02-0'00 2 II 1*QOI 0°443 
mean values of A(@*) in column 4 of Table ITI increase fairly steadily in the way we 
should expect. For R, we have adopted the value 1-997kmsec™', the mean 
value from the first four groups in the table: these values come from a region of the 
disk unlikely to be affected by scattered light, and where A(6*) changes only slowly. 


The individual values are each based on 30 observational equations and they agree 
well together. 





km sec* 


0-6 


Ate") 








| 
0-4 O+3 0.2 
Cosine 0* 
Fic. 3.—The observed limb effect. 
Filled circles—east limb. 





1.4. The observed limb effect.—With this value of R, ar the appropriate f 
factors we have derived a value of A(6*) for each of our 261 points. The resulting 
limb effect curves for the east limb and the west limb are shown in Fig. 3. For 
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each limb the values of A(@*) from the three plates have been grouped together in 
intervals of 0-02 in cosine @* to give a mean value for the interval. Comparing 
Fig. 3 with Fig. 2 we see that the velocity-field effects are much reduced by this 
combination of the observations. From cos #*=0°5 to 0-2 the east and west limb 
curves agree well together, but the anomalous feature of the new measures is the 
steadily increasing difference between the two curves from cos 9*=0:2 to 0-0, 
which amounts to 0-2 km sec~ at the extreme limb. 


2. The effect of scattered light on the observed curve 

Further in on the disk, the differences between the east and west limb observa- 
tions in Fig. 3 might be attributed to local velocity fields persisting in the mean from 
three different disk positions. But from cos 6* = 0-2 to the extreme limb we cover 
1*4 x 10°km on the Sun, which corresponds to 190” at the disk centre, so that in 
view of their estimated sizes, we must be averaging out field effects. The plates 
18, 19, 20 were taken in that order over an interval of 40 minutes, during which the 
sky conditions, although reasonably good, were noted to be deteriorating, and we 
find that the departures of the west limb curves from the east limb curves increase 
in the same order. There seems little doubt that the peculiarities in the observa- 
tions for points further out than cos 6* = 0-2, must be the result of scattered light. 

We saw in Section 1.2 that on the east limb solar rotation and limb effect give 
shifts in opposite directions, and combine to give a velocity (of approach) which is 
almost constant from cosine 6*=0-2 outwards, with a slight decrease at the 
extreme limb. With such a distribution scattering of light gives very little change 
in the velocities. On the west limb solar rotation and limb effect both act in the 
same direction, both increasing outwards, and there is a sharp rise in the resultant 
velocity as cos 6* approaches zero. At the limb, and especially in the presence of 
limb darkening, such a distribution is distorted by scattering to give a less steep 
gradient for the velocity, disk position curve. 


- 2.1. Determination of the scattering function.—In order to test this explanation 
and to arrive at the true limb effect, we need to know the scattering function 
appropriate to our observations. We have obtained this from velocity measures 
in the ‘scattered light’ plates which were obtained from time to time during the 
observing period. These plates were obtained with the arrangement described 
by Miss Hart (11): a rotating sector was placed immediately in front of the slit 
to cover the lower half, and the solar image was guided so that the upper limb fell 
over the sector. The limb intensity was thus reduced to 0-035 of its normal value, 
while the adjacent sky spectrum was recorded with full intensity using the upper 
half of the slit. The slit was always radial to the solar image. In these spectra 
velocities could be determined in the sky spectrum for 10” to 70” beyond the limb, 
and these velocities enable us to determine the scattering function, independently 
of any photometric procedure, by measurements precisely similar to those used 
in the determination of limb effect. 

Three scattered light spectra were measured and reduced to give velocity shifts 
in exactly the same way as the six limb spectra. The measurements were made at 
three points on the disk and nine or ten points in the sky spectrum, the furthest 
measured point being 76” from the solar limb. Fig. 4 shows for two of the plates 
the measured velocities as a function of p/p,; the shifts are corrected for Earth’s 
motion, assuming that the scattered sky light requires the same correction as the 
limb points on the solar disk. The smooth curves show, for each of the plates, the 
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velocity distribution on the disk as calculated from a combination of solar rotation 
v, and limb effect A(@*) at each point. vv, was obtained using the value of R, as in 
Section 1-3 assuming, for our present purpose, that we have solid body rotation: 
as a best approximation to the true value, A(@*) was taken as observed on the east 
limb. The scattered light spectra are all for the west limb, approximately 30° 
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Fic. 4.—Velocity shifts at the limb and in the sky spectrum. 
Half filled circles—continucus line: Plate 95. 

Half filled circles—dotted line: Plate 96. 

Open circles: calculated points using the deduced scattering function. 


from the equator, and the two velocities are additive. We seek then a scattering 
function which, from the unsmeared velocities of the smooth curves, will produce 
the distribution shown by the observed points. The plates 95 and 96 were 
regarded as approximately typical for sky conditions during the observing period, 
while the third plate referred to more than usual scattered light. We see that the 
velocities observed on the disk fall below what we should predict in the absence 
of scattering by amounts comparable with the difference between the two limbs 
shown in Fig.3. The third plate showed that the departure from the predicted 
curve becomes greater as the sky conditions became worse. 

To find the effect of scattering on the measured velocities we assume that the 
velocity observed at any disk position is a weighted mean velocity of the light 
scattered to that position; the velocity contributions from neighbouring points 
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being weighted according to the disk intensity at those points. This is an assump- 
tion that has been tested observationally by Miss Hart(xz1). Our observed points 
are all fairly close to the limb and we therefore neglect the curvature of the image. 
If I(x) is the intensity at a point whose radial distance from the limb, is x and u(x) 
is the line-of-sight velocity at that point, we assume that J(x) and u(x) remain 
unaltered perpendicular to the x direction for the fairly short range with which we 
are concerned. Ifthe scattering processes have the effect that, for distances mea- 
sured in the x direction, the light at any point on the disk is redistributed in distance 


s according to the function S(s) the velocity a(r) measured at a point distant r from 
the limb will be 


a I(x) u(x) S(r — x) dx 
a(r)= 1 
| I(x) S(r —x) dx 





(4) 


where 2p, is the solar diameter. 

The general form of the scattering function S(s) is fairly well known; there is a 
region 2"~5” on either side of the peak within which the scattered light has a 
Gaussian distribution, followed by a ‘tail’ within which the light distribution 
falls off much more gradually. Our measurements of velocities in the scattered 
sky light make it fairly easy to find the form of the ‘tail’, but the Gaussian ‘core’ 
is more difficult, since we rely on the few and rather widely scattered limb velocities. 
The form of the ‘core’ is not, however, an important matter in considering the 
modification of the velocity distribution on the disk, since the effects come 
chiefly from the more distant scattered light. For J(x) in (4), the limb 
darkening, we have used the intensity distribution given by ten Bruggencate, 
Gollnow and Jager (5) which is corrected for scattered light, combined 
with the observations of Moll, Burger and van der Bilt (13) for points 
further in on the disk. Plates 95 and 96 were used together in deriving the 
scattering function; the sky velocities for these plates show no systematic differ- 
ences, and for u(x) in (4) we have therefore used the mean of the continuous 
curves shown in Fig. 4 for these two plates. Our numerical integrations in the 
expression for a(r) have all been made by summations using intervals of 0-002 in 
p/pp. The scattering function we deduce is a Gaussian curve (half, half-width 
3"°5) out to 5” from the peak value, after which it decreases more slowly to 1 per 
cent of the peak intensity at a distance of 30”. We have carried the ‘tail’ out to 
110” from the peak and then assumed that just under 1 per cent (0-78 per cent) of 
the total light is scattered on each side beyond this. In deducing the light 
arriving at any point on account of scattering, we have assumed that this small 
fraction of light, scattered from great distances, has the mean properties of light 
from the centre of the disk. 

We first of all tested the deduced scattering function by applying it back again 
to the plates 95 and 96, using again the mean predicted u(x) curve for these plates. 
The open circles in Fig. 4 show the velocities we then expect to observe in the 
scattered light plates. We see that the agreement with the disk point observations 
is quite fair, and that in the sky spectrum we can account for the observations very 
well out to 50 or 60"; the departures after this are presumably due to the artificial 
simplification of the function for large distances. 

2.2. Application of the scattering function: the final limb effect curve.—As a 
first step we have so far assumed that the east-limb observed velocities are affected 
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to a negligible extent by scattered light. We must now apply our scattering func- 
tion to test this assumption, and if necessary proceed to the true limb effect by 
iteration. On the east limb we have accordingly calculated ar) in equation (4) 
over the whole observed range in p/p», using for u(x) the mean of the observed 
velocities from the east limb spectra, after each has been corrected for the Earth’s 
motion. u(x) had to be extended beyond the range of the observations for these 
calculations, but this is a fairly straightforward matter. It is found that scattering 
has very little effect on the velocity distribution, the differences between a(r) and 
u(x) not exceeding 0-02 km sec~! at any part of the observed range. 

We next proceed in a similar way for the west limb. For u(x) we now use the 
sum of A(@*) as given by the east limb curve, and a mean value for the three spectra 
(1-2) of the line of sight velocity due to solar rotation. We find that the velocity dis- 
tribution is now considerably modified by scattering. At the extreme limb there is 
a difference of 0-18 km sec~! between u(x) and a(r), the difference decreases quite 
rapidly as we move in from the limb, so that at p/py=0-98 the difference is only 
o-04kmsec~!. For limb effect, however, this last 2 per cent of the way across the 
disk is of considerable interest; it corresponds to cos 6* =0-:20—0-00, which is 
just the region of theoretical importance. If our scattered light explanation is 
correct this calculated u(r) should now agree with the west limb observed velocities. 
The agreement between calculation and observation is shown in Fig. 5, in which the 
continuous line shows the limb-effect we should expect to observe for west limb 
disk points according to our calculations for the effect of scattered light. We see 


that the calculated curve gives quite a fair representation of the limb effect as 
measured on the west limb. 


centre to limb change in solar wave-lengths 471 





km. sec™* 


0.5 


A(e*) 








| L 


0.5 0.4 0.3 02 0.1 0.0 
Cosine @* Limb 





Fic. 5.—The calculated and observed limb effect for disk positions on the west limb, 


Ideally we should now invert the order of our calculations to pass from the 
observed a(r) to the unsmeared velocity distribution u(x) on each limb, The 
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calculations above show, however, that such a procedure would give limb-effect 
curves for the two limbs in substantial agreement, since light scattering accounts 
adequately for the differences between the two observed curves. In view of our 
ignorance of the exact scattering function appropriate to the actual limb spectra, 
we shall do best to accept now the limb effect as given by the east-limb spectra, the 
errors we predict are only of the order of 0-02 km sec™', and we cannot correct for 
these with any reality. 





km. sec™ ] 
0.6r- 








| | 
0.5 0.2 0.1 0.0 
Cosine 6° Limb 
Fic. 6.—The final limb effect curve, as observed on the east limb. 
Dotted line: 1948 observations. 


Taste IV 
The solar limb effect curve for equatorial regions close to the limb 
/Po A(@*) km sec P/Po A(@*) km sec=* 
1-000 0°526 0"950 0'263 
0'995 0°442 0°925 0°206 
0990 0°407 0"g00 o'1s9 
0985 0°381 0875 o-118 
0°980 0°358 den et 
o’O75 9°339 — sail 
0970 o°321 0°000 0°000 





The east-limb observations from Fig. 3 are shown again in Fig. 6, and in 
Table 1V we give, as a function of p/p, the values of A(@*) which are obtained 
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from a smooth curve drawn through the observed points in the figure. 
The 1948 limb-effect curve is also shown on Fig. 6 for comparison with the new 
determination. ‘There is fair agreement between the two series of observations 
over the range in common, though the earlier curve appears to be displaced by 
0-05 to o:10kmsec~ from the new measures. Such a displacement could arise 
from errors in the centre wave-lengths which are quite comparable with their 
estimated uncertainty. From Section 1-2 we have now an estimated uncertainty 
of 0-04 km sec! for the centre spectra (since field effects will be the same for all 
three lines). The 1948 centre measures were for 11 spectra taken on three days 
and the stated probable errors indicate an uncertainty of apparently the same 


magnitude as we find now: together these uncertainties might well account for the 
observed displacement. 


Discussion and conclusion. 


We can most easily compare our observations with theoretical predictions by 
considering the extreme limb values which have now been obtained for the first 
time. According to Freundlich’s hypothesis the centre-to-limb shift of a solar 
line depends upon the path length of the photon in the radiation field of the solar 
atmosphere. On this theory Freundlich and Forbes (10) predict a limb effect 
A(m/2) of approximately 1-3 km sec~!, which is more than twice the observed value ; 
moreover the form of the A(6@*) curve, which we can now follow to the extreme 
limb, is quite at variance with the path-length interpretation. 

On the radial current theory the red shift of a solar line, (solar wave-length 
minus vacuum-arc wave-length), arises from the combination of the Relativity 
shift and shifts due to radial streaming motions in the solar atmosphere (18). At 
the extreme limb any radial motions will be perpendicular to the line of sight and 
we must therefore observe the exact Relativity shift. From Table IV we see that 
the solar wave-lengths at the extreme limb already exceed the centre values by an 
amount equivalent to 0-53 kmsec~", which, from our general knowledge of centre 
wave-lengths, strongly suggests that the extreme limb red shift is well over the 
Relativity value. Centre red shifts have not been measured in this investigation, 
but according to the relation between centre red shift and line strength which we 
have obtained previously (2), we should expect a shift of 0-45 km sec~! or 0-009 A, 
though admittedly there is considerable scatter in the observational relation. We 
have more precise information since all three lines appear in the I.A.U. standard 
solar wave-lengths, and vacuum-arc wave-lengths have been measured for them in 
the Allegeheny Observatory—Bureau of Standards work (7). These wave- 


TaBLe V 
Centre red shifts of measured lines 
Wave-length El. Row. I1.A.U. Solar wave- A.O.-B.S. Vac. are Ay 
(A) No. length wave-length (A) 
6297°8 Fe 5 ‘799 ‘790 0-009 
6301°5 Fe Y -508 “501 0°007 
6302°5 Fe 5 "499 "493 0-006 


lengths and the resulting centre red shifts A, are shown in Table V: the mean 
centre red shift is 0-007 A. - We can make some comparison from our own mea- 


sures with both series of wave-lengths. Using oxygen lines as standards the mean 
wave-lengths of the three solar lines from our 10 centre spectra average 0-002 A 
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shorter than the I.A.U. values. Our vacuum-arc lines, measured by the method of 
circular channels agreed with the A.O.-B.S. values at 6020 A and were 0-002 A 
shorter at 6500 A, these differences would lead to a mean A, of o-006A. If we 
take o-0065 A or 0-31 kmsec~ as a likely value for the centre red shift, then with 
the extreme limb effect of 0-53 kmsec™, the solar red shift at 6* = 72/2 is o-84km 
sec~!, which is well over the Relativity value (0-636kmsec~'). This excess is 
quite inexplicable on the radial current theory. 

In conclusion we must point out, that although the discrepancies between 
observation and theory which we find are marked enough to be decisive, our ‘ final’ 
limb effect curve is only a small part of the work that is needed to arrive at the correct 
explanation of solar red shifts and limb effect. The chief use of the present 
investigation is probably in showing the difficulties in finding the centre to limb 
change in solar wave-lengths. In the first place the programme must be planned 
to allow for the existence of velocity fields on the Sun, and also, from the most 
recent work (17), for organised tangential motions on the solar surface. Secondly, 
allowance must be made for the effect of scattered light, which affects those regions 
of the disk which are most important theoretically. We find that at the extreme 
limb errors amounting to o-2kmsec~' occur at the west equator, under quite 
fair seeing conditions; at the two poles the estimated error would be o-1 kmsec=}, 
while it is practically zero at the east equator. With the possible exception of east- 
limb equatorial disk points, the appropriate corrections must be made in order to 
find the true limb effect. The limb effect may well be different for different solar 
latitudes, and vary with element, line strength, excitation potential and other 
physical properties of the absorbing atoms; all but the most outstanding of such 
effects are likely to be undetectable unless the real nature of the observational 
problem is taken into account. 


I should like to thank Professor Plaskett for his generous help throughout this 


investigation, and also Dr A. B. Hart for many helpful discussions on the problems 
of scattered light. 
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Oxford: 
1959 April 16 
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Summary 


The technique of interpretation of solar granulation observations is 
examined with special reference to the stochastic nature of the physical 
phenomenon, and to the absence in any optical image of sufficient 
information to allow the intensity distribution in the object to be explicitly 
inferred. Many of the points made are elementary, and so not essentially 
novel, but the need for explicit discussion and illustration of each one has 
been evidenced by the publication in the recent literature of results which 
do not have the significance attributed to them. Application is made to the 
simultaneous velocity and brightness observations of Richardson and 
Schwarzschild and of Plaskett. It is concluded that the evidence for corre- 
lation between velocity and brightness is weaker than has been supposed. 
Harmonically related periodicities in two sets of the Oxford observations 
may possibly be connected with magneto-hydrodynamic effects, but the 
suggestion of some authors that there are significant periodicities over large 
regions of the Sun is not confirmed. The investigation has been planned as a 
contribution towards bringing techniques of interpretation into line with 
current advances in the observation of granulation, and the p 
developed for carrying out the required calculations on the EDSAC II 
computer may form a basis for the rapid analysis of future observations. 





1. Introduction 


It has become increasingly apparent that for many purposes the assumption 
of plane stratification of the outer layers of the Sun is no longer a sufficiently 
good approximation. The formation of spectral lines, for example, cannot be 
treated accurately by a theory which ignores the inhomogeneities which are 
observed as solar granulation; and other phenomena, such as the structure of 
the chromosphere, and possibly the supply of energy to the corona, may be 
essentially dependent on processes of which granulation is the outward and 
visible sign. 

Granulation has accordingly received increasing attention in recent years. 
The work has three obvious aspects; namely, the observations themselves, the 
interpretation of the observations, and the incorporation of the information 
thus gained into theoretical solar models. The present paper is concerned with 
the second step, the reduction of the raw observational data to a form suitable 
for theoretical purposes. 

The initial difficulty is, of course, that even the largest granular elements 
are near to the limit of resolution set by ordinary daytime seeing. The problems 
which this creates are enhanced by the circumstance that conditions on the Sun 
are so different from those well known to us that we have, to begin with, no model 
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in terms of which to interpret granulation. A second circumstance which creates 
special problems is the stochastic nature of the phenomenon. 

The first of these circumstances is typical of astronomy generally. In ordinary 
landscape and portrait photography the existence of previous knowledge normally 
makes it possible to interpret correctly detail which is near to the limit of the 
optical resolution. This ability is specially developed in interpreters of air 
reconnaissance photographs, and it is closely associated with the fact that the 
observer in everyday living has the opportunity to handle and examine at close 
quarters the objects represented. This opportunity does not exist in astronomy, 
and accordingly it is only when all significant details are well resolved that it is 
safe to assume the reality in the object of the structures which suggest themselves 
to an observer examining the image. Moreover, it is difficult to establish that 
this condition is fulfilled, since to trust the visual impression that it is so is to 
assume just what requires to be investigated. When the structures in the object 
are not fully resolved, it is obviously misleading to treat measurements made on 
the image as if they had been made on the object itself, and explicit account has 
to be taken of the specific way in which the optical system transmits information 
about the object. 

An optical system of given aperture completely suppresses some information 
about the object, and other information is presented in the image in a 
metamorphosed form. The most clear and compact way of expressing these 
relationships is in terms of the two-dimensional Fourier transforms of the object 
and image distributions (Fellgett and Linfoot 1955, Linfoot 1958). Some 
spatial frequencies in the object are not reproduced in the image, and others are 
reproduced with changed amplitude; if aberrations are present, their phase 
may also be changed. Once the effect of optical imaging on the spatial frequencies 
has been formulated, its effect on other functions of the image brightness, for 
example the auto-correlation function, can be obtained by using the ordinary 
techniques of Fourier theory, and this will be illustrated in what follows. 

The difficulties of interpretation are much reduced when it is permissible to 
assume that the object structures are geometrically simple, or that they can be 
interpreted in terms of structures known on the Earth. For example, star fields 
can be taken as effectively consisting of point sources, and the major features 
of the lunar surface can probably be interpreted in terms of a rocky landscape. 
In other cases, where no such previously known pattern is applicable, the 
reservations implied by the theoretical relations correspond very closely to 
practical difficulties. The surface features of Mars are a well-known example 
of this, and the effort to interpret granulation in terms of terrestrial experience 
has traditionally led to comparison with rice grains or willow leaves. A further 
example may be provided by the comparisons between different surveys of radio 
sources. The apparatus used in the various surveys is responsive to different 
ranges of spatial frequencies in the brightness distribution over the sky, and 
indeed these ranges do not overlap in general. The resultant records have been 
interpreted in terms of the distribution of point sources which most simply 
explains each one: it is not obvious that the distributions so found ought to 
agree with each other, except for the brightest sources which are effectively 
resolved individually. 

The physically stochastic nature of granulation implies that techniques of 
interpretation appropriate in more deterministic circumstances are not necessarily 
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well adapted to granulation studies. It may be useful to map a star-field in 
detail, because the stars remain accessible to subsequent observation, but there 
is no reason to believe that the presence of a ‘‘ granule’’ in the solar disk at a 
particular place and time ordinarily tells us, in itself, anything of general interest, 
or allows any prediction to be made that can be checked by subsequent observation. 
It would in any case be impracticable, even if it were desirable, to record from 
instant to instant all the changing details of the granular pattern over the Sun. 
Interest centres in those properties of granulation which may reasonably be 
supposed to be characteristic of the phenomenon as a whole, and which enable 
predictions to be made which can later be confronted with further observations. 
It is accordingly necessary in the reductions to eliminate, as far as possible, those 
parameters which pertain only to the particular sample observed, and to retain 
those which are statistical averages in the sense that, as more and more data are 
included, they may be expected to approach some limit which is independent 
of the particular set of data used. Of course this separation cannot be made 
cleanly on the basis of any finite set of observations, and it is necessary to consider 
the uncertainties which result from the use of a limited sample. An example is 
discussed in Section 2.4. 

Several authors have previously discussed the precautions that need to be 
taken in the interpretation of observations of granulation. Uberoi (1955) has 
put forward serious objections to some of the methods that have been used, and 
indeed there seems little that can be urged against the substance of his criticisms. 
Wlerick (1957) has presented detailed illustrations of how the method of calculation 
used can introduce artefacts into the interpretation of the data. 

These questions are taken up in Section 2 of the present paper, where several 
of the more important topics, including some also treated by Uberoi, Wlerick 
and others, are discussed theoretically and illustrated. Numerical and photo- 
graphic experiments are presented which show what happens when random 
data are subjected to processes that have been used in the interpretation of 
granulation. It is evidently necessary to know this before one can assess the 
reality of structures which appear to be present in actual granulation data treated 
in the same way. The implications of this discussion lead to the choice, in 
Section 3, of methods of calculation aimed at interpreting the brightness and 
velocity changes observed by Richardson and Schwarzschild (1950) and by 
Plaskett (1954). The results of these calculations, which constitute the main 
computational part of the paper, are presented in graphical form in Figs. 7-11 
and their significance is briefly discussed in Section 4. 

It is not the purpose of the present paper to urge any particular view as to the 
exact nature of solar granulation, but some comment on this question seems 
appropriate in order to avoid misunderstanding. The appearance of a discrete 
granular structure is common when a ‘‘random’’ brightness-distribution is 
photographed with limited resolution and at high contrast, as has been pointed 
out by Uberoi and is further illustrated in Section 2.1. It follows that this 
appearance in images of the solar disk is insufficient in itself to prove that any 
such structure exists on the Sun, and in particular granulation photographs 
reproduced at high and unstated contrast are useless as evidence on this point. 
Uberoi concluded that all the measurements discussed by him were consistent, 
so far as they went, with a uniform turbulence spectrum. The older observations 
did not, in fact, establish that granulation was anything else than a random 
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brightness distribution with statistically uniform spatial power spectrum. From 
the purely observational point of view, the fact that the patterns seen in images 
of granulation have suggested discrete convection cells to modern observers is 
on exactly the same footing as the fact that in the last century they suggested 
comparison with rice or willows; the latter comparisons are evidently to be 
taken no more literally than is the appearance of bright rings round stars seen in 
a small telescope. From a more general point of view, the convection theory 
acquires special status from the theoretical prediction that the region of the Sun 
immediately underlying that from which most of the visible light comes will be 
convectively unstable, and from the independent observational evidence that the 
formation of spectral lines is affected by turbulence. It does not follow that the 
convection will be in discrete cells. Ordinary thermal convection at a given 
spatial wave-number occurs when the lapse rate exceeds the adiabatic rate by an 
amount which is, roughly speaking, an increasing function of wave-number, 
so that convection begins at nearly the lowest wave-number which is geometrically 
possible. For a layer confined between rigid horizontal planes, the circulation 
accordingly begins in cells having dimensions comparable with the distance 
between the planes; in an atmosphere the cell size depends on the scale-height. 
Magnetic or coriolis forces may modify the relative ease of convection at different 
wave-numbers. If convection is only just initiated, only small velocities develop 
and the motion is almost entirely in discrete cells corresponding to the most 
favourable wave-number. ‘The amount of heat that can be transported in this 
way is small, and this motion prevails only over a narrow range of temperature 
gradient. It is accordingly hydrodynamically rather unlikely, though possible, 
that solar granulation represents this sort of convection. If the driving gradient 
much exceeds the minimum needed to initiate convection, the layer becomes 
unstable over a wider range of wave-numbers, and the motion at low wave-numbers 
becomes fast enough to be turbulent, so that energy is fed from lower to higher 
wave-numbers until dissipated by viscosity. The spectroscopic evidence for 
microturbulence in the regions where lines are formed suggests that this is indeed 
what happens on the Sun. Under these conditions it remains true that the bulk 
of the energy of motion resides in the lower wave-numbers present. When the 
resolution in an image is sufficient to show only these wave-numbers, the pattern 
seen, because it lacks higher wave-numbers, resembles the pattern of discrete 
cellular convection in which higher wave-numbers are absent from the actual 
motion. One difference which is observable without attaining higher resolution 
is that in cellular convection the pattern can persist for a long time relative to the 
time of circulation in a cell, whereas in turbulent convection the large-scale 
pattern changes more quickly. 

Recent observations of granulation appear to show genuine resolution of 
non-random features. The visual observations of Thiessen (1955) were made 
with adequate controls to show when real features had been resolved, the only 
reservations being those inherent in the subjective nature of visual work. An 
important increase in photographic resolution has been attained in the series of 
observations made at high altitude by Blackwell, Dewhirst and Dollfus (19572, b) 
using a manned balloon of the traditional type, and by Schwarzschild (Rogerson 
1958) using an unmanned plastic balloon. Schwarzschild’s photographs, in 
particular, generally confirm Thiessen’s conclusions. Valuable additional 
information has been obtained from ground level by Résch (1956); and in the 
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massive and carefully planned photographic observations of Leighton (1957) 
who has studied the statistical time-changes in the granular pattern. 

The structures thus resolved and made available for study accord well with the 
idea that convection is the underlying cause. The relationship of the observed 
structures to convection processes in the Sun is necessarily a matter of detailed 
quantitative study. There are some difficulties in supposing that we see down 
to the actual convective region, and Whitney (1958), for example, discusses this 
problem and makes alternative proposals. For the reasons given above, the 
structures now resolved may be expected to be themselves composed of smaller 
elements, including those corresponding to the spectroscopic ‘‘ microturbulence ”” 
(Fig. 1(c)). A detailed knowledge of the properties of optical blurring and of 
the numerical processes used to interpret the data seems likely to be even more 
essential to the correct interpretation of the more sophisticated observations 
now becoming available than it has been in the past. 


2. Properties of optical blurring and of data-reduction methods applied to granulation 


2.1. Granular structures produced from random objects.—We have already 
referred to the discrete granular structures often seen when a field of randomly 
varying brightness is photographed with limited resolution and at high contrast. 
An effect of this kind is well seen when random noise is viewed on a television 
monitor screen. The screen appears covered with a closely-packed array of 
bright spots, changing from instant to instant. If the monitor is de-focused, the 
number of spots (and their contrast) decreases, while their size increases so that 
they continue to fill the whole of the scanned area of the screen. This behaviour 
is closely analogous to the so-called ‘‘5 sec-’’ ‘‘2 sec-’’ and “‘1 sec-”’ granulation 
of the Sun viewed under different seeing conditions. ; 

Skumanick (1955) concluded that the granulation images measured by him 
showed no significant black-white asymmetry. Since the probability distribution 
of the brightness variations at each point tends towards a symmetrical Gaussian 
form in an image where the largest individual object structures are not resolved, 
there may nevertheless be real asymmetry in images of higher resolution than 
were used by Skumanick. Dr D. E. Blackwell has very kindly allowed me to 
examine the original negatives which he obtained from an altitude of 18000 ft. 
When the photographic calibrations, which are available for these negatives, are 
taken into account, it seems quite definite that the images can be interpreted as 
bright granules separated by relatively narrow dark lanes, but not the other 
way round. This is of course inherently probable, since if, for example, 
turbulence on the Sun gave rise to statistically symmetrical temperature fluctua- 
tions, then the emission of visible light would have a skew probability distribution. 

Another characteristic of granulation images is that they have a reticulated 
appearance whereby the cells seem to form ordered patterns over a distance of 
several cell diameters. 

I have tried to reproduce these two properties of granulation experimentally, 
starting from fine photographic grain which can be taken as ‘‘random’’ for 
present purposes. A skew intensity distribution can readily be obtained by 
using the non-linear photographic. characteristic. The reticulated appearance 
is associated with a spatial spectrum in which the structural power in the object 
increases with frequency. When combined with the instrumental cut-off, this 
may cause the spatial spectrum of the image to be more or less strongly peaked 
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at some frequency. A rising structural power-spectrum can be obtained photo- 
graphically by superposing a positive and a negative taken of the same object 
but with different resolutions. Fig. 1 (a) shows the result of repeatedly copying 
photographic grain using these techniques. If the figure is held at some distance, 
a marked granular structure appears, whereas close examination shows that 
this is not in fact the predominant structure present. In Fig. 1 (6) the attempt 
has been made to reproduce a blurred image of Fig. 1(a) showing granular 
structure. It is not to be expected that with ordinary uncontrolled darkroom 
techniques it would readily be possible to match quantitatively the appearance 
of actual solar granulation. It can be said, however, that the original print of 
Fig. 1 (b) shows asymmetry and reticulation in the sense aimed at; and indeed 
several astronomers to whom it was handed without comment implicitly 
identified it as a negative of solar granulation. This shows that it is possible, 
starting from a ‘“‘random”’ grain distribution, to reproduce at least some of the 
characteristics of solar granulation ‘‘ sans aucune intervention de la main 
humaine ’’ (Janssen 1896). 

It may be mentioned here that the available prints of the excellent granulation 
picture obtained by Schwarzschild show remarkably large granules (> 1"). 
The sharpness of the apparent granule-boundaries, and their high contrast, 
provide immediate evidence that these pictures are of higher resolution than 
has been attained previously, and the large size of the granules is especially 
apparent in these parts of the published pictures which are sharpest according 
to these criteria. This suggests that in some previous pictures the granular 
structure may be broken up by multiple imaging of the kind often observed in 
star images under bad seeing conditions. It underlines the incorrectness of 
judging resolution solely by the size of the finest structures seen. 

2.2. The non-significance of guessed frequencies.—Experience of the necessary 
calculations shows that only in the simplest cases can the frequency-components 
of an observed function or its correlation function be found by inspection. 
A guessed ‘‘ quasi-frequency ’’’ may mean nothing at all, or it may refer to the 
cut-off above which frequency components have greatly reduced amplitude. 
If, as often happens, this cut-off is instrumental, the guess then gives no 
information about the phenomenon under investigation (see Section 2.5). 

As an example of an effect of this kind, the function sine x = sin (7x)/7x has 
an ‘‘obvious’’ oscillation with frequency 0-5. But its Fourier transform is in 
fact 


1, |t|<: 


rect (t) eS 


(9, | >: 


the apparent frequency turns out to be the cut-off. Sinc?x has an ‘‘ obvious’”’ 
oscillation at unit frequency; its transform is in fact 


tri(t)= {1 Ith Iel< 
lo, \t|>1 
which has a cut-off but no discontinuity at t= 1. 

2.3. Apparent serial correlation caused by blurring.—lf a quantity f(x,y), 
say the brightness over part of the solar disk, is observed optically, the finite 
resolution of the observation necessarily introduces correlation between the 
observed values at neighbouring (x,y). The presence of such correlation does 
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(a) Object made from photographic grain by repeated photographic copying under such conditions 
(see text) as to give a spatial power spectrum which is an increasing function of frequency over a 
range of frequencies, and an asymmetric probability distribution of brightness from point to point 
of the object. It may be regarded as a (negative) model of a possible brightness distribution on the 
Sun. 


(b) Blurred reproduction of object (a). This reproduction may be regarded as analagous to a 
(negative) photograph of the solar granulation blurred by seeing. It has a reticulated appear- 
ance, and is covered with “ granules’? which show on the negative as dark spots separated by 
narrower “‘ intergranular’’ spaces. This granular structure is not a predominant feature of the 
object shown in (a), but the two pictures appear more similar when viewed from some distance. 


Fic. 1.-—Resemblances to granulation in photographs derived from random objects (a) and (b), and 
in cloud formations (c), 
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(c) Stratus cloud layer showing small-scale turbulence within larger polygonal structures. 
If this layer were photographed with limited resolution it could give the appearance of smooth 
structureless granules. The mechanism of optical contrast is of course quite different here 
from on the Sun. 


Fic. 1. (contd.)—Resemblances to granulation in photcgraphs derived from random objects (a) and (b), and 
in cloud formations (c). 


Peter Fellgett, On the interpretation of solar granulation 
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not therefore of itself indicate that f is similarly correlated. The observed quantity 


is in fact 
fow=[ { ; 


w(E,n)f(x—£,9—n) dé dy (2.31) 


the convolution of f with the weighting or ‘‘ apparatus ’’ function w representing 
the resolution of the observation. From the convolution theorem of Fourier 
transforms it follows (see Section 2.5) that if f were random, in the sense of being 
itself serially uncorrelated and therefore having a uniform power spectrum, the 
Observations f*w would have the same auto-correlation function as w. This 
result is often known as Campbell’s theorem. 

The situation may be pictured by noting that the observation f##(x,, y,) 
contains contributions from all the f(x,y) in the neighbourhood of (x,,¥,), 
each with weight w(é,7) depending on the distance (€,7) between (x,y) and 
the point of observation (x,,y,). Evidently the weighting of the principal 
contributions is not much changed if the point of observation is moved to (x., y2) 
such that w(x,—%*2,¥,;—Y¥2)~w(o,o). Consequently the observed values few 
at (x,,y,) and (x9, y_) will be correlated in these circumstances whether or not 
f(x,y) is itself serially correlated. 

An example of this effect may be found in the observations of the variable 
component of radial velocity across the solar disk reported by Hart (1956). 
The auto-correlation function of this component falls to one-half for a shift of 
about +0°-5 heliographic longitude, corresponding on average to +4"°5 are. 
This correlation-distance agrees with the spread due to seeing, which was 
estimated as o-2 on the Oxford scale, corresponding to a tremor disk about 5” 
in radius (Plaskett 1952). That the observed correlation is due to seeing, rather 
than to any structure on the Sun itself, is the natural inference from Miss Hart’s 
Figs. 1 and 2, which show the variations in velocity plotted against angle, and 
her Fig. 4 which shows them plotted against distance along the solar surface. 
The angular plot has the appearance of a statistically stationary series, whereas 
in the distance plot the fluctuations appear to be much sharper at the end where 
the solar surface is seen more normally than at the other end where it is seen more 
foreshortened. Of course this conclusion is not in any way inconsistent with 
genuine velocity fluctuations having been observed, nor does it call into question 
Miss Hart’s discovery that the larger features of the velocity field persist for 
longer times than had been previously realized. 

2.4. Peaks in auto-correlation functions caused by the finite size of the available 
sample.—Suppose that an observed quantity, such as the deviation of velocity 
or brightness from its mean value, takes the -value f(x,y) at a point (x,y) on 
the solar disk, and f(x+é,y+) at a point distant (¢,7) from (x,y). Then 
as (x,y) explores the region covered by the observations, (£,y) remaining 
constant, the product f(x, y).f(x+€,y+7) takes many values depending on 
the details of the particular field studied. These variations are of little interest 
in themselves, and in taking the normalized auto-correlation function for shift 
(é,n) we form the normalized mean of these products. This mean may be 
written 


Cé,n)= Si (2.41) 
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where the summations are over all suitable values of f(x,y) available, say N in 
number. The function C defined in this way has useful properties with respect 
to the values of f(x,y) on which it is based; for example, its Fourier transform 
is the normalized squared modulus of the Fourier transform of these values 
(see equation (2.55)). But if another set of values of f(x,y) were taken, the C 
derived from it would generally differ from the previous one, because of failure 
of the variations in the products f(x,y).f(x+é,y+7) to ‘‘average out”’ 
completely. These differences are certainly properties of the sets of data in 
question but, from the point of view of characterizing granulation, they are 
accidental and are no more significant than are the variations in the individual 
products themselves. 

For (€,7) such that C(£,7) would be practically zero when averaged over 
a very large set of data, the accidental correlation over a sample containing n 
independent values has by equation (2.41) a root-mean-square (r.m.s.) 
expectation n—?, 

Similarly, if N values of f(x,y) are measured, each with independent r.m.s. 
error ¢, the error in the derived C is uncorrelated from point to point and has 
r.m.s. value «eN~"2(2f2 + «2)¥2/( f2 + €2), 
C(o,0) by e/( f?+€). 

In practice, errors such as those made in setting a cross-wire, statistical error 
due to photographic grain, and rounding errors in calculation, are usually almost 
independent for each measured value. Measurement is often made at a 
sufficiently close interval in (x,y) that the values of f(x,y) are fairly strongly 
correlated over a range of several times the measuring interval; say R times 
this interval. C(o,o) then stands out from the smooth run of C by an amount 
«*{e2 + f2}-1 and the error of measurement € can be estimated (see Section 3 for 
an example of this). The scatter of r.m.s. amount «N~?(2f2 + %)¥2/(f? +?) 
in the other values of C is not usually important if N is reasonably large. 

These results are obvious enough theoretically, but unless they are taken 
into account the appearance of auto-correlation curves derived from observations 
may be, at first sight, deceptive. Such curves often have peaks which are much 
larger than is reasonably consistent with the effect of observational error 
eN~'?(2f? + €*)"?/(f2+%) or with the expected fluctuation N-? in the mean 
of N independent products. Moreover, these peaks are not represented by single 
isolated points, but are apparently evidenced by the consistent trend of a number 
of successive points, and so may seem “‘convincing’’. 

All these putative criteria are in fact incorrect. We have seen that the 
statistical significance of a feature is properly judged with respect to an r.m.s. 
expectation m~"?, The r.m.s. amplitude N-"?, which would be expected if all 
the N observed values were independent, may be considerably smaller than n-*?, 
since in two dimensions x ~ N/R®, and in one dimension (f(x, y) replaced by f(x)) 
n= N/R (the definition of R, which was left loose above, can be chosen so as to 
make these approximations into equalities). The effect of the observational 
error € is usually smaller still. Moreover, the consistent trend of the points 
constituting a feature is merely a consequence of the fact that if f(x, y) is smooth 
then, by equation (2.41), C(€, 7) will also be smooth and any feature in it necessarily 
involves the smooth run of some R neighbouring calculated values. 

A good example of these effects has been presented by Zweig (1956), who 
shows in his Fig. 2 how a secondary peak in the auto-correlation of a sample 


The error of measurement augments 
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of photographic grain disappears as more data are added; this peak would have 
been taken as significant if the measured values on which it was based had been 
treated as mutually independent. Wlerick (1957) has demonstrated that peaks, 
other than the one at the origin, do not generally repeat in auto-correlation curves 
derived from different samples of actual granulation data. 

Fig. 2 of the present paper shows the result of a numerical experiment 
performed on the EDSACII computer to illustrate these relationships. 
A random series of numbers was convoluted (equation (2.31)) in turn with 
three functions w such that long runs of the convolutions f*w would have the 
auto-correlation functions displayed in the left-hand column of the figure. The 
respective auto-correlation functions calculated from a run of N= 100 values 
of f*w are shown in the right-hand column. The convoluting functions w used 
were derived in the manner described in Section 2.5. The random numbers 
were from a set of some 12000 generated on the EDSAC I computer and 
available on punched tape in the EDSAC library. They have a Gaussian 
probability distribution with zero mean and variance 100, and their properties 
have been subjected to statistical tests. The values of N/R are about 30, 25 and 
50 respectively for curves (6), (d) and (f). The r.m.s, amplitudes, in the region 
away from the peak at the origin are respectively 0-19, 0°26 and 0-16. (These 
values are based on many more ordinates than it was convenient to plot in the 
figure.) These auto-correlation curves were created by a numerical process the 
properties of which can be determined, and so there is no room for doubt that 
those features of the right-hand set of curves which fail to show up in the left-hand 
set are an effect of the finite size of sample, and do not represent statistical 
properties in the usual sense. 

The auto-correlation of the brightness distribution in solar granulation given 
as Fig. 3 by Stuart and Rush (1954) shows similar effects. It fluctuates beyond 
its first zero with an r.m.s. amplitude of about 0-2. As N is 96 and n is probably 
about half this value, the expected amplitude is about 0-15, and it is unsafe to 
attach any significance to the peaks occurring in this region. Similarly Fig. 5 
of Miss Hart’s (1956) paper shows an auto-correlation of solar velocities which 
has an r.m.s. amplitude beyond the first zero of 0-06; N= 140 and n is 30-40. 
The expected amplitude is therefore about 0-17, so the peaks beyond the first 
zero are actually fortuitously small. 

2.5. Artificial periodicities caused by smoothing and subtraction of running 
means.—I\n this section we consider for shortness one-dimensional functions, 
but the results are similar in many dimensions and if desired the variables can 
be regarded as multi-dimensional vectors. Lower- and upper-case letters will 
be used to denote functions which are Fourier pairs, and the Fourier operators F 
and F* will have the meanings implied by 


F[v(x)]=V(t)= ts o(x) exp {amixt} dx, (2.51) 


F*(V(t)]= | a V(t) exp {— 2mixt} dt (2.52) 


= v(x). (2.53) 


Equation (2.53) is the Fourier inversion formula (Titchmarsh 1932), which may 
be expressed by the operator equation FF*=F*F=12. We shall need also the 
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Fic. 2.—The effect of size sample on secondary peaks of derived auto-correlation functions. 
(a), (c) and (e). Auto-correlation curves corresponding to very long runs of three statistically 
stationary series. (b), (d) and(f). The corresponding auto-correlation curves calculated from 
samples, each comprising 100 values, of the series to which (a), (c\ and (e) refer. The right hand 
set of curves, derived from the limited samples, show secondary peaks which are not present in the 
left hand set and which therefore have no statistical significance. The heights of these peaks are 
discussed quantitatively in the text. 
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Fourier product theorem (ibid.) 


F[u*o] = Flu]. Flv] (2.54) 
where the asterisk has the meaning given in equation (2.31). In the special case 
v(x)=u(—x), (2.54) gives 

F[C(u]] = [Fu]? 


C[u) = F*[|Ul*) (2.55) 


where C[u] is the (unnormalized) auto-correlation function defined by 


C[u}= i) * w(é).u(é—x) dé (2.56) 


and |UJ* is called the power-spectrum of wu. 

Suppose that interest centres in the sharper features of an observed quantity 
f(x) rather than in its slower variations. The attempt may be made to remove 
the slow changes by subtracting a running mean; that is to say, replacing f(x) 


Ki g(x) =f(x) —F(e) (2.57) 
where I +a 
flay= 5 [" fe—Oaé, 


I 


=1[” f(x—&)rect (é/a) dé, 


a 
= ~ fe rect (€/a) ; (2.58) 


‘rect’? has the same meaning as in Section 2.2. Then by equation 


F[/]=F[/]. sinc (at), 
F[g] =F]. {1 —sine (at)}. . (2.59) 

These relations are illustrated in Figs. 3-5. Fig. 3 shows (a) an assumed real 
spatial spectrum F[f]; (6) the spectrum sinc(at) for a chosen value of a; 
(c) 1—sine (at); and (d) the result of multiplying the original spectrum by 
1—sinc(at) in accordance with equation (2.59). It will be seen that the 
resultant spectrum F[g] has a fairly pronounced maximum for a certain value 
of the frequency argument t. This is equivalent to saying that the function g 
calculated by equation (2.57) has a quasi-periodic structure which is a property 
of the method of calculation. If f represented (for example) a brightness 
distribution, and if the modified function g were convoluted with radial velocity 
observations (say) from which a similar running mean had also been subtracted, 
then according to equation (2.54) the spectrum of the convolution would be 
even narrower than that of g and the corresponding velocity function, and the 
artificial periodicity would be correspondingly enhanced. 

The assumed real spectrum F[f] shown in Fig. 3 (a) refers, as it stands, to 
the observation of a symmetrical and fairly compact solar feature. A more 
realistic illustration of effects occurring in granulation measurements is obtained 
when curve (a) of Fig. 3 is re-interpreted as the square-root of the mean power 
spectrum of the observed quantity. Since eqn. (2.59) is true separately for 
each feature included in the observations, it is true also in the r.m.s. sense. 
Accordingly the effect of subtracting a running mean from a stochastic function 
is represented by Fig. 3, and the conclusions drawn above from this figure remain 
true, in the new interpretation, 
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(2.54), 
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Fic. 3.—The distortion of frequency-distribution by the subtraction ef running means. 
(a) Assumed amplitude spectrum (Fourier transform) of an observable distribution, such as bright- 
ness or velocity along a line of the Sun. (6) Mean spectrum of the running mean of a function 
having a uniform spectral density (see text). (c) Mean spectrum which results from subtracting a 
running mean from a function having a uniform spectral density. (d) The spectrum which results 
from subtracting a running mean from a function having the spectrum shown in (a). 

The subtraction of the running mean has produced a peaked spectrum, and the observed function 
after subtraction of the running mean will accordingly show apparent periodicities which are artificial, 





No. 5, 1959 On the interpretation of solar granulation 487 


In a numerical experiment to illustrate these effects, carried out on the 
EDSAC II computer, a sample out of a set of random numbers (see Section 2.4) 
was convoluted in turn with three weighting functions w according to the numerical 
equivalent of eqn. (2.31). Let the sample used be the numerical representation 
of a function u(x). Then by eqns. (2.54), (2.55) and (2.53) 


F[uew])=U.W 
F[C[uew)}]=|U. WP 
C[uew) = F*(|UP. |W). (2.591) 
If the spectrum U is uniform, it follows by a second application of eqn. (2.55) 
that €[uew]=€[w] which justifies the statements made earlier to this effect. 
It also follows (Woodward 1953) that functions which are random (that is, 
entropy maximized), subject to having mean-squared spectra proportional 


respectively to the squared ordinates of curves (a) and (d) of Fig. 3, can be 
generated by convoluting random numbers with functions w, and w, given by 


F[w,(x)] =the ordinates of curve (a) Fig. 3, 
F[w,(x)] =the ordinates of curve (d) Fig. 3. 


Fig. 4 shows (a) the result of this operation with w,; (6) with w,; and (c) with 
a function w,(x)=w,(2x); the scale of seconds of arc in the figure is of course 
purely notional. The curves displayed on the left-hand side of Fig. 2 are in fact 
the auto-correlation functions of w,, w,, w3, and the curves on the right-hand 
side of this figure are the respective auto-correlation functions of curves (a)-—(c) 
Fig. 4. 

The process of generating curve (a) Fig. 4 is analogous to the blurring of a 
random distribution when it is observed with resolution corresponding to a 
‘spread’ or ‘‘ apparatus ’’ function w,, and it will be seen that some appearance 
of periodicity has already been produced by this blurring (cf. Section 2.2). 
The appearance of periodicity is enhanced in curve (6), which is of course simply 
curve (a) with a running mean subtracted; this is especially noticeable at the 
left-hand end. In curve (6) the maxima (except those corresponding to 
fluctuations of less than 1/10 of the general amplitude) always occur at positive 
values of the function, and the minima correspondingly at negative values; whereas 
this is not everywhere true of curve (a). The effect of convoluting two portions 
of curve (0) is shown in Fig. 5. The periodicity has here become very obvious 
and, in the analysis of actual observations, this might well be taken as physically 
significant if it were not realized that its appearance is a simple property of the 
numerical process used. 

The putative periodicity in the velocity-brightness cross-correlation function 
shown in Fig. 1 of the paper of Stuart and Rush (1954) is apparently due to 
effects of this kind. These authors are aware of the falsification which occurs 
when a running mean is subtracted, and seek to avoid it by varying the interval a 
in accordance with the trend of the observed values; a better method will be 
described in Section 3. These considerations do not of course invalidate Stuart 
and Rush’s conclusion that the velocity-brightness correlation is greater in the 
higher spatial frequencies. 

We saw in Section 2.4 that the auto-correlation curves shown on the right 
of Fig. 2 are subject to statistical errors due to the limited size of sample on which 

















Fic. 4.—Apparent structures in smoothed random distributions. 


(a) Plot of a sample of random numbers after a smoothing process (convolution) has been applied. 
(b) Identical to (a) except that the function used in the convolution is such that a running mean is 
subtracted as part of the smoothing operation. (c) The convoluting function used here has the 
same shape as that in (a) but only half the width. The curve is accordingly similar to (a) except 
that it has twice the *‘ resolution’’. 


The smoothing process has caused some appearance of quasi-periodicity in curve (a), and this has 
been enhanced by the subtraction of running mean in (b). All the curves were made from the same 
sample of random numbers. The features marked A and B are analogous to features of granulation 
photographs to which significance has sometimes been attached. The changes in their appearance 
with variation of resolution and numerical processing is discussed in the text. 
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they are based, and are significant only in so far as they agree with the corresponding 
curves on the left of the figure. We now see that the negative peak in curve (c) 
is an artifact caused by the subtraction of running means, and the width of the 
peak at the origin in each curve is likewise simply a property of the blurring process, 
as is most obviously demonstrated by curve (¢) which (from the definition of w3) 
corresponds to twice the resolution of the other two. Indeed, these three curves 
cannot (by the very way they were derived) tell us anything about the random 
numbers which, in these experiments, played a part analogous to the object- 
distribution in real observation. The conclusion is that none of the curves of 
Fig. 2 gives any explicit information about the fact that these numbers had zero 
mean serial correlation. Analogous curves produced from granulation observa- 
tions may similarly lack significance. A comparison of curves (a) and (0) does, 
of course, give information about the serial correlation of the numbers used in 
the experiment. This emphasizes the importance of including in granulation 
observations an accurate determination of the blurring function corresponding 
to curve (a); clearly, this determination must be independent of any assumptions 
about the phenomenon under investigation. 

If the curves in Fig. 4 represented scans across images of granulation, then 
in a high contrast reproduction of the image corresponding to curve (@) it might 


PeMWrn Aeros I dX Af. 
i TAAANS 


Fic. 5.—Artificial periodicities created by numerical processes. The curve shows the cross- 
correlation function between two parts of the function of which Fig. 4(b) is a sample. It is analogous 
to cross-correlation functions between solar velocity and brightness observations from both of which 


running means have been subtracted. The marked quasi-periodic appearance is totally artificial; 
the curve has been derived from a random set of numbers. 





appear that the feature matked B had a size of 0”*5. Significance has often:been 
attached to the size of the smallest features apparent in granulation images, as 
indicating either the resolution achieved, or the size of the granulation structure. 
We see in curve (6), which has the same resolution as (a) but from which a 
running mean has been subtracted, that the size of feature B has been reduced to 
about 0-3 and A is nearly as small, although it appeared much broader than B in 
curve (a). Curve (c), which corresponds to twice the resolution, shows A as 
a complex structure extending over more than 1”, and B is no longer identifiable. 
Moreover, both A and B are in fact spurious in the sense that complete 
resolution would reveal randomness devoid of any coherent structure. 
Caution is evidently required in interpreting granulation observations in this 
manner. 

2.6. The indirect equivalence of one- and two-dimensional spatial spectra.—The 
observed brightness and velocity fields on the solar disk are two-dimensional 
manifestations of effects which vary with depth, as well as perpendicular to the 
line of sight. ‘The relating of these observations to the three-dimensional 
velocity- and radiation-fields in the Sun would need detailed optical and 
hydrodynamic treatment, and so lies outside the scope of the present paper. 


35 
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Our concern here is with the relationship between the two-dimensional variation 
of a quantity over the solar disk and the information obtained by focusing the 
Sun on to a spectrograph slit, or from scanning a solar image in a microphotometer. 

Some general discussion of scans of two-dimensional distributions has already 
been given in connection with analysis of photographic grain (Fellgett 1953)*, 
and the discussion which follows will be mainly concerned with the special case 
in which the observation can be represented as f(x) = f(x, 7), where » is a constant. 
That is to say, f, represents the observation secured when an image of the 
distribution f(x, y) is viewed through a long straight slit of negligible width lying 
along the line y=, or when a microphotometer spot of negligible size scans 
along y= 7. 

Fig. 6 shows a Fourier component of f(x, y) having frequency t, the wavefronts 
of which make an angle @ with the line of the one-dimensional scan. This 
component is seen to give rise to a frequency ¢,=tsin@ in the one-dimensional 
spectrum of f,(x). As t,<t, the total amplitude of the Fourier component of f, 
having frequency ¢, is the sum of contributions from all the two-dimensional 
components of f(x, y) which have frequencies t>t,. Accordingly, the determina- 
tion of the two-dimensional spectrum F[f(x,y)] from the one-dimensional 
F[/,(x)] involves the solution of an integral equation, one form of which is given 
by Uberoi (1955)+ as his eqn. (6). Uberoi has emphasized that the two- 
dimensional sinusoidal contrast transmission factors (Fellgett and Linfoot 1955), 
which are used to characterize the resolution of an optical system, are of course not 
directly applicable to F[ f,(x)]. But even if the mistake is avoided of interpreting 
the one-dimensional spectrum of a scan as if it were a section of the two-dimensional 
spectrum, the difficulty remains that experimental or sampling error may so 
dominate the solution of the integral equation as to prevent a usefully accurate 
estimate of F[f(x,y)] being made from the one-dimensional observations f,, 
especially in the lower-frequency range. Dr E. H. Linfoot has pointed out 
to me that this difficulty is avoided if the small scanning spot, commonly used 
in this kind of measurement, is replaced by a long slit perpendicular to the 
direction of scan, as employed in measuring optical spectra. With this method 
of scanning (which is of course quite different from the case of a fixed spectrograph 
slit masking the image) the observed quantity may be written 


g(x) = | dy f(x,y) (2.61) 


where the integration is over the length of the slit. 

In general, when f(x,y) is scanned with a mask having transmission 
w(x —a,y—b), where (a,b) is at any instant the position of a chosen reference 
point of the mask, the observed quantity is 


$(a,b)= | | c fx, y)wo(x—a,y—b) dedy. (2.62) 


Equation (2.54) then shows how the spectrum of f has been modified by the 
scanning process, and the effect of any practical shape of slit or other aperture 
can be determined in this way. 


* The Fourier operator in this paper is incorrect by a factor 27 in the exponent. 
+ In equations 16 to 19 of Uberoi’s paper, and in the next unnumbered expression, the factor 
printed 


~=( 2-8) -G)')" seen eS 
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The function f,(*) and g(x) have particularly simple properties which it is 
more convenient to derive by a method which is special to these limiting cases. 


vy 
wave LENSTe YE tin @ 


Fic. 6.—Projection of a two-dimensional Fourier component onto a line scan (see text). 


The effects of the finite length and width of the slits used in practice can then 
be assessed by reference back to eqn. (2.62). Consider the spectrum of f(x,y), 


namely 


Film y=Ft)=[[" ferexptanitas+yt)}dedy. (2.63) 


Along the s-axis t=o0, 


F(o)=["_ dvexp{anixy}[" dyf(x.y) (2.64) 
and, if the scanning slit is so long that infinite limits can be used in equation 


(2.61), 
F(s,0)=F[g(x)}. (2.65) 
It follows that F[f(x,y)] can be simply determined from scans of a long narrow 
slit across f(x,y) in different directions, without the need to solve an integral 
equation. If f is symmetric so that F(s,t)=F(r), r?=s*+#, then only a single 
scan is needed. These favourable properties may be thought of as arising because 
any localized feature in f(x,y) remains well localized with respect to the set of g 
taken over all directions, so that ‘‘sharpness’’ is preserved; and because the 
whole of f(x,y) is covered in each scan. 
The whole set of f,(x), for all n, also has these properties of localization and 
completeness, and F[f(x,y)] can be simply derived from the set. Numerical 


two-dimensional Fourier transformation is commonly done in this way by using 
equation (2.63) in the form 


Fif(x,y)]= i) dn exp {2mint} dx f(x,n) exp {2mixs} (2.66) 


35* 





492 Peter Fellgett Vol. 119 


where the integrations are nominally infinite and in practice are taken over the 
available sample of f(x,y). ‘This equation may be written in the more compact 


form 

F(s,t)=F(FLf,)) (2.67) 
where the inner Fourier operation is with respect to the dummy variable x, 
and the outer one with respect to 7. When, however, only a single f, is known 
it is evidently impossible to derive F(s, t) explicitly unless further information 
is available which in effect enables the whole of f(x,y) to be inferred from the 


sample f,. Integrating both sides of equations (2.66) or (2.67) yields (after a 
limiting process) 


| dt F(s,t)=F[f,],  =0. (2.68) 


This equation is symmetrical with equation (2.64). In many applications where a 
single scan f, is available, it is permissible to regard this scan as being along the 
x-axis. The equation shows that knowledge of this f, alone fixes only the t-integral 
of the two-dimensional spectrum F(s,t); this function may be pictured as the 
result of projecting or ‘‘squashing flat’’ F(s,t) on to the t-axis. If f is 
symmetrical, equation (2.68) takes the form 
rdr 

Finl=2[ Panay 1-2. (2.69) 
This is a special case of the conclusion drawn earlier, from Fig. 6, that the value 
taken by the one-dimensional spectral function F[f,] at a given frequency is the 
weighted sum of the values taken by the two-dimensional spectrum F(s,?) at 
equal and higher frequencies. When (2.69) applies, F(s,t) can be explicitly 
evaluated from f, but (as has been mentioned above) often only with low 
accuracy. 

Indeed if f is random, subject to having a statistical mean power-spectrum, 
the amplitude of a Fourier component of f, or of any scan of it, has a Gaussian 
probability distribution with zero mean and variance equal to the mean power 
(Shannon 1948). The squared amplitude is therefore an estimate of this power 
having variance equal to twice its mean value. It follows that the spectral density 
of a stochastic phenomenon can be estimated accurately only by considerable 
averaging, either with respect to frequency or over different scans. This is so 
even when errors of measurement are negligible; and conversely errors of 
measurement only a little smaller than the quantities measured need not greatly 
increase the variance of the estimates of the principal components. For point 
scans, the relative uncertainty can be so much increased in the transformation 
from the spectrum of f,(x) to that of f(x, y) (using the inverse of equation (2.69)) 
that it is not usually worth while to attempt this transformation unless very many 
scans f, are measured. It may then be better to arrange them in a raster so that 
the two-dimensional spectrum can be obtained directly using equation (2.66). 

Slit scans, as well as being more convenient (cquation (2.65)), avoid this serious 
loss of accuracy; g(x) is a much better sample of a stochastic function f than 
is f,(x). Lengthening the slit extends to lower frequencies the frequency-range 
over which equation (2.65) is a good approximation to equation (2.62). It also has 
the practical advantage of increasing the total fluctuating signal, asymptotically as 
the square root of the slit-length. This enables the amplifier gain to be reduced, 
so that amplifier noise and detector dark-noise have less effect. The signal-to-noise 
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ratio in the light itself is not increased, however, since noise uncorrelated at the 
two ends of the slit increases in the same manner as the signal. Photon noise, 
photographic grain, and the scintillation of the granulation behave in the same 
way. The mean brightness increases directly as the slit-length, so that if for 
example the r.m.s. brightness fluctuation were 0-05 of the mean for 1” slit length, 
it would be 0-002 for a slit 600” long. The mean brightness does not affect the 
amplifier, since zero-frequency response is not required and the detector will 
normally be coupled to the amplifier through a condenser capable of passing 
alternating signals of order 10~* V while blocking a steady potential difference 
some 10° times larger. (The situation would be quite different if it were necessary 
to measure the mean brightness with a relative accuracy of 10-*, which would 
require the gain-stability and dynamic range of the amplifier to be enhanced by 
negative feedback.) The low contrast does imply that direct slit-scans of the Sun 
will be sensitive to changes in the overall brightness of the sunlight, which will 
presumably be due more to transparency than to seeing effects for a disk 30’ in 
diameter. If these changes cannot be sufficiently separated in frequency from 
the desired ones by a suitable choice of rate of scan, it may be necessary to use 
balanced detectors; for example the difference in the light passed by two halves 
of a slit may be measured. A brightness discrimination of 10~* has proved 
practicable in the balanced photocells often used to measure small optical 
deflections, even where zero-frequency response is required. Photon noise, as 
modified by the detector, is especially relevant to the measurement of velocity. 
If the r.m.s. velocity variation were say 0-5 km/sec for a 1” slit, it would be 
0-02 km/sec r.m.s. for a 600” slit. I am grateful to Mr D. W. Beggs for information 
that the Cambridge Solar Magnetograph gives an r.m.s. noise equivalent to 
0-01 km/sec for a slit only 60” long, indicating that the 12-inch aperture of this 
instrument is sufficient for the measurement of Doppler-granulation. The 
ratio of this r.m.s. noise to the velocity amplitude accords as to order of magnitude 
with the accuracy of measurement attained in the point-scans discussed in 
Section 3, as the present discussion leads us to expect. 

The disadvantages of the point-scan f, are somewhat mitigated when the 
purpose is to determine the two-dimensional power spectrum of a function 
f(x,y) which is statistically stationary as well as isotropic. In photographic 
measurements, it has the advantage that the non-linear relation of density to 
exposure is more easily allowed for than in other methods of scanning. The 
auto-correlation functions of f(x,y) and f,(x) may be written 


C(é,n)y=CLf= i) de dy f(x,y) .f(x+éy +n), 


CO), = CL = [ defen) fle +67) 
The domain of integration is as in equation (2.66). For 7=o0 


C(é,0),= dy { dx f(x,y). fle+é.y), 
- dy C(E)y, 


oe C(é), dy, 
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where the last step follows from the stationary property. Hence knowledge 
of f, is sufficient to fix C(f,0). By isotropy C(é,n)= C(C), €+7?=@, and by 
equation (2.55) the required power spectrum is |F(s,t)/>=F[C(é,»),]. In the 
practical use of this method tke heavy work of two-dimensional Fourier 
transformation is avoided by forming 


Feoyh=F| (CEn)ean | (2.692) 


-"[2[° CO ata || (2.693) 


The proof of these relations follows by comparison with equation (2.64). 


3. The data and calculations 

The foregoing discussion illustrates that a central problem in the interpretation 
of solar granulation is to determine, in a way which does not introduce artifacts, 
which frequencies are most strongly present in the observed velocity and 
brightness fields, and which frequencies contribute most strongly to the 
correlation between velocity and brightness. It is first necessary to give a 
definite meaning to these questions in a way which is acceptable to intuitive 
notions. 

Let v(x) and b(x) be respectively the deviations from their mean values of 
the Doppler velocity, and of the brightness, at a point x of the Sun; as before, 
x may be regarded as either one- or multi-dimensional, and integrations will be 
over the available data unless otherwise indicated. The total power, or un- 
normalized variance, of v is defined (since d=0) as 


C.= | Peas; (3.1) 


the definition for } is of course similar, here and in what follows. The correlation 
between 5 and v is defined as 

Co= fee) b(x) dx. (3.2) 
The ordinary correlation coefficient is obtained when C,,, is normalized by division 
by (C,C,)'*. We have seen that difficulties may arise when the attempt is made 
to modify the observed quantities v(x) or b(x) by subtracting running means, or 
in other ways, for the purpose of finding the contribution to C, or C, from 
different frequency ranges. The problem becomes tractable when account is 


taken of the fact that C,, C,, are the respective values for =o of the auto- 
correlation function 


Cle] =C(é)_= | v(x) oe +8) as (3.3) 
and of the cross-correlation function 
Cle, b]=C(eo= J o(x) d(+ €) dx. (3.4) 
By equations (2.51)-(2.55) 
F[C(E)o, 0] = V(t). B*(¢), (3-5) 
C(E)e,o=F*[V(t). B*(t)], (3-6) 


Cor=COno= | _ VIB(e)de (3-7) 


a VBedt+ |" VBeat-+ [~ VB de, (3.8) 
ty ty 


-—-@® 
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where the asterisk denotes the complex conjugate. Equation (3.8) expresses 
the correlation C,,, as a sum of contributions from different frequency ranges. 
Each of these contributions is the integral, over the appropriate range of 
frequency t, of the cross-correlation spectrum F[C(€),»}. 

Accordingly, we may without inconsistency define 


{ V(t)B*(t) dt 
ty 


as the contribution to the correlation C,,, from the frequency components of the 
cross-correlation spectrum lying in the interval t, to t,. The formal substitution 
v(x) = (x) shows that in the same manner 


f \V(t)P at 
4 


may be called the contribution to the power of v from the frequency components 
of the power-spectrum of v lying in the interval t, to t,. It will be seen that, in 
agreement with intuitive ideas, the contribution to C,, from any frequency 
interval is large and positive if V and B, the frequency-spectra of v and 3}, 
both take large values and have the same phase in this interval; it is large and 
negative if V and B are large but of opposite phase. In a similar way, the 
contribution to C, or C, is large for those frequency intervals in which the 
power-spectra |V[* or |Bf* take large values. If C,, C, or C,, is real, the 
contribution from —t, to —f, is the complex conjugate of that from #, to ty. 
The contribution from the interval t, to t, >t, will then be understood to mean 
the contribution from the double range t, < |t|<#, (a similar convention exists 
in electrical engineering). 

These definitions, by giving a precise form to questions concerning which 
frequencies contribute most to C,, C, or C,,, make it possible to calculate these 
contributions from the experimental data. It was shown in Section 2 that any 
‘‘smoothing’’ or ‘‘unsmoothing’’ process distorts the frequency spectra 
(equation (2.54)). Therefore the only safe procedure is to calculate the frequency 
spectra numerically from unmodified raw data, and to pick out explicitly any 
frequency ranges that show themselves to be of special interest. 

The data used are the 340 brightness and velocity measures obtained by 
Richardson and Schwarzschild (1950) at Mt Wilson and the three sets of 116 
measures obtained by Plaskett (1954) at Oxford on 1949 September 11-4306, 
1949 September 11-4319, and 1950 May 11-7277. ‘These data will be referred 
to as Mt.Wilson, Oxford September A, Oxford September B and Oxford May 
respectively. The auto-correlation functions and power spectra of the Mt Wilson 
data extend slightly the calculations already published and discussed by Frenkiel 
and Schwarzschild (1952), and by Uberoi (1955). The calculations on the Oxford 
data, and the calculation of the cross-correlation function and its spectrum for 
both the Oxford and the Mt Wilson data, are new. The data represent one- 
dimensional scans; when the one-dimensional correlation and Fourier calcula- 
tions on them were completed it appeared that little additional information would 
be gained by deriving spatial spectra of the two-dimensional granulation on the 
assumption of isotropy (Section 2.6). A much larger sample would be needed 
before this could be done with useful accuracy. 
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The auto- and cross-correlations were programmed and run on the EDSACI 
computer by Mrs M. Mutch. The Fourier transforms were calculated by 
Miss C. P. Wilks on Lipson Beevers strips (Lipson and Beevers 1936). This 
combination of machine and semi-hand methods was very economical in relation 
to the sources available when the calculations were made. It has been rendered 
obsolete, however, by the much faster and more convenient computers which 
have since become available. A programme which I wrote for the EDSAC II 
computer enables the Fourier transform of 1000-ordinate data to be calculated 
at a rate of 20 output ordinates per minute, and further economies would be 
possible by improved programming. When new data on solar granulation are 
published it should be possible to do the calculations on them with little more 
time and labour than is needed to convert them to digital form. 

The correlations were evaluated in the numerical form 


N-M-1 


a,b, mm 


C= r= M 
a,, 





ee -* V-M-1 
r= M 


where a,(a,)ay_, are the N ordinates of one set of data at equal intervals of 
the argument, 5)(b,)by_, the corresponding ordinates of the other set; m is 
the value of the shift and is calculated in the range -MemeM. If a,=6, 
then ¢,, is an auto-correlate. This form of c,, was suggested to me by 
Dr M. V. Wilkes. For use in connection with Fourier methods, it has the 
advantages that the normalizing factor and the statistical weight are both 
independent of m. The disadvantage is that the weight is less than it could be 
for small values of the shift m. The auto-correlation in this form is statistically 
symmetrical, but not necessarily symmetrical for any particular sample of data. 

The results of the computations are presented graphically in Figs. 7-11. 
The Mt Wilson data are at intervals of 0”-45 arc= 328 km on the Sun, and the 
correlations are calculated up to a maximum shift of M = 60 of these intervals. 
The Oxford data are at intervals of 1”-04 arc, corresponding to 759 and 756 km 
on the Sun for 1949 September and 1950 May respectively, and the correlations 
are for Me 30. The total intervals covered by the Mt Wilson and Oxford data 
are thus 111500km and 88000km, and the maximum shifts correspond to 
+19700 km and +22700 km respectively. Only the cosine components of the 
Fourier transforms are shown; as mentioned above, the sine components are 
statistically zero. 


4. Discussion 


4.1. The brightness auto-correlation and power spectra.—T he present brightness 
auto-correlation curve derived from the Mt Wilson data, the corresponding curve 
of Frenkiel and Schwarzschild (1952, Fig. 1), and that of Uberoi (1955 a, Fig. 3), 
differ among themselves by more than might be expected from the differences 
in numerical techniques. 

Each of the four brightness auto-correlation curves shown in Figs. 7-10 
has a peak centred on zero shift with width roughly corresponding to the expected 
instrumental resolution. For the Oxford results, this peak is superposed on a 
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wider one (particularly marked for September A) which may be of instrumental 
origin associated with slow changes in sensitivity or black-level. 

We saw in Section 2.6 that an ordinate of the unsmoothed power spectrum 
of a ‘‘ random ”’ function has an r.m.s. uncertainty equal to +/2 times the expected 
value of the ordinate. For the Oxford data, the power spectra have been found 
by Fourier transformation of correlation functions determined over a total range 
of 60 units of shift. The total length of the available data is 116 units, and 
accordingly the power spectra found have a resolution of only about one-half 
that which corresponds to the whole run of the data. This implicit smoothing 
gives a ratio of signal to r.m.s. noise of about unity for each ordinate. For the 
Mt Wilson data, the resolution is about one-third the value corresponding to 
the length of the data, and the ratio of signal to r.m.s. noise is about 1-2 for each 
ordinate. The o® ordinates of the brightness-spectra have no physical 
significance; they simply reflect the mean level of the brightness compared 
with the reference level chosen. When these facts are taken into account, the 
brightness power-spectra appear to be more or less steadily falling functions of 
frequency and present no very remarkable features. The independent brightness- 
correlation curves of Uberoi (1955, Fig. 3) and of Frenkiel and Schwarzschild 
(1955, Fig. 1) also agree in showing no significant secondary peaks; the latter 
curve is based on about 2000 independent measures and hence has much greater 
weight than the brightness data used in the present work. 

4.2. The velocity auto-correlations and power spectra.—The o™ ordinate of 
each of the Oxford velocity auto-correlation curves is seen to be larger than would 
correspond to smooth interpolation between the neighbouring values. This 
reflects errors of measurement (‘‘noise’’) uncorrelated from point to point of 
the velocity data. It causes the mean of the ordinates of the power spectrum at 
large wave-numbers to be non-zero. The actual value of this mean is about 0-4, 
corresponding to a ratio (see Section 2.4) 


N = 0° ee 
San°°* W'S 


where S=the power of the noise-free data and N=the power of the noise 
fluctuations. This estimate agrees well with the ratio S/N=1-4 estimated by 
Plaskett (1954, p. 256). The Mt Wilson velocity spectrum has a mean value 
for large wave-numbers of about 0-05, corresponding to S/N=20; the value 
implied by Richardson and Schwarzschild is 25-4. 

Apart from the o ordinate, the Oxford May velocity auto-correlation 
appears similar to the Oxford brightness auto-correlations discussed in 
Section 4.1. The September A and September B correlations, however, have 
a markedly periodic appearance, which shows up as peakiness in the corre- 
sponding power spectra. The likelihood of these curves is quite small on the 
assumption that the periodicity is accidental. The fact that the periods differ, 
by a factor of nearly 2: 1, between the two exposures made within about a minute 
of time of each other is somewhat opposed to the view that the effect is instrumental. 
Its reality is also supported by the circumstance that similar periodicities show 
(though weakly) in the corresponding brightness auto-correlations, and 
that they are prominent and have consistent sign in the cross-correlation 
spectra. 
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Fic. 7.—Mt Wilson granulation observations. 
(b) Power spectra and cross-correlation spectra corresponding to the functions displayed in (a) (see text). 
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Fic. 8.—Oxford September A observations; plotted as in Fig. 7. 
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Fic. 8.—Oxford September A observations ; plotted as in Fig. J. 





Peter Feligett 


or 


OXFORD Ser8 


Mes 0 
J 





r 
za000Km 





20.000 


ee 
wD Ad 





oz 


See Lr [\ 
ia a 





= 





-0°6_) 


(a) 
Fic. 9.—Oxford September B observations; plotted as in Fig. 





No. 5, 1959 On the interpretation of solar granulation 


OXFORD Seer B 











ae 
L. T 


16.000 6ooo 4000 2000 





e 


VaevV 
ZL. vita pee 
1 j 
T _ T = 
Ku 


BaV 


1° 
i <—— 
T T 


> aa 





a ® 
L— 


L 





(6) 
Fic. 9.—Oxford September B observations ; plotted as in Fig. 7. 
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(a) 
10.—Oxford May observations; plotted as in Fig. 7. 
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. 10.—Oxford May observations; plotted as in Fig. 7. 
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Fic. 11.—“ Incorrect ’’ cross-correlations. 
(a) Cross-correlation between brightness from one Oxford spectrogram and grey + from a different 
one; top, May brightness—September A velocity; centre, September B brightness—September 
B velocity; bottom, September B brightness—May velocity. 
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Fic. 11.—* Incorrect”’ cross-correlations, 


(6) The cross-correlation spectra corresponding to the functions shown in (a). 
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4.3. The brightness—velocity cross-correlations and correlation spectra.—In the 
table which follows, the values calculated in the present work of the brightness- 
velocity cross-correlation coefficient for zero shifts are compared with the values 
derived by the observers themselves. 


Published value Present estimate 
Mt Wilson — 0°304 —0°344 
Oxford : September A —0o'26 —0°05 
September B — 0°29 — 0°42 
May +0°07 —o18 
For the Oxford data, the decrease in weight by a factor of about two at the origin 
(see Section 3) has altered the values found for the correlation coefficients so 
markedly as to suggest that they have little statistical significance. 

Plaskett states that the probability of getting a correlation arithmetically 
greater than 0-26 between two sets of 116 values taken from an uncorrelated 
parent population is less than 1 in 460. ‘This is not directly relevant, however, 
because the velocity and brightness data are each serially correlated, as their 
auto-correlation curves show (compare Section 2.4); and in fact 36 of the 
183 ordinates of the three Oxford cross-correlation curves, or about one-fifth 
of them, exceed 0-26. As a further experimental check on the probability of this 
value, three cross-correlations were computed (Fig. 11) between the brightness 
on one spectrogram and velocity on a different one. ‘The proportion of ordinates 
exceeding 0-26 is about one-quarter, and about one-sixth of the ordinates exceed 
0°30. This confirms that the values shown in the table could well have arisen 
by chance. 

Moreover, the ‘‘false’’ cross-correlation curves do not present any features 
that obviously distinguish them from the ‘‘correct’’ ones. Indeed one of the 
most ‘‘convincing’’ curves is September A brightness-September B velocity, 
which gives a correlation coefficient for zero shift of —o-41. Although in view 
of the persistence of velocity fields found by Hart (1956) this correlation may 
not be entirely accidental (but the change in the spectrum of the velocity field 
between the two exposures suggests that it is largely so) this circumstance 
nevertheless further weakens confidence in the reality of the putative correlations. 

Similarly the negative correlation of 0-344 found for the Mt Wilson cross- 
correlation curve at zero shift is arithmetically equalled or exceeded by 8 of the 
121 ordinates of this curve, and this value therefore appears to be too small to 
have much significance in itself. However, the negative peak at the origin does 
appear to be significantly sharper than the other peaks. Examination of the 
four spectra of ‘‘correct’’ cross-correlations in the light of this remark seems 
to lead to a fairly consistent picture, in the following way. 

The sign of the contribution to the cross-correlation from the lower frequencies 
in the correlation-spectra shows no apparent systematic trend, but the higher 
frequencies make a predominantly negative contribution for the Mt Wilson data, 
Oxford September A and Oxford September B; their average contribution is 
small for Oxford May. This is in agreement with the observation of Stuart 
and Rush (1954) that the negative brightness—velocity correlation for the 
Mt Wilson data can be approximately doubled by removing the lower frequencies. 
Moreover, in all the three cases where marked peaks in the velocity spectra for 
Oxford September A and B coincide with reasonably large values of the 
brightness spectrum, there is a negative peak at the same frequency in the 
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corresponding cross-correlation spectrum. The likelihood of these three 
contributions all being negative by chance is one-eighth. The ‘“‘false’’ 
cross-correlation spectra do not share this systematic property. 

The lack of consistent correlation at low frequencies need not be a property 
of the Sun, since (as was mentioned in Section 4.1) the observations at these 
frequencies may be strongly affected by slow drifts in instrumental characteristics. 
The differences between the correlation coefficients found by different observers 
may be partly due to this cause. 

If the spectrograph has an asymmetric spread function (sometimes called 
‘‘apparatus function’’), the non-linearity of the photographic response will 
evidently cause the centroid of the recorded lines to be shifted by an amount 
that depends on the intensity of the line and the local brightness. The effect, 
which is analogous to the so-called magnitude equation in astrometry, may cause 
an apparent correlation between brightness and velocity, independently of whether 
such correlation exists on the Sun, unless the solar and telluric comparison lines 
used have exactly equal strengths. This may be another source of error in the 
correlation coefficients which have been found. There is nothing to suggest 
positively that the Mt Wilson or Oxford data are affected in this way, but since 
even small brightness-correlated shifts of instrumental origin would seriously 
falsify the measured correlation it seems prudent to envisage that such effects 
may be present unless there is evidence to the contrary. A desirable check in 
any future observations would be to illuminate the spectrograph slit with a 
defocused image of the Sun, and to place in the optical train a screen (such as 
a grainy negative) arranged so that the brightness across the slit has granule-like 
variations which are known not to be associated with any real velocity variations. 


5. Conclusions 


There appears to be no reliable evidence of consistent peaks in the power 
spectra of either the observed brightness or velocity distributions. On the 
contrary, both kinds of spectrum fall on average with increasing frequency in a 
manner apparently dominated by instrumental resolution. The likelihood of 
the observed values of correlations between velocity and brightness occurring 
by chance is much greater than has been supposed. The evidence for real 
velocity—brightness correlations is weak for the Mt Wilson data, and very weak 
for the Oxford data, if it is based solely on the value of the correlation coefficient 
at zero shift. The spectra of the cross-correlation functions do, however, show 
consistently non-positive average contributions from the higher frequencies. 
For this reason, and on the ground of physical likelihood, it seems permissible, 
though by no means necessary, to assume that velocity and brightness are indeed 
negatively correlated. 

The strongly peaked velocity spectra shown by the two Oxford September 
exposures are rather surprising, but not inconsistent with the excellent 
Doppler-granulation pictures secured by McMath et al. (1956). Some of these 
pictures displayed at the I1.A.U. Dublin Congress suggest that the velocity 
variations over limited regions of the Sun may sometimes show strong sinusoidal 
components. The Oxford September A and B spectrograms were exposed 
within about a minute of time of each other and on nominally the same region 
of the Sun. If the effect is real and not merely accidental, one would like to know 
whether the changes between the two pictures represent time- or space-variations, 
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and whether the 2:1 frequency ratio represents a progressive change in 
frequency or a switch from one harmonic to another which might conceivably 
indicate magneto-hydrodynamic effects. 

Apart from this tentative suggestion, there seems no reason to disturb Uberoi’s 
conclusion that all the measures so far made can be satisfactorily described in 
terms of a uniform turbulence spectrum on the Sun. In order to investigate 
effectively the non-random features which seem to be present in the latest 
photographs (see Section 1) it will be necessary to secure statistically more 
massive data than have been used in the past, and especially to evaluate the 
distributions over the solar disk in two dimensions so far as is possible. The 
brightness data discussed in the literature have been obtained by making 
point-scans across photographs. ‘The evaluation of these in two dimensions 
would require no new observational procedure ; the need is for improved methods 
of measurement and interpretation (Section 2). The method of two-dimensional 
auto-correlation by double transmission suggested by Uberoi (1955a) and by 
myself (1953, p. 278) should reduce the computations to manageable proportions. 

The problem of observing two-dimensional velocity fields, and especially 
simultaneous velocity and brightness fields, with the necessary resolution appears 
formidable. Possibly the best that can be done is to secure slit-scans, the 
advantages of which over point-scans have been discussed in Section 2.6. The 
method would be to scan an image of the Sun across the slit of a spectrometer 
in which the averages over the slit-length of the brightness and of the Doppler 
velocity are simultaneously measured by photoelectric means. It seems the more 
appropriate to call attention to this method of observing because its disadvantages 
are not quite self-evident. 
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ON THE H-FUNCTIONS FOR ISOTROPIC SCATTERING 
D. W. N. Stibbs and R. E. Weir 
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Summary 


It is shown that the H-functions for isotropic scattering in a semi- 
infinite plane-parallel atmosphere may be readily and accurately obtained by 
direct quadrature of an integral form. Methods are described which enable 
the accuracy of the calculations to be maintained in the neighbourhood of 
singularities in the integrand and its first derivative. The final form has been 
programmed for the IBM 704 computer, and a total of 160 H-functions have 
been obtained for values of the albedo 


™ =0 (0°01) 0°90 (0°005) 0°95 (0°001) 0°99 (0°0005) I°0. 
Good agreement has been found in checks involving the zero-order moment 
and the first and second moments. The H-function obtained for the 
conservative case m= 1 agrees to within one unit in the sixth decimal place with 
that computed by Placzek from a modified form of the Wiener-Hopf 
integral. 

Tables of the H-function are presented for 4 =0 (0°05) 1-0 and selected 
values of mw, the tabular values being rounded to the sixth decimal place. To 
enable H-functions to be obtained for any required value of w, the complete 
set of computed functions has been approximated for ~=0 (0°05) 1°o by 
polynomials in w using Chebychev polynomials and the method of least 
squares, several ranges being used to cover the entire range in w from o to 1. 
The H-functions computed from these polynomials are correct to within one 


unit in the fourth decimal place. The first-order moment has also been 
obtained in polynomial form. 





Introduction.—The H-functions were originally introduced by Ambartzumian 
and later developed extensively by Chandrasekhar. They are solutions of the 
non-linear integral equations of the form 


AE (py \ dt 
H(u)= 3+ nH) {EP Hud (1) 
ontp 
where ’(y) is usually a polynomial of even order in yp satisfying the condition 


[(Y@du<s. 
0 


These functions are of considerable importance in the theory of radiative transfer 
in planetary and stellar atmospheres. The detailed numerical application of the 
mathematical theory depends, however, upon the availability of the functions in 
sufficient tabular detail. This requirement has been partly met by the heroic 
computations of Chandrasekhar and Breen (1, 2) for isotropic and anisotropic 
scattering, supplemented later by similar calculations by Harris (3). Valuable 
though these tabulations have been, there is still a need for more detailed numerical 
work on the H-functions, particularly in a form which would enable the required 
function to be rapidly and accurately computed for any value of the particle albedo. 
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Accordingly, we have re-examined the question of the computation of the 
H-functions with the object of providing for current, and we may hope, future 
requirements for the case of isotropic scattering. The work presented in this 


paper forms part of an extensive investigation on the H-functions for atmospheres 
with coherent and non-coherent scattering. 


1. Previous methods of solution 
When the scattering is isotropic with particle albedo w, we have - 
‘Y(«) = constant = 4. 
The integral equation for the H-function is then 


1 , 
H(y)=1+ foun) [AED ay (2) 
oetp 

1.1. The iterative method.—The integral equation for the H-function may be 
written in an alternative form (2, Section 38) 

See scr ty’ H(p'), , 

Ha) =(1-w) + beef EO ae (3) 
Chandrasekhar and Breen solved for H(y) by iteration in the integral equation (3), 
starting with an approximate form for the H-function given by the third or fourth 
approximation in the method of discrete ordinates. The accuracy of the H-function 
at any stage in the sequence of iterations was judged by a comparison between 
the theoretical value of the zero-order moment and the value obtained by 
quadrature of the computed function. 

The same procedure was followed by Harris (3) who used a digital computer. 
Starting with an approximate function H(u#)=1+2 for all w, and using 
Simpson’s rule for the numerical quadratures over » in the range 0(0°05)1°0, 
Harris ceased iterating when the difference between the theoretical and computed 
zero-order moment became less than 0-00005. However, as mentioned by Harris, 
it is well known that the derivative of H(u) at «=o is logarithmically infinite. 
This singularity should be removed before making the quadratures otherwise 
it is open to question whether a satisfactory check on the zero-moment is a 
criterion for the accuracy of the H-function or whether it represents a 
compensation between errors in the »-quadrature and errors in the function 
itself. That there are such errors is indicated by a comparison which Harris 
made between his results for the conservative case w=1 at 4=0°05 and o-10 
which are in error by +3 in the fourth and +7 in the fifth decimal place 
respectively when compared with the very accurate values given for this case 
by Placzek (4). Apart from the uncertainties introduced by the presence of 
the singularity at »=0, the H-function for w=1 is almost a pathological case 
when handled by the iterative method. 

1.2. The Wiener—-Hopf integral.—In the conservative case, that is when w = 1, 
the H-function is given by the well-known Wiener—Hopf integral (§) 


Hu) = (1+) exp {2 (oC ao}. (4) 





cos* @ + u* sin? 6 
Introducing some elegant transformations, Placzek (4) reduced the integral to 
a form which is well suited to accurate numerical quadrature, namely 
\. s) 


4g SED T (8 Oca (uean) 
HW)= Gamer {2 avr 
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Placeck calculated the function #()« }//(m) to 7 decimal places for values of 


pw o(o't)i-o, and then obtained 4() to 5 decimal places for intermediate values 
of » ~ O(0°01)1o by interpolation, 

lhe computational simplicity of obtaining this H-function for each value 
of » from only one quadrature suggests that the use of a suitable integral form 
for M(j) for any @ should have some advantages over the method of repeated 
iteration in the non-linear integral equation. 
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Integral form for the H-function for non-conservative isotropic scattering 
Chandrasekhar (2) has shown that the solution of the integral equation 
(1) tor (ja) can be represented as a complex integral in the form 


‘= 
In Af (ja) =| In 7 (w) No rt (6) 
2nt ie w fe 


where 


0 * pe 


' "V(n)d 
(a)et ast | (+) dp 
bollowing Kourganott and Busbridge (§, Section 29.5), the complex integral in 
equation (6) can be transformed into a real integral by the substitution @ «7 cot, 
With ‘(ja)e der we obtain 


PUcore) «1 ~~ wl cote, 
and 


— 


“re? In {1 wf cot @) 
in M(u e rd ae ™ 
2 (#) aR cos® @ 4 p® ain? é : (7) 


When aw ~ t this is simply the form given by Kourganott and Busbridge (§, 
equation (29.39)). ‘The corresponding form of equation (7) for a problem in 
non-coherent isotropic scattering has been given by Busbridge and Stibbs 
(6, Section 4). Apart from the use of an approximation to the integral form 
in (6) to obtain starting values of the H-functions for iteration in the integral 
equations of the problem, it appears that no detailed application has yet been 
made of the integral form for non-conservative isotropic scattering. 


;. Analysis of the integral form 


We propose to evaluate the function //(,) by means of the integral form 


H(y, ar) @ exp [1 (ja, @)), (8) 


where 


oe ee x In (1 ~ cot) 44 (9) 
mJo cos*é4 p*sin® 6 
for 
Og nel, Oon@msai 

In general, the integral J(u, a) can be computed numerically by means of some 
suitably accurate quadrature formula. However, the accuracy achieved will not 
be satisfactory unless any singularities in the behaviour of the integrand are first 
removed, ‘There are, in fact, two cases which require special consideration, 
Vhese arise from the fact that both 0H/a and 0H/dm are infinite when po 
and w= 1 respectively, 

3.1. The case when p is small and ai. In this case, the integrand in 
equation (9) increases to a sharp maximum near the end of the range of integration 
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and falls to zero at 0 w/a. In the limit as » tends to zero, the maximum goes to 
infinity. ‘This maximum may be eliminated by integration by parts, Since 
dé 1 
« =tan™' (stané@ 
cos®) + pt aint d MM sit saan 


and the integrated part vanishes, we obtain 


Ku.) = = {"" Hepa (10) 


where 


{(@) ine - wb cor0)} tan” (jp tan 6), (11,1) 


6 cosec® 0 — cot é 
one Oe * (tan @). (11.2) 





Although the function /(@) is now a monotonic increasing function of @, the first 
derivative at @ = w/2 is large when » is small, ‘This is due to the factor tan-(p tan @) 
in ((@). ‘The remaining factor is zero at @ 0, and increases monotonically to 
the value (w/2)er at @«/2. Hence we write 


g(0)  (m/2)em tan (tan @), (12) 


and consider the function 


wa 
Har) © {"" (700) ~ (0) d0-+ © {™ g(0) A, 


«1 (uy, ) + Lyn, @). (13) 

lhe function J,(4,@) can now be obtained with satisfactory accuracy by 

numerical quadrature since the function /(@)~g(@) vanishes at Go and w/a, 

and is well behaved throughout the entire range. ‘The second function J,(y,@) 

in equation (13) can be evaluated to any desired accuracy by means of the 
following formulae* for the integral 


© t+ 


! seek + Se 
L(y, @) jer} } In in =e Seal og P< TO, (14.1) 


nial’ .% soa nan 


‘The series in equation (14.1) is obviously most convergent near so and the 
series in (14.2) near = 1-0, the convergence being equal when p= (1 —)/(1 +), 
that is when p= \/2—1, Hence equation (14.1) may be used for p¢4/2—1, 
and equation (14.2) for p> /2—1 if required, 

3.2. The case when w is near unity.—In this case, the integral 1\(4,@) in 
equation (13) requires further manipulation in order to obtain satisfactory 
accuracy in the quadrature, ‘The difficulty arises from the behaviour of f(@) 
for small values of 6. Expanding in powers of @ 


(0) = §yer(O® — 20° — hut + ...)/[1 er + herd...) (15) 
so that as 6-0, ((@)->ap if a & 1 but f(@)-0 if w #1. 





* The formulae may be obtained by differentiation with reapect to » pte oasis 
@ to give 0/,/Q.< — der (in p(t»). Equation (14.1) then follows by —— 
over ». Equation (14.2) is obtained by expanding In » as a power series in (1-~)(t +a) 
performing the p-integration, 
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When w = 1, there is no difficulty in the quadrature for J,(, a) as the functions 
f(@) and g(@) are well-behaved, but the function f(@) is discontinuous for small @ 
when w>1. This behaviour may be removed by introducing a function A(@) 
which, for small 0, is the difference between {(@) when w #1 and w=1, namely 

h(0) = —2p(1 —w)/[1-w +406]. (16) 
We now have 

wi2 a2 
Luar) = = [™ {f(0)—e(0)—W(0)} 0+ = [oye 
TJo 7TJo 
=1;(u, a) +1,(u, a). (17) 

The function J,(,@) may be obtained with satisfactory accuracy by numerical 
quadrature for all » however small, and for all w however close to unity. The 
function J,(.,~) is an elementary integral which gives 


tine) (ES) olga) 


To the first order, this function represents the difference between J(u, a) and 
I(u, 1-0). Hence, from equations (8) and (18) we obtain the following approxima- 
tion for w near unity, 


Heed mon[ = (HSA) Ilqtal}} 


This expression may be reduced to the approximate form quoted by Harris (3) 
by expanding the inverse tangent and exponential functions, and retaining terms 
to the first order in (1—a)'*. We thus obtain 
H(u,a)~ H*(u,a), 
where 
H*(p, w) I 
= ; (20) 
H(p,t-0) 1+pr/[3(1-w)] 
‘To sum up, the required H-function given by equation (8) may be obtained 
with good accuracy over the range 0 <1 < 1, 0<_a@ <1 by means of equations (13), 
(14), (17) and (18), which give the final form 





Hu,)= = ["* {f(0)—g(0)—W(0)} 40 + I(wver) +1), (22) 


the integral being evaluated by numerical quadrature with f(@), g(@) and h(@) 
given by equations (11), (12) and (16), the functions J,(u, a) and J,(y, a) being 
given by equations (14) and (18) respectively. 

3-3. Moments of the H-functions—The zero-order moment, and the first and 
second moments with respect to » of the computed H-functions may be used 
as an overall check on the accuracy of the computation, using Chandrasekhar’s 
theorems (2, Section 38) that if 


x, = [. H(p,a)u" du (n>0), (22) 


%) = 2a [1 — (1 — a)! ?] (23) 


fora,” + 3(1 — ar)? = 1. (24) 
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In applying these checks, it is necessary that the quadrature over » should be 
as accurate as possible. Sufficient accuracy can be achieved by first removing 
the singular behaviour of 0H/ dp at p=o. 
It follows from equations (8), (13) and (14.1) that in the limit as » tends to 
zero, H(pu, a) ~ F(u,a) where 
F(u,@)=1— doing. (25) 
Accordingly, we write equation (24) for the general moment in the form 


1 1 
2,= | (H(ua)— Fae) u dat i) Fluve)a da, (26) 
0 
performing the quadrature over the first integral, and using the fact that 
a S I w 
[Fee du= + (27) 


It should be stressed that although the function F(u,a) has the same behaviour 
as H(p,m) near »=0, and is effective in improving the ~-quadrature, it is only 
a good approximation to the H-function itself for very small values of p. 





4. Numerical results 

The scheme described in Section 3 for the computation of the H-functions 
and their moments has been programmed for the IBM 704 computer. The 
floating point mode has been employed in which the mantissa has 27 binary digits 
or about 8 decimal digits. Simpson’s formula has been used for the quadratures 
in preference to the formulae of Newton—Cotes or Gauss which were found to 
be no better for this computation. The time required for the computation of 
one complete H-function, using 200 divisions in the range of @ and 40 divisions 
in » was about 45 seconds. A total number of 160 H-functions* have been 
computed for 4=0(0°025)1-0, the values of w being 

a = 0(001 )o-09(0°005 )o-95(0°001 )o-9g(0-0005 ) 1°. 

4.1. The H-functions.——A selection of the H-functions is presented in Table I 
for the restricted interval in 4 =0(0°05)1-0, the values of H(u,a) being rounded 
to the sixth decimal place. With 4 =0(0-05)1-0, the greatest error in the computed 
zero-order moment a, is found to be o-o000051 when w = 1-0, the greatest error 
in the first and second moment check being 00000016 also ata =1-0. However, 
with 4 =0(0°025)1-0 the error in the computed zero-order moment is reduced to 
0-0000001, and in the first and second moment check to o-o000008. This suggests 
that most of the inaccuracy revealed by moment checks for 4 =0(0°05)1°0 is due 
to the interval size in the »-quadrature. Since there may be some cancellation 
of errors from the quadratures over @ and p, the very satisfactory checks involving 
the moments of the H-functions for the smaller interval size in » are only taken 
to indicate, perhaps conservatively, that the tabulations are probably correct to 
the rounded sixth decimal place. Support for this view is to be found in a 
comparison with the accurate values computed by Placzek (4) for the case w= 1 
by means of equation (5), the modified Wiener—Hopf integral, which shows a 
maximum discrepancy at »=1-0 of 1 unit in the sixth decimal place. On the 
basis that the functions involved in the quadratures for other values of @w are 
better behaved than for the case w = 1 while the moment checks are equally good, it 
would appear that we have achieved six-figure accuracy throughout the entire 
range in » and w. se 

* A copy of the complete set of H-functions is available in the Library of the Society. 
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Taste I 
H-Functions for isotropic scattering 
0°05 8 fe) O'ls 0°20 o°25 


1000000 1000000 1000000 1000000 1900000 
1003854 1007809 1011872 1°016053 1°020360 
1°006089 1°012378 1018881 1°025618 1'032608 
1007776 1°O1§841 1024221 1°032943 1°042039 
1°009135 1018644 1°028559 1°038917 1°049762 
1°010271 1°020992 1°032204 1°043954 1°056298 
1°O11243 1°023005 1'035337 1°048295 1°061945 
1°012087 1°024759 1038073 1°052094 1066900 
1°012830 1°026306 1°040491 1°055459 1°071297 
1013492 1027684 1042648 1058466 1°075236 
1'014085 1°028922 1044589 1061176 1078791 
1014621 1'030041 1046346 1063634 1082019 
I‘O1§107 1°031060 1°047947 1°065875 1°084968 
1°O15552 1°O31991 1°049412 1°067929 1087674 
1°015§960 1°032846 1°050759 1:069820 1090167 
1°016336 1033634 1°052003 1071567 1°092474 
1016683 1°034364 1°053155 1°073186 1094615 
1'017006 1035042 1'054226 1°074693 1096607 
1°017306 1035673 1055223 1'076098 1°098468 
1017586 1°036263 1056156 1°077412 I*100209 
1°017848 1°036815 1°057029 1078644 1101843 
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o°25 1069298 1°083033 1097592 1°113087 *129653 
0°30 1076365 1091646 *107899 "125258 143889 
0°35 1:082580 1°099243 “117017 -136060 156568 
0°40 1088109 1°106017 "125169 145745 -167972 
0°45 1093072 ‘112111 *132519 154500 *178307 
"50 1°097559 ‘117631 "139192 162466 ‘187735 
0°55 ‘101641 122662 145285 169755 1°196381 
060 "105374 *127272 "150876 176456 1°204347 
0°65 *108805 "131513 *156030 -182644 I*211717 
0-70 111971 "135432 "160799 188379 1'218559 
o"75 114902 *139066 "165227 193712 1°224932 
080 117626 *142446 *169351 1-198686 1°230885 
0°85 120164 "145599 “173203 1203337 1°236460 
0°90 122536 "148549 "176810 1-207698 1°241693 
0°95 124758 "151314 *180196 1°211796 1°246617 
1°00 126844 "1§3913 *183380 1'215654 1°251259 
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First moment 0°553121 0°564221 0°576212 0°589240 0°603484 
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TaBLe 1—continued 
0°55 0-60 0°65 o"70 o"75 


I°000000 1*000000 1°000000 1*000000 1I-oc0000 
1'049655 1°055317 1°061295 1067655 1'074483 
1°081576 1°091348 1°101788 I*113032 1°125274 
‘107143 1°120442 1'134770 =: 1°§0344 «=: "167476 
128729 1*145164 1°162991 1°182516 1°204176 
147462 1°166733 1°187758 1*210934 1°236832 
164003 1°185868 1'209842 1°236419 1°266306 
178790 1'203044 1°229757 1°2§9517 1°293174 
192133 1218600 1°24.7867 1°280619 1°317848 
"204263 1°232789 1°264447 1*300018 1°340641 
"215359 1°245807 1*279709 1°317945 1*361796 
"225560 1°257808 1'293824 1°334581 1°381510 
'234980 1'268920 1*306929 1*350079 1*399943 
°243713 1°279245 1°319139 1°364562 1°417230 
'2§1837 1'288870 1*330550 1*378135 1°433485 
259418 1'297870 1°341244 1°390889 I 

‘266510 1*306307 1*351291 1°402900 1°463278 
‘273163 1°314235 1*360750 1°414236 1°476973 
‘279418 1°321701 1°369675 1424956  — 1°489958 
‘285311 1°328746 1°378112 1°435110 1°502288 
1290874 1*335406 1°386102 1°444745 I*S14015 


0°619179 0°6 36632 0°656264 0°678667 0°704720 
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080 0°85 0"90 0°925 0°950 


1°000000 1°000000 1°000000 I°000000 I-000000 
1'081914 1°OQO171 1°099678 I*lO§17§ I°111508 
1°138808 1°1§412%0 1°172143 1°182774 1°19§232 
1°186640 1°208619 1°234918 1°250664 1°269354 
1°228639 1°257012 1°291434 1°312301 1°337337 
1°266322 1°300863 1°343271 1°369262 1°400740 
1°300588 1°341089 1*391350 1°422458 " 
1°332034 1°378302 1°436280 1472486 

1°361089 1°412941 1°478496 1°519773 

1°388080 1°445338 1°518327 1-564638 

1°413262 1°475756 1°§56033 1607335 

1°436842 1°§04409 1°591826 1°648067 1719155 
1°458989 1°531471 1°625880 1°687006 1°764711 
1°479848 1°557092  1°658344 1°7242906 = 1-BoBs93 
1°499539 1°581400 1°689347 1*760061 1°850917 
1°518169 1604504 1*719000 1°794411 1891784 
1°535827 1626502 1°747401 1°827441 1'931286 
1°§52595 1°647478 1°774637 1°859238 1969501 
1°568542 1°667507 1*800787 1°889877 2'006§03 
1°583731 1686658 1825919 1°919428 2°042357 
1°598219 1°704991 1°850098 1°947955 2077123 
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TaBLe I|—continued 
0975 0980 0985 


1*00000 ’ 1000000 1°000000 1°000000 
1°114406 1°119418 1°121367 1°123542 
I*20101§ "207516 1211146 1°21§132 1°219616 
1'278122 288052 1°293631 1°299786 1°306750 
1°349185 "362688 1°370314 1°378763 1°388366 
1°415753 "432958 1°442720 1°453574 1°465962 
1°478704 "499723 1°§11700 1°§25062 1°§40368 
1°5 38604 "563535 1°577797 1°593757 1612105 
1°595848 624776 1°641386 1°660026 1°681526 
1650732 ‘683729 1°702742 1°724137 1°748892 
1°703485 1°740613 1°762077 1°786293 1°814396 
1°754293 1°795603 1°819560 1846656 1878193 
1°803307 1°848842 1°875329 1°905 360 ° 
1°850657 1°Q00451 1°929501 1962513 

1°896454 1°950536 1:982176 2018211 

1940796 1°999186 2°033439 2°072534 2118516 
1983767 2°046481 2°083366 2°125554 2°175204 
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2°065900 2°137284 2°179480 2°227934 2°285329 
2°105192 2180915 2°225782 2°277402 2°338686 
2°143380 2°223439 2°270984 2°325788 2°390993 


0°923547 0°949349 0°964480 0°981748 1002057 
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*224875 1*227978 
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“000000 1°000000 1*000000 
129278 1°131483 1°136575 
231624 1236318 1°247350 
"325613 1°333076 1°350834 
414626 1°425121 1°450351 
"500126 1°513902 1°547326 
“582909 1600205 1°642522 
634185 1°647509 "663467 1684510 1°736403 
"707503 1°723233 “742128 1°767136 1°829275 
778913 1°797153 1819124 1°848307 1°921349 
‘848601 1°869448 1°894628 1°928189 2012778 
‘916709 1°940257 1°968772 2°006905 2°103677 
‘983355 2°009690 2°041661 2°084554 2°194133 
2048633 2°077838 2°113379 2°161213 2°284214 
2°112624 2°144775 2°183995 2236945 2°373974 
2°175394 2°210565 2°253569 2°311804 2°463459 
2°237003 2°275263 2°322151 2°385832 2°552704 
2°297503 2°338918 2°389784 2°459069 2641738 
2°356941 2°401571 2°456505 2°531546 2°730587 
2°415358 2463260 2°522349 2°603291 2°819271 
2°472792 2°524021 2°587346 2°674331 2°907809 
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4.2. Polynomial approximations.—The H-functions are frequently required 
for values of w other than those tabulated. To facilitate the computation of the 
functions for intermediate values of w, polynomials of the form 

H(p, 7) = Ay + Aj(w —wo) + A,(w —w)* +..., (28) 
have been obtained for each value of » by the method of least squares using 
Chebychev polynomials. For convenience, wy was taken to be the mid-point 
of the range of w for which H(z, a) was being approximated. The degree of the 
polynomial was determined by specifying a tolerance on the maximum permissible 
magnitude of the residual at any tabulated point, the polynomial of lowest degree 
within the specified tolerance being adopted. A tolerance of 0-000025 was 
chosen in order to ensure that the polynomial approximation would give 
H-functions accurately to four decimal places. It is, of course, impossible to 
fit a polynomial over the whole range of w due to the increasing curvature of the 
H-function with increasing w, culminating in the singular behaviour of 0H/dw 
at w=1°0. Accordingly, the range was subdivided into the ranges w =0 to 9-4, 
0-4 to 0°6, 0°6 to o°8, and o°8 to og before fitting the polynomials by means of 
equation (28). With this subdivision, the tolerance was satisfied by quartics 
when » = 1-0 and by polynomials of lower degree for smaller values of pu. 

In order to obtain a sufficiently accurate polynomial representation in the 
range w =0°g0 to I°0, it is necessary to remove the singular behaviour of 6H/dw 
before approximation. The form 

H(y, a7) — H*(u, 2) = By + B(1—w) + B(1—w)* +..., (29) 
where the function H*(, a) is given by equation (20), was found to be satisfactory 
and was applied with the tolerance 0-000025 in the separate ranges w=0-90 
to 0°95, 0°95 to O-gg, and 0-99 to 1-0. 

The coefficients for the polynomial approximations (28) and (29) are given 
in Tables II and III. These polynomials reproduce the tabular values of the 
H-functions with an accuracy better than 0000025. It has been verified by some 
additional computations of the H-functions from the integral form that the errors 
at intermediate points do not exceed the required tolerance. This indicates that 
the polynomials are satisfactory functional forms in the chosen ranges, and that 
the H-functions calculated from them are only in error by round-off in the fourth 
decimal place. 

The first-order moment «,(a) of the H-functions has also been approximated 
by polynomials for the same subdivision of the range in w and the same tolerance, 
the coefficients being given in Tables II and III. The singular behaviour of 
dx,/a at w=1 has been removed by means of a function «,*(e) which is an 
approximation to the first-order moment of the function H*(y,a). Using 


a,°(e)= | (1 —py/[3(1 ew) Au, 1-0) dy 


= 0 (1-0) — a9(1-0)4/[3(1-@)] (30) 
with «,(1°0)=1°154700 and a,(1-0) =0-820352, the first-order moment has been 
approximated in the ranges w =0-go to 0°95, 0°95 to 0-99, and 0-99 to 10 by the 
polynomial 

a, —a,°=B,+B,(1-w)+B(1-—w)+.... (31) 
Only the first-order moment «, has been obtained in this way since the 
second-order moment a, may be obtained from «, by means of equation (24), 
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Taste II 


Polynomial approximations to the H-functions and their first moments 
Hp, w) 


a, (a) 


} =A,+A,(w-—w,)+A(w—a,)?+... 
o<weso4, Wy=0'2 


Ay A, A, Ay; 
0°08484 0°0257 O'O14 
0°13718 O'O514 0°030 
0°17805 0°0760 0°047 
021186 00991 0°065 
0°24071 0°1208 0081 
0°26581 0°1382 07098 
0°28796 0°1566 O'rl4 
0°30771 0°1738 0130 
0°32548 0°1899 O°145 
0°34158 0°2050 O°159 
0°35626 o"2192 0°173 
0°36970 0°2325 0°187 
0°38207 0°2451 0*200 
0°39350 0°2569 o'212 
0°40409 0°2681 0°224 
0°41395 0°2787 0°235 
0°42314 0°2887 0246 
0°43174 0°2682 0°256 
0°43980 0°3073 0:266 
0°44737 0°3158 0°276 
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TaBLe [l—continued 
o6<aw<o8, w,=0'7 


A, A; 
O°1315 0096 
0°2340 0°204 
0°3258 0°319 
0*4100 0°423 
0°4883 0°539 
0°5616 0°655 
0°6306 o"'771 
0°6956 0886 
0°7571 0°999 
0°8154 
0°8707 
0°9233 
0°9735 
1°0213 
1:0669 
11106 
1°1523 
1°1923 
1°2306 
1°2674 





0°4808 


o8<w<o'9, wy=0'85 


Ay A, Ay 
1°OQ017 0°1756 0°250 
I°I54it 0°3290 0°542 
1*20861 0°4757 0°864 
1°25700 06178 1°209 
1°30085 0°7561 1573 
1°34107 0°8908 1°952 
1°37830 1°0220 2°254 
1°41294 1°1498 2°636 
1°44534 1°2742 3°024 
1°47576 1°3954 3°416 
1°50441 1°5133 3°810 
1°53147 1°6281 4°206 
1°55709 1°7398 4602 
1°58140 1°8485 4°997 
1°60451 1°9542 5°390 
1°62650 2°0572 5°781 
164748 2°1573 6°170 
166751 2°2547 6555 
1°68666 2°3496 6-937 
1°70499 2°4419 7314 


0°77438 0°8713 2°369 
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Taste III 


Polynomial approximations when w is near unity 


Hp, w)—H*(p, w) 


a, (ar) — x, *(a) 


\ = Bo+Biu-—w)+BA1 —am)?+ eee 


o'90 [ BW <O'9S 


Bo B, : B, 
— 0°00037 
—0°00026 
— 0°00043 
— 0°00063 
— 0°00085 
—0°00109 
— 0°00134 
— 000160 
—0°00187 
—0°00214 
— 000242 
—0°00270 
— 0°00298 
— 0°00326 
— 0°00354 
— 0°00382 
— 0°00409 
— 0°00437 
—0°00464 
— 000490 


0°00371 


0°95 <M <0'909 


B, 
—0'0681 
— 01209 
—0°1603 
—0°1946 
—0°2254 
— 02532 
—0'2976 
—0°3248 
— 0°3504 
— 0°3745 
—0°3971 
—0°4186 
—0°4388 
—0°4579 
—0'4761 
—0°4931 
— 0°5093 
—0°5246 
— 0°5393 
—0°5528 


1°459 
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Tasie I1]—continued 
o'99 SW < 1°00 


ad Bo B, B, 
0°05 — 0°00000 —0°074 oe 

o"10 — 0°00001 —O°124 

O15 — 000002 —o'165 ee 

0°20 — 0°00000 —o216 14 

0°25 — 000001 —0°255 "9 

0°30 — 0°00001 —0o'291 2°5 

0°35 — 0°00001 —0°326 3°0 

0°40 — 0°00001 —0'360 3°7 

0°45 — 0°00001 —0°392 43 

0°50 — 000002 —0°423 50 

0°55 — 0°00002 —0°453 5°7 

0-60 — 0°00002 —0°483 6°4 i 
0°65 — 0°00001 -- 0536 13°4 —415 
0°70 — 0°00001 — 0566 14°9 — 466 
o'75 — 0°00001 —0°597 16°5 —§21 
0-80 — 0°00001 —0°626 181 —575 
0°85 — 000001 —0°656 19°8 — 633 
090 — 000001 — 0684 21°5 " —693 
0°95 — O°00001 —o°713 23°3 ~—755 
1'0O — 0°00002 —0°740 25°! —819 


a, —a,* 000001 1°564 —152 —444 


and the zero-order moment a, is given explicitly in terms of w by equation (23). 
The value in having «, «, and a, readily calculable for any @ lies in the fact 
that the flux of radiation can often be written down in terms of these moments. 

Conclusion —The scope of this paper has been limited to a discussion of the 
H-functions which arise in the case of the coherent isotropic scattering of radiation 
in a semi-infinite plane-parallel atmosphere, and to the presentation of the results 
of calculations ina form suited to the need for rapid computationof the H-functions 
and their moments for any value of the particle albedo. The methods described 
in the paper may, however, be extended with further elaboration to the scattering 
functions (6) which arise in the theory of the formation of multiplet lines by 
non-coherent isotropic scattering. 
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THE DWARF BINARY HD 16157: AN INTERIM REPORT 


David S. Evans 


(Communicated by H.M. Astronomer at the Cape) 
(Received 1959 April 11) 


Summary 


Spectroscopic observations of this remarkable star have been carried on 
since 1951 using the equipment of the Radcliffe Observatory, Pretoria. It 
now seems tolerably certain that it is a spectroscopic binary, of which the 
primary component is a late K-type dwarf, having an orbital period of 
1°56145 days, and a mass function with the very low value of 0-00871 solar 
masses. On a few plates evidence of the spectrum of the secondary compon- 
ent is thought to have been found. The parallax has been measured as 0-083. 
The spectra resemble those of YY Geminorum in showing Balmer lines 
and Calcium lines in emission. The star also resembles YY Geminorum 
in the fact that it is also a light variable showing changes in the light curve 
from season to season which are even more pronounced than those shown by 
the latter. It seems difficult to account for these in terms of the ordinary 
eclipsing binary model, and the data which have been obtained so far are 
presented with the object of encouraging independent observations 
elsewhere. 





The star HD 16157 (=CZC 1209), R.A. 2" 32™-5, Dec. —44° o1’ (1950), 
was included in the Cape radial velocity programme with the Radcliffe reflector 
because of its high parallax. This amounts to 0’-083, and is the mean of two 
determinations (1), the relative parallax being given as 0”-084. In the Henry 
Draper Catalogue, the spectral type is given as K5, but a re-determination from 
the Cape radial velocity plates gives K7 V, or possibly even a shade later, but 
not as lateas Mo V. In addition to the typical absorption spectrum, the following 
emission lines are found: Balmer series, alpha to zeta; Cal, H and K. Both 
the radial velocity and the magnitude of the star are variable. 

With a few exceptions, the spectra show only single lines, the derived velo- 
cities ranging from +5 to +8okm/s. The velocity determined from the 
absorption lines is always essentially identical with that derived from either group 
of emission lines, and deviations are confined to cases where the absorption 
spectrum is underexposed and the velocity derived from it inaccurate. The 
results of the radial velocity measures are shown in Table I. Here the Julian 
dates have been corrected for light-time to the Sun. The measured velocity 
from the absorption lines is accompanied by its probable error and the number 
of lines measured. The number of Balmer emission lines measured is also 
indicated. ‘The adopted velocity for each spectrum is, usually, the absorption 
line velocity corrected to the system of the Mount Wilson General Catalogue of 
Radial Velocities in the way described in a previous publication (2). Excep- 
tionally, when the absorption spectrum is weak, the emission lines have been 
taken into the adopted velocity to give greater weight. These cases are indicated 
by a colon, and the resulting velocity assigned half the normal weight for a good 
observation. All the spectra were obtained using the c-camera (49 A/mm) with 
exposures ranging up to 3 hours in extreme cases, the normal exposure being 


about 60-90 minutes. 





Nd. 5, 1959 The dwarf binary HD 16157: an interim report 


Tasie I 
Results of radial velocity measures 


Measured velocities Adopted 
ae " ——S ee 
Absorption lines H&K ~e M.W.C. 
System 
4250°596 0°2095 + 3°0+4°8(7) +18 +41(2) +11: 
4261°547 "2224 +37°7 2°8(4) +25 +2913) +31: 
4297°472 "2301 +230 1°3(12) +21 +32(3) +218 
5444°327 “7101 +73°8 1:1(14) +68 +67(4) +72°9 
5704°580 "3842 +50°2 1°6(12) +54 +5214) +483 
5708°573 "9412 +17°7 1°8(13) +19 +21(5) +158 
5711°471 “7976 +47°7 2°5(12) +50 +4914) +458 
5855°301 ‘9105 +259 2°1(13) +29 +32(5) +240 
5858-298 82905 +46°3 3°0(12) +50 +57(3) +46°7 
6084°615 7702 +54 § (13) +58 +56(4) +52°9: 
6087-603 6836 +73°1 2°2(13) +76 +745) +712 
6090°580 *§902 +80°4 4°5(12) +82 +84(4) +79°5 
6158-446 "0533 —14 (3) e 6+ = 2 
6172°299 9254 +196 1°3(13) +26 +2444) +17°7 
6181-266 6682 +83°0 1°5(14) +84 +84(5) +8171 
6202°271 "1205 + 64 1°7(13) ° + 544) + 45 
6204309 "4257 +72 (3) +59 +58(2) +62: 
6212°324 5586 +81°5 1°5(13) +80 +84(4) +796 
6444°639 *3400 +43°3 1°4(13) +41 +46(5) +414 
6446°561 “$711 +82°8 2:1(12) +83 +84(4) +809 
6470°537 9261 +25°0 1°6(14) +21 +18(4) +2371 
6479°486 "6572 +78°7 2°2(14) +77 +81(4) +768 
Secondary —27K -29 S 
6486-437 “1090 + 79 1°8(r1) + 6 + 73) + Go 
6515°412 6651 + 82-1 +1°1(13) +80 +8144) +802 
Secondary —24(1) —25°4S 
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NN 


Check observations : 
6612-291 “7095 +71°8+3°1(5) eae +56(2) +65: 
6614°253 9660 +149 2°6(13) +12 + 443) +13°0 
Secondary ?+95K +93 S? 


Trial and error have yielded a period of 1-56145 days from the spectroscopic 
data, (P-!=0-64043) and phases have been reckoned, using this period, from 
J.D. 2430000°0. The last two observations in the table were not included in the 
solution for the orbital elements, being obtained subsequently when it became 
necessary to make as certain as possible that the correct elements had been 
obtained. For the derivation of the orbital elements the following normal points 
were used, the period being taken as known :— 


Taste II 


Velocity Phase Weight Computed V 
+ 3°80 0°1025 + 5°54 
+21°40 "2230 + 18°69 
+ 48-28 “3748 +49°77 
+ 80-10 5699 +79°79 
+79°37 6635 +7706 
+72°05 6969 +72°84 
+47°58 “8049 +4915 
+20°15 + 20°03 


our QAuut wu 
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The last column represents the residuals of the observed normal points from the 
velocities computed on the basis of the derived elements. These were found in 
the usual way by computation of corrections to preliminary elements, the results 
being :— 
P=1-56145 days, y= + 41°94 km/s, K =37°77 km/s 
e€=0°0461, w= 63°-83, do= 0°7635 


where ¢, is the phase corresponding to passage through periastron. The two 
check observations fit these elements well, the residual, O-C, having the values, 
respectively, —5 km/s and +2km/s. The adequacy with which the data are 
represented can be judged from Fig. 1. 
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Fic. 1.—Radial velocity observations of HD 16157. Filled circles: absorption line measures 
of primary. Vertical crosses: emission lines attributed to secondary. Open circles: halfweight 
observations of primary (see text). Diagonal crosses: check observations. 


Table I contains three entries marked ‘‘Secondary’’. These refer to 
plates on which one or more of the emission lines have faint secondary com- 
ponents. In the best case three such components can be detected. Careful 
examination of the absorption spectra has led to suspicions that, in certain cases, 
duplicity of some lines, notably the neutral calcium line at 4226 A, can be detected, 
but there is nothing which inspires a feeling of certainty. In view of this, it is 
best to adopt an attitude of reserve towards these faint components, and not to 
assume too readily that they belong to the secondary component of the binary. 

We can begin by computing from the derived orbital elements the mass- 
function of the system, which can be expressed in terms of the unknown masses 
and inclination : 

m’® sin? i 
(m+ m’)® 
where m and m’ are the masses of the primary and secondary respectively. 


= 0-00871 solar masses 
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If, for the moment, we assume the secondary to be both considerably fainter 
and less massive than the primary, we can temporarily ignore its existence. 
Then, if we take from the photometric data (dealt with in detail below) an 
apparent visual magnitude of the order of V = 8-76, and a value of B— V = 1-39, we 
find from the calibration by Evans, Menzies and Stoy (2), a type of about Kg V, 
and M,= 8-47, leading to a photometric parallax of 0-087 (compare the trigo- 
nometrical value 0”-083 quoted above and the estimated type of K7 V). This 
system of values is quite self-consistent and will suffer no serious change if we 
assume, as we legitimately can from the experience with the spectra, that the 
secondary is at least two magnitudes fainter than the primary. 

From the discussion of the mass-luminosity relation for late-type dwarfs 
given by Eggen (3), we can infer a mass for the primary of about 0-47 solar units. 
Introducing this value for m into the mass function we find, for different values 
of sini, the following corresponding quantities : 


sin 7 m’ 

I°000 O'I50 
"900 «=*170 
‘B00 = 196 
“700 = 230 
‘600 +285 
"500 “365 


One can hardly assume that the fainter component is the more massive and one 
is, therefore, driven to the conclusion that sini is, probably, greater than about 
o-5 and that the mass ratio is greater than unity and less than about three. This 
is a range into which the data will fall if we decide to take the faint emission com- 
ponents as belonging to the second star of the system. If we do this we easily 
find K’ = 73°3 km/s, giving a mass-ratio of 1-94, and we have 


asin i = 810,000 km a’ sini = 1,571,000 km 
m sin® i = 0-1462 m' sin® i = 0°0753 


Moreover the measured value of the faint emission component on the last check 
plate is + 93 km/s while the computed value is +98 km/s. The evidence is not 
strong but is, at least, not discouraging. 

The assumption of m=0-47 then gives m'=0-24, sini=0-678, i=42°-7, 
a=1 200000 km and a’=2310000km. This set of figures creates difficulties 
for us. We shall, presently, have to consider whether this can be an eclipsing 
binary system, and, if this is to be possible, either the stars must be large enough 
to be almost in contact, or the inclination of the orbit must be near go°. It seems 
impossible to secure these conditions. According to Allen (4) a dwarf with 
M,=8:-5, should have a radius of about 63 per cent of that of the Sun, that is, 
about 440000 km, which, unless the secondary has a radius grossly inflated in 
comparison with what might be expected, is totally inadequate to bridge the gap 
of about 3 million km between the star centres. Even if we took the bull by the 
horns and postulated that i=go° we should still have about the same distance 
between the star centres. This follows because the expression for @sini, 
depending, as it does, only on constants derived from the orbit of the primary, 
remains unchanged and our assumption thus gives a=810000km. We must 
also accept the consequences of the assumption, that is, a mass ratio of 3-1 and 
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hence a value of a’ =2540000km. ‘Thus, if we throw out the evidence from 
the secondary, we can force the system into being an eclipsing binary but, if the 
star radii are normal, any eclipses there can be will occupy only a very small 
fraction of the period. 


The point of the foregoing discussion will become clear as we turn to a dis- 
cussion of the photometric data. These are set out in Table III. 


Taste III 
Photoelectric observations of HD 16157 made with the McClean refractor 


(B— V)=constant = 1°39. 
(U—B)c=refractor UV colour =constant = 2:19. 


Julian date Phase V Julian date eee 
2430000 + 2430000 + 
5759°457 0°5290 8-92 5815°390 0°3502 
5773°395 4554 8-93 5815°433 ‘3778 
5776°404 "3824 8-91° §815°481 “4085 
5777°409 ‘0260 8°77 5818-333 ‘2347 
5779°359 ‘2749 8-86 5819°348 "8846 
5779°476 "3498 8-88 5820°403 "5602 
5780°387 "9332 8-76 §821°319 "1467 
5783°337 "8225 8°75 5821-359 "1724 
5783°407 "8673 8-76 5821°430 2178 
§785°312 0874 8-76" 5828-319 ‘6290 
5785°357 ‘1162 8-79 5828-357 6534 
5786-405 7874 8-78 §829°326 ‘2739 
5788-352 "0343 8°77 5829°383 "3104 
5789°332 6619 8-81" 5837°354 "4144 
5791°342 "9492 8-76 5837°413 "4522 
5792°406 ‘6306 8°84 5841°358 ‘9783 
5793°341 ‘2294 8°83 §841°441 "0315 
5793°433 ‘2883 8-88 5842°344 ‘6097 
5796°294 "1206 8°74 5842°394 "6417 
5797°371 "8103 8-77 5843°347 ‘2519 
5798388 “4616 8°93 5847°343 ‘8107 
5799°348 "0764 8-77 5848°345 "4523 
5800°378 7361 8-76 
5801°290 *3202 8-go0° 
5801°319 *3387 8-92 6105°477 "1016 
§801°345 "3554 8°94 6142°379 ‘7311 
5801-367 *3695 8-95 6161-388 "9031 
5801°388 3829 893° 6162°400 “5511 
5§801°406 ‘3944 8-94 6168-409 "3989 
5804°289 "2408 8-86 6169°426 "O501 
5804°366 ‘2901 8-90 
5807‘291 +1634 8°84 
5807°357 "2056 8-84 6486°541 "1413 
5811°370 7757 8-75 6497°501 "1264 
5812°298 *3700 8-94 6532°470 "5181 
§812°397 "4334 8-95 6546°446 "4673 
5812°459 "4731 8-91 6553°391 "9144 
§813°358 ‘0489 8-76 6570°345 *7706 
5814°299 “6515 8-83 6583°324 “0814 
5814°331 ‘6720 8-81 6584°335 “7288 
5814°372 "6983 8-78" 6589°306 "9119 
5814°447 7463 8-76 6595°305 ‘7532 
5814°493 ‘7758 8-76 6598-307 *7538 
§814°532 ‘8007 8-76 6599°308 "3948 
§815°335 0°3150 8-92 6604°303 0°5938 
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In this Table the Julian dates have been corrected to the Sun. The phases 
have been computed as were the phases of the radial velocity observations. The 
magnitudes were ‘usually determined by comparison with the nearby star CZC1204, 
but the observations have also been reduced independently by comparison with 
standard E-region stars, and it is clear that the comparison star is not itself a 
variable. It has been found convenient to group the observations according to 
the season (1956-7, 1957-8, 1958-9) in which they were made. 





8.7 


~” 
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SEASON 1956-7 























x 6 
Fic. 3.—Visual magnitudes of HD 16157. 


These results are plotted graphically in Figs. 2 and 3. The observations 
for the season 1956-57 are not too intractable; we see a single very broad 
depression in the curve, which might be attributable to an eclipse. The depth of 
the minimum is about 0-18 magnitudes, and is pointed, not flat, at the minimum 
phase. The position of the minimum does not coincide precisely with the phase 
of passage of the velocity curve through its mean value. However, the phase 
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difference, approximately 0-07, is not so large that it might not, conceivably, be 
explained away by some special hypothesis. The real difficulty presented by 
this light curve is the breadth of the depression, which lasts for something like 
0-7 of the period, and the absence of any secondary eclipse. "These difficulties 
may not be insuperable. We can imagine that the radius of the primary star has 
its normal value and that the absorption spectrum is produced in its outer layers. 
We can further imagine that it has, in addition, a very extended tenuous atmo- 
sphere, which is the source of the emission spectrum, and that it is the radius of 
the latter which counts so far as the production of the eclipse phenomena is con- 
cerned. Finally, if there is such an extended atmosphere surrounding a star of 
small mass, the gravitational field controlling its outer parts will be weak so that 
we could imagine that this atmosphere would not be entirely stable. _ If it suffers, 
possibly erratic, changes of size and form the eclipses will sometimes be early, 
sometimes late, leading to a broadening of the light curve averaged over a season 
of observation. All this is speculative but, even at this stage, we are faced with 
phenomena which do not lend themselves to easy explanation in conventional 
terms. It is not worth while to try to make a detailed analysis of this light curve, 
but, if we assume that the eclipse curve has been broadened by some such effect 
as that suggested above, and that the eclipse is only just total, possible para- 
meters are L=8s5 per cent, L’=15 per cent, r/r’=5 and 1= 48°, the orbits being 
assumed circular. This keeps the depth of the secondary eclipse below the 
range where it would have been certainly detected observationally. 

However, although these parameters may well describe the state of the system 
as it was during that season, the photometric observations obtained during the 
next two seasons show that it has certainly not exhausted its capacity to surprise 
us. ‘These photometric observations are plotted in Fig. 3. Those obtained in 
the season 1957-58 are rather few in number, largely because of the mistaken 
idea then entertained that all that would be necessary would be to define the 
curve more precisely. Instead, in that season the indications are that the un- 
eclipsed brightness has dropped from V =8-76 to V =8-80, that the range of 
variation has been reduced, and that there is a phase-shift of the minimum. 
These are relatively minor changes, however, in comparison with those which 
have taken place by the next season (1958-59). ‘The uneclipsed brightness is 
now V =8-71: the star has been observed as faint as V = 9-03, and, if the rather 
speculative light-curve can be trusted, the minimum may go as faint as 9-08. 
There is no constant phase anywhere, and the phase of minimum is now such as 
to correspond to an eclipse of the primary star by the secondary. 

What are we to make of all this? First, can we be sure that the period used 
in discussing the radial velocity observations is the correct one? The check 
observations, which post-date the latest photometric observation, were secured 
in order to verify this and, as far as can be judged, they rule out the possibility of 
such an error. Clearly more observations of all kinds are desirable. So far it 
seems that the light curve changes rather slowly and that from one season to the 
next it is relatively stable, but it can only be coincidence that we have found it so. 
Sooner or later we may hope to catch the system in the act of changing from one 
kind of light curve to another. 

Any conclusions must be of the most tentative kind, but it does seem that the 
most probable explanation must be sought in the hypothesis of an extended un- 
stable atmosphere which partakes to a greater or less extent in the motion of the 
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component stars. ‘The system would then not be a true eclipsing binary system 
at all, for the actual star disks would not eclipse each other. The possibility of 
the existence of such an atmosphere might be of importance in considering the 
evolution of late-type dwarfs, but it will still be necessary to discover whether 
such a phenomenon is a characteristic only of binaries, or whether it can occur in 
the case of single stars. One would naturally feel more certainty if a precedent 
for behaviour of this type could be found. The closest analogy seems to be 
that of YY Geminorum, where Kron (5) has found evidence of erratic changes 
in the light curve, far less pronounced, it is true, than in the present case. His 
explanation was based on the hypothesis of the occurrence of ‘‘starspots’”’ 
analogous to sunspots, but this explanation does not seem a possible one for this 
system, if only because we have to explain the extraordinary duration of the 
eclipse phenomena. 


I am much indebted to colleagues working at the Cape and Radcliffe Obser- 
vatories for valuable assistance. Dr P. A. Wayman has, at my request, obtained 
some of the spectra at crucial phases. Dr A. J. Wesselink has re-examined all 
the data to ensure that another value of the period could not be found. Professor 
R. H. Stoy has critically examined all the photometric data and re-reduced them 
to ensure that some less bizarre account of the system was not possible. 


The Royal Observatory, 
Cape of Good Hope: 
1959 April. 
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VISUAL AND FAR-RED GRADIENTS AND COLOUR 
TEMPERATURES OF y CASSIOPEIAE, II 


D. R. Barber 
(Received 1959 April 29) 


Summary 


A further series of spectrophotometric gradients of y Cass at mean wave- 
lengths 5120 A and 7ogo A, and corresponding colour temperatures, are 
listed for the period 1951 April to 1958 December. Probable errors of ¢ are 
+0°03 (visual), and +011 (far-red) respectively. 

Although during the six-year interval, 1952-57, $5320 fluctuated between 
extreme values of 0-56 and 1°23, the mean level showed no significant change 
(Ad = +0°03) when compared with its average value for the preceding six- 
year period. 

By contrast, the mean value of $7999 for the same period, 1952-57, 
decreased by as much as 0°38; and there commenced towards the end of 1951 
a cyclic fluctuation of the far-red gradient whose amplitude was ~ 0°75. When 
plotted, these gradient values fit a curve yielding periods of ~ 1-0 yr. in 1951- 
52, and ~1°8 yr. in 1956-57. These changes, undoubtedly real, appear to be 
connected with short-period variations of infra-red intensity in the spectrum 
of y Cass associated with a variable continuum to the red of the Paschen series 
limit. 





Since 1938 routine colour temperature observations of y Cass have been made 
at Sidmouth using the 12-inch McClean prismatic camera; and results for the 
years, 1938-41, and 1945-50 have already been published (1, 2). 

The observations reported here cover the period 1951 April to 1958 December, 
and the gradient values are based on 26 spectrograms of y Cass and the Greenwich 
standard, 5 Cass, photographed in blue, and red regions (mean wavelength, 
5120 A) on Kodak I-E plates; and 16 spectrograms of the same two stars in red, 
and infra-red regions (mean wavelength, 7og0 A) on Kodak I.R.E.R. plates. 
The 12-inch McClean prismatic camera was again used for both series. 

Photometric calibration, and reduction of the spectrograms followed the 
same procedure as that described in Paper I(2). The probable error of a single 
determination of gradient remains the same as that found previously, namely, 
+ 0°03 at 5120 A, and +o0°11 at 7ogo A*. 

In Table I the observed visual and far-red gradients, together with equivalent 
colour temperatures, are listed in chronological order (with Julian date). Those 
entries marked with a colon indicate a lower weight due to uncertain correction 
for atmospheric extinction, difficulty in setting the true continuum level on the 
microphotometer tracing, or imperfections of the photographic plate. 


* These values are derived from direct inter-comparison of Greenwich standard stars. 
Twenty-five stars were used in the visual series; and thirteen stars brighter than 3™-1, including 
the adopted primary standard star, § Cass, formed the far-red series. The P.E. of a single visual 
gradient was found from a discussion of the residuals by the method of least squares. Because of 
the smaller number of far-red gradients a least squares solution was not applied. Instead, the P.E. 
was obtained directly, using mean residuals from repeated gradient comparisons of each Greenwich 
standard star with 5 Cass. 
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Tasze I 


Gradients and colour temperatures of y Cassiopeiae 


J.D. ¢ T(deg.) 
5120 7°99 5120 7°90 
2433 
751 1°08 ian 16 000 
942 1-08 5 16 000 
2434 
O15 1°08 si 16 000 
030 19 100 
043 see 
094 > 50 000 
150 " 16 400 
250 ‘ 27 200 
296 ea 
301 ae. 24 400 
360 s 29 000 


oe 


396 . ous > 50 000 
440 nu 

475 see 13 300 
627 P 14 000 
629 : eS 12 200 
658 . 17 900 
661 . aa 24 400 
686 _ 33 000 


me 


8 
3 
7 
6 
8 
7 
° 
4 


75° 
837 
2435 
007 
046 
o61 


387 
391 
428 


485 
778 


072 


165 
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‘The tabulated results show clearly that for both visual and infra-red regions 
the dispersion in the gradient is always well above the experimental error. 
Comparison of the gradient data for the two wavelengths reveals also important 
differences in the nature of the variation in each spectral region; for, whereas the 
visual gradient showed quite irregular changes of amplitude +0-45 with no 
significant alteration in the mean gradient over the whole period of eight years, 
the far-red gradient exhibited a remarkably persistent sequence of quasi-periodic 
fluctuations of amplitude about seven times the probable error, and average 
period equal to 1-4yr. Maximum and minimum values of $7999 were 1°34, and 
0°37 respectively. The latter corresponds to an equivalent colour temperature 
approaching infinity that is obviously not a measure of the true stellar temperature. 
This is confirmed by a direct comparison of the monochromatic intensities in 
visual and far-red spectral ranges which suggests that the very low gradient 
values measured at 7090 A are due almost entirely to a decreased infra-red flux 
at 7846 A, thus indicating a strong continuous absorption below the Paschen 
series limit. 

If the above explanation is correct, it may be inferred also that the short- 
period variation in ¢y999 is due to transitions from strong absorption to intense 
emission, and vice-versa, in the region of the Paschen continuum. Some con- 
firmation of this is given by the independent observations of Fellgett (3) and Mlle. 
Lunel (4) using Pb-S cells. They report abnormally strong infra-red radiation 
between 0-7 and 2-7» present in the spectrum of y Cass in 1950 March, and 1953 
September respectively. Sidmouth spectrograms taken on close dates also show 
some continuous emission, superimposed on the water-vapour bands beyond 
7846 A, and associated with a rising far-red gradient. 

In addition to the changes in the infra-red region just described, there 
appeared also towards the end of 1951 a very marked discontinuity in the average 
value of 2999; and this is shown most clearly when the mean values of far-red 
gradient for the consecutive six-year intervals, 1946~51, and 1952-57, are directly 
compared. ‘Table II facilitates this comparison by listing gradients and colour 
temperatures for the two periods. In the last column of the Table is entered 
the gradient difference, Ad o99, denoting the departure of the observed gradient 
value from that computed by extrapolation of the visual gradient, $51.9, for the 
same epoch assuming a grey-body energy distribution in the spectrum of y Cass. 


TABLE II 
Comparison of long-term average values of gradient, and colour temperature for the years 
1946-51 and 1952-57 
6-year $5120 Prov0 Ts120 Tr090 Adzo90 
Period (in 10° deg.) (O—C) 
1946-51 0°89 (33) 1°10 (18) 23°0 21°4 +003 
1952-57 0°92 (26) 0°72 (15) 21°4 > 50°0 —0°38 


(The numbers in parentheses indicate the total observations in each group) 


The observed decrease of ~0-40 in the far-red gradient for the second six- 
year period appears in striking contrast to the constancy (within the error of 
measurement) shown by the visual gradient over the same period of time. 

As already stated, this feature can be accounted for by a large, though 
probably temporary, decrease of spectral intensity at 7846 A associated with a 
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varying Paschen continuum. It does not appear likely that Fellgett’s suggested 
modification of the infra-red intensity in the spectrum of y Cass by surrounding 
nebulosity (3) can contribute significantly in this instance. 


Acknowledgment.—I am greatly indebted to Professor H. H. Plaskett for his 
most helpful comments on the MS. of this paper. My thanks are due also to 
Mr S. Archer for his valued assistance in the preparation and reduction of the 
microphotometer tracings from which the results presented here were obtained. 


Norman Lockyer Observatory, 
University of Exeter, 
Sidmouth, Devon: 
1959 April 28. 
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GLOBULAR CLUSTERS, II]. THE SPECTRAL TYPES OF 
INDIVIDUAL STARS AND OF THE INTEGRATED LIGHT 


T. D. Kinman 
(Communicated by the Radcliffe Observer) 


(Received 1959 March 4) 


Summary 


The problem of classifying the spectra of individual stars and spectra 
taken in integrated cluster light is discussed. ‘Types are given to the spectra 
of both kinds previously discussed in Paper I and a comparison is made 
with Mayall’s and also Morgan’s classification for spectra of lower dispersion. 
Mayall’s types show small systematic differences from both the Radcliffe 
and Morgan types: a correction is therefore applied to Mayall’s types to give 
integrated types for 63 clusters on a uniform system. 

The effect of the metal/hydrogen ratio on the metal line strengths in the 
spectra of the cluster giants is discussed and an attempt is made to calculate 
this ratio for different clusters from the present spectra. It is concluded 
that the metal/hydrogen ratio ranges from greater than 1/10 to less than 
1/100 the solar value for different clusters. For a cluster such as NGC 5272, 
it is shown that if this ratio has ~1/10 solar value then the cluster giants 
should have nearly normal spectra for their colour and luminosity as is 
observed but that the metal lines in the spectra of the subgiants should be 
weakened. An attempt is made to predict integrated spectral types for 
NGC 5272 and (with less accuracy) for five other clusters. Differences in the 
integrated spectra are largely due to the differences in the metal/hydrogen 
ratio and to the different relative populations of the horizontal branch. The 
predicted types agree with the observed types within the expected accuracy in 
support of the present analysis. 

A correlation is found between the integrated spectral type of a cluster 
and its concentration to the galactic plane. Clusters of the earliest spectral 
type whose stars have a very low metal abundance form an extended spherical 
distribution with little concentration to the plane, while those of latest type 
whose stars have the highest metal abundance show a strong concentration 
to the plane: clusters of intermediate type, which are the most numerous, 
show an intermediate distribution. The less definite correlation between 
spectral type and the properties of the RR Lyrae stars is also discussed. 





Introduction.—In the first paper of this series (16) details were given of the 
spectra of thirty globular clusters and of the radial velocities derived from them. 
These spectra are classified and compared with previous observations in the 
first section of this paper. The second section deals with the derivation of 
metal abundances from the spectra of the cluster giants and the third section is 
an analysis of the integrated spectra. ‘The correlation of the cluster spectra with 
galactic distribution and with the properties of the RR Lyrae stars in the clusters 
is discussed in a final section. 

1.1. Spectral types for individual stars.—Spectral types taken from the present 
spectra will be uncertain for the following reasons: 
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(1) The stars are faint so that with a dispersion of 86 A/mm at Hy the spectra 
are narrow: for about one half of the stars the spectrum width is 0-25 mm while for 
the faintest stars the width is o-1omm. in some cases the spectrum does not 
extend to the blue beyond A4226. 

(2) The stars are intrinsically red and often further reddened by absorption. 
Prismatic spectra of such objects have a steep density gradient with wave-length ; 
in most cases the range suitable for obtaining types does not exceed 250 A (usually 
MA 4400-4150) and is sometimes less. Grating spectra would be preferable for 
this purpose. 

(3) The spectra are peculiar and do not in general match the spectra of MK 
stars (taken with the same slit width, dispersion and emulsion) which were used 
as standards. Comparison with spectra of other high velocity galactic stars 
proved unsatisfactory, partly because too few of such spectra were available and 
partly because the luminosity classes were not in general comparable. 

The luminosity class may be found for a few spectra by observing the ratio 
4077/A4045. In the majority of cases however, the ratio A4290/A 4325 (cf. 
Popper (23)) was used and also the appearance of the red side of the G band. 
In all cases stars were found to be luminosity class II-Ib in agreement with the 
absolute magnitudes given by colour magnitude diagrams (e.g. (2)). For 
this luminosity class the CN band at A 4215 was in general found to be 
weak, Z 

In these spectra the usual criteria involving hydrogen to metal line ratios (e.g. 
Hy/A 4325) give an approximate type denoted by (H/m) which is earlier than the 
type (m,/m,) derived from the ratios of different metal lines. This difference 
has been interpreted as being due to a metal deficiency in the cluster stars com- 
pared to nearby stars. Attempts have been made to find the type (m,/m,) from 
the ratio Cr1 4 4254/Fe1 A 4250 (cf. Deutsch (6)) but the lines are poorly resolved 
in these spectra; the ratio (A 4226/A 4325) has proved helpful if the luminosity 
class is taken into account. A Hartmann spectrocomparator, kindly loaned by 
the Royal Observatory, Cape, was used for comparing the spectra. 

Following Deutsch (§), it is found that the clusters may be assigned into 
three groups: A, Band C. This classification is one of convenience, however, 
and does not necessarily signify that a continuous range of properties does not 
occur in cluster stars. One cluster in each group is common to both Deutsch 
(38 A/mm) and the present observations and these clusters have served to some 
extent as prototypes. The following properties have been found for each group 
from the Radcliffe spectra. 

Group A: Type (m,/m,)~K2—K4~ Type (H/m); Spectra roughly match MK 
standards although metal lines may be slightly weaker. (Prototype 
NGC 6121.) 

Group B: Type (m,/m,)~G8; Type (H/m) definitely earlier (~Gs5), 
(Prototype NGC 6656.) 

Group C: Type (m,/m,)~G5; Type (H/m) much earlier (~FS8). 
(Prototype NGC 7078.) 

Different stars in the same cluster are all found to belong to the same group, and 

the classification into these groups is given in Table Ia. 

1.2. Spectral types for integrated spectra.—Although the width (o-50-0°25 mm) 
and the available wave-length range (AA 4400-4100) of these spectra make them 


more suitable than the stellar spectra for estimating types, their composite nature 
38* 
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TaBLe la 
Classification of spectra of globular cluster giants 
Group A Group B Group C 
NGC NGC 


c 


104 3201 
288 4372 

Radcliffe 6362 5139 
Radcliffe 6752 

¢ ) 6809 


Radcliffe 


6121 





Deutsch 


7 


5272 6656 
Deutsch 5904 

: 6093 4147 
6205 Deutsch 4 5024 
6218 6254 
6779 
L 7089 








makes this estimation difficult. For the sake of uniformity, Morgan’s (21) 
method of classifying McDonald globular cluster spectra (150 A/mm Hf-He) 
has been followed and three types have been given to each cluster by comparison 
with MK standards. (CH/Hy) is a fundamental type obtained by comparing 
the G band with Hy; (H) is a hydrogen type from the absolute strength of Hy 
while (Fe) is a metallic line type from the absolute strength of the Fe lines in 
the region AA 4250-4400. ‘These types are given for sixteen clusters in Table Id 
and types corresponding to (CH/Hy) types are given for five clusters from 
Newtonian spectra in Table Ic. Different Radcliffe spectra of the same cluster 


match well and suggest that the cluster light was adequately integrated on all 
occasions. Mayall’s (19) types for Crossley spectra agree best with the funda- 
mental type (CH/Hy) for both Radcliffe and McDonald spectra with an r.m.s. 
_ difference of +0°2 to 0-3 class for a cluster (Fig. 1). There is some evidence 
for systematic differences between the low dispersion (430 A/mm) Crossley 


TaBLe Ib 
Spectral types for integrated Casscgrain spectra of 16 globular clusters 
(CH/Hy) NGC (CH/Hy) (H) 
; 6397 Fs Fs5 
6541 F6 Go 
6715 F7 F7 
6752 F6 F7 
6864 Go 
7078 F6 
7089 F8 
7°99 F7 


TABLE Ic 


Spectral type for the integrated Newtonian spectra 

of 5 globular clusters 

NGC (CH/Hy) 

1261 F8 

5139 F7 

5927 G2 

6139 F8 

6584 F7 
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types and the more recent types on higher dispersion. In particular the latter 
show a smaller range largely because the Crossley types are systematically too 
early with respect to the latter for the earliest types. Morgan and Mayall’s esti- 
mates from the same Crossley spectra also show differences of the order of 0-2 to 
o°3 of a class for a cluster. Mayall’s 36” refractor spectra (130 A/mm at Hy) 
agree best with the mean of the (CH/Hy) and (Fe) types with an r.m.s. difference 
of +0-3 class for a cluster. For NGC 6715 however the Lick type differs from 
the Radcliffe type (3 spectra) by og class. Only one cluster (NGC 7078) was 
observed in common with Morgan who obtained F3: F6: [0], compared with the 
Radcliffe F3:F6: <Fo for the (CH/Hy), (H) and (Fe) types respectively. 
From these small samples it is judged that the Radcliffe and McDonald (CH/Hy) 
types form a consistent system to within 0-2-0-3 class and that the Crossley 
spectra may be converted to this system using the curve in Fig. 1. Table V 
gives a list of the (CH/Hy) type for 63 clusters from either the Radcliffe, 
McDonald or corrected Crossley type. 











1 
G5 





Fic. 1.—A comparison of spectral types. The ordinate is the Radcliffe (filled circle) or 
McDonald (open circle) (CH|Hy) type or Morgan’s type for the Crossley spectrum (cross); the 
abscissa is Mayall’s type for the Crossley spectrum of the same cluster. 


Figures 2a and 2b show the variation of the (H) type (filled circles) and the 
(Fe) type (open circles) with the (CH/Hy) type (abscissa) for the Radcliffe and 
McDonald spectra respectively. The general trend of both sets of observations 
agree although they have only one cluster in common. The (H) type is some- 
what later, and the (Fe) type earlier, than the (CH/Hy) type. These differences 
are greatest for the earliest (CH/Hy) and progressively decrease as we move to 
later (CH/Hy) types until for (CH/Hy)~ G3-5 all three types are roughly equal. 
The only difference between the Radcliffe and McDonald types is for the (Fe) 
type for clusters with a (CH/Hy) type F5-8, where the McDonald types are 
systematically 0-4 class earlier than the Radcliffe type: this is somewhat larger 
than the likely error noted above. 
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To summarize: the spectra of the cluster giants may be divided into three 
groups: Group A, for which the spectra are practically normal when compared 
with normal K type giants (9 clusters); Group B, for which the lines are definitely 
weakened (11 clusters); and Group C, in which the lines are very much weakened 
(7 clusters). It must be emphasized that this classification is only one of con- 
venience and does not necessarily signify grouping as distinct from a continuous 
range of parameters. 





T ce 














Fic. 2.—The ordinate is the (H) type (filled circle) or the (Fe) type (open circle). The 
abscissa is the fundamental type (CH|/Hy): (a) Radcliffe and (b) McDonald integrated spectra. 
The integrated spectra show a range in type Gs5—F2 for the (CH/Hy) type, 
- Gs5—Fs5 for the (H) type and Gs5~< Fo for the (Fe) type. A rough relation holds 
between the (CH/Hy) types and the (H) and (Fe) types although Radcliffe and 
McDonald spectra of different clusters differ systematically for the (Fe) type when 
the (CH/Hy) type ~Fs5—F8. Plates 9 and 10 show representative spectra both 
of individual cluster stars and in integrated light; spectra of standard stars are 
also included for comparison. In the next section we consider the interpretation 
of the cluster giant spectra and in the third section the interpretation of the 
integrated spectra. 


2. The interpretation of the spectra of cluster giants 

2.1. The colours and magnitudes of cluster giants.—No colours are at present 
available for the stars from which the present spectra of cluster giants were taken. 
Since however these stars were amongst the brightest visually in each cluster, we 
may assume that they belong to the brightest stars of the red giant sequence. 
Fig. 3 is a composite colour magnitude diagram derived from Arp’s data (2) of 
the bright end of this giant sequence for six clusters. These clusters have been 
combined into three groups of two clusters according to the classification of their 
giant star spectra. ‘The clusters are NGC 5904 and 6205 of Group A, NGC 6524 
and 7089 of Group B and NGC 6341 and 7078 of Group C. These mean curves 
give absolute visual magnitude (M,) assuming (with Arp) that for all clusters 
the RR Lyrae stars have My = —o-10; the vertical lines give the r.m.s. scatter 
of individual stars about the mean curve. 
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Representative spectra of individual cluster giants (a) and MK standards (6). 


T.D. Kinman, Globular clusters, Il. Types of individual stars and integrated light. 
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Representative integrated cluster spectra (a) and MK standards (6). 


T. D. Kinman, Globular clusters, I. Types of individual stars and integrated light. 
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Arp (loc. cit.) suggested that the position of the giant sequence in a colour- 
magnitude diagram depended upon whether the cluster stars had nearly ‘‘ normal ’’ 
spectra or spectra in which the metal lines were appreciably weakened; this 
effect is clearly shown in Fig. 3. The different extensions of the mean curves of 
the three groups may have significance but could be due in part to the smaller 
number of stars observed in clusters of Groups B and C compared with Group A. 
Considering all stars with C.].> +1-20, Group C is o™21+0™-25 visually 
brighter than Group B which is not significant. For the same range of colour, 
however, Groups B and C together are o™-39 + o™-06 visually brighter than Group 
A which on a Student’s ¢-test is highly significant. At C.I.= + 1-3, the absolute 
visual magnitudes of the stars of Groups B and C are ~0™-5 brighter than those of 
Group A, if the RR Lyrae stars in all groups have the same absolute magnitudes. 





-3.-5 T T T ’ T 
Cc 
8 
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Fic. 3.—Composite colour magnitude diagrams from Arp’s data. 
Group A (NGC 5904 and 6205); 
Group B (NGC 6524 and 7089); 
Group C (NGC 6341 and 7078). 


Arp (loc. cit.) also suggested that the mean period of the RR Lyrae stars of 
type ab in a cluster was correlated with the spectra of the giant stars; those 
clusters whose giants have nearly ‘‘normal’’ spectra (Group A) having the 
shortest period RR Lyrae stars. Some support for this suggestion is given in 
Section 4.2. On the pulsation hypothesis (15) (25), this has been interpreted as 
being due to the RR Lyrae stars of Groups B and C being ~ o™-3 visually brighter 
than those of Group A. If this be so, then the My for Groups B and C giants is 
~ o™-§ brighter than those of Group A at C.I.= +1-30. The absolute magnitude 
of field RR Lyrae stars determined by Alania (1) from Pavlovskaya’s proper 
motions is M,y= +0-2+0-2. Eggen,* however, from a study of the common 
motions of groups containing these variables, gives the provisional value My = +0°5 
for RR Lyrae stars of the period associated with clusters of Group A. 

The calculation of the luminosities and effective temperatures of cluster 
giants from their magnitudes and colours is set out in Tabl¢é II. Col. 1 gives the 
colour index (C.I. scale) and col. 2 gives log Terr interpolated for these colours 
from the data given by Morgan and Keenan (20) and Johnson and Morgan (23) 
for luminosity class II stars. These temperatures can only be regarded as 
approximate because only meagre data are available for such stars. Cols. 3 and 4 


*I am grateful to Dr O. J. Eggen for giving me this information in advance of publication. 
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give the corresponding M, from Arp’s data for clusters representing Group A 
and Groups B and C respectively. Cols. 5 and 6 give the same My corrected for 
the assumption that the RR Lyrae stars of clusters of Groups B and C have 
My, = +0°2 and those of Group A have My= +0°5. Kuiper’s bolometric cor- 
rections (18) are applied to these magnitudes to obtain the bolometric magnitudes 


Taste II 
Assumed magnitudes, colours, temperatures and luminosities for cluster giants 





My Mpboi log L/L ro) 





Assumption I 2 2 2 





Cl log T Group Group | Group Group | Group Group | Group Group 
ae a BC A BC A BC A BC 





+o-7 3°694 | —o'§5 —1'2 |} +o'r -—oO'9 |] —O'2 —I1°2 1‘97 2°37 
3°674 | —1°0 —1'5 | —0o4 —1'2 | —O8 —16] 2°21 2°53 
+o9 3°655 | —1'5 —1°8 | —o-Qg —1'5 | —1°4 —2°0 | 2°45 2°69 
3°638 | —2-0 —2:3 | —1'4 —2°0| —2°0 —26] 2°69 2°93 
+11 3622 | —2°93 —26] —1°7 —2°3 | —24 —3°0 | 2°85 3°09 
3600 | —2°5 -—29 | —1'9 —26]| —2a'7 —3°4 | 2°97 3°35 
3°580 | —26 —31 | —2°0 —28 | —2-9 —3°8 | 3°05 3°41 


























Assumption 1: All RR Lyrae stars have My=—o-10 
Pi 2: RR Lyrae stars of Group A have My= +0°5 
RR Lyrae stars of Group BC have My= +0°2 


+07 +08 +09 410 +11 +412 41-3 CLL. 
1 i i qT T 
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3.70 3.65 3.60 
log Te 

Fic. 4.—Comparison of models of globular cluster giants. Upper line of hatched area: 
cluster giants of spectral Groups B and C; lower line of hatched area: cluster giants of spectral 
Group A. 

Curve 1: Hoyle and Schwarzschild (A~!=o); Point 1a: last model of 1 with A~!=7°5 x 10~*. 

Curve 2: Haselgrove and Hoyle (A~-'=o0). 

Curve 3: Kippenhahn, Temesvadry and Biermann (A =7-5 x 10~*, a=1). 

Curve 4: Kippenhahn, Temesvdary and Biermann (A-!=7°5 x 10-*, a=2). 

Curve 5: Kippenhahn, Temesvdry and Biermann (A-'=5 x 107", a=1). 
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in cols. 7 and 8, and cols. g and 10 give the corresponding relative luminosities 
expressed as log (L/L) (assuming Mo; = 4°73). The final column gives the 
normal spectral type for the colour for a star of luminosity class II. The main 
errors present are due to uncertainties in the RR Lyrae absolute magnitudes and 
also in the bolometric corrections. Log (L/L¢) may be uncertain by perhaps 
+ 0-2 and in addition errors are undoubtedly present in the effective temperatures. 
Fig. 4 is a plot of log (L/Lq) against log Terr or C.1.; the upper limit of the 
hatched area corresponds to the clusters of Groups B and C and the lower limit to 
those of Group A. 

2.2. Stellar models of cluster giants.—Stellar models have been calculated for 
cluster giants on the assumption that they have about 1-2 solar masses and an iso- 
thermal helium core which is unmixed with the outer hydrogen envelope. Hoyle 
and Schwarzschild (11) give a sequence of models (curve 1 in Fig. 4) which for 
zero metal abundance (A~!=o) are more luminous than cluster giants but which 
for normal solar metal abundance (A~!~ 10~*) correspond to normal late 
giants. ‘They concluded that cluster giants must have a metal/hydrogen ratio 
of A-!=7-5 x 10-* (point 1a). Hazelgrove and Hoyle (8) repeated these cal- 
culations assuming A~!=o and found a similar sequence (curve 2) but termi- 
nating at a lower luminosity: their identification of this limit with the observed 
end of the giant sequence seems dubious because of the difference in effective 
temperature at this point between the model and the cluster stars. Kippenhahn, 
St. Temesvary and Biermann (17) have repeated Hoyle and Schwarzschild’s cal- 
culations and find that their models depend critically on assumptions made about 
a which is a function of the mixing length in the convective envelope of the star: 
Curves 3 and 4 are their sequences for A~'=7-5 x 10~* assuming a=1 and a=2 
respectively while curve 5 is for A~t=5x 10-7 anda=1. All models agree that 
a lower metal abundance produces a higher luminosity at the same effective tem- 
perature, which is in qualitative agreement with the observation that Group B and 
C cluster giants are brighter than those of Group A of the same colour. It does 
not seem possible to draw further conclusions at present, however, because of the 
uncertainties in both models and the observed colours and luminosities. 

2.3. The metal line strengths in globular cluster giants.—The present spectra of 
cluster giants are not suitable for spectrophotometric analysis. We note, however, 
that the normal spectral type for a luminosity class II star of C.1. ~ +1°3 is about 
K1~—K2 and this is about the same type as we found for the giants of Group A. 
The giants of Group C, on the other hand, have metal lines whose absolute 
strengths correspond to about F8. Using the equivalent widths of strong Fe1 
lines in G and K giants found by Pannekoek (22) and Wright (31) and also 
the photometric atlas of Hiltner and Williams (1x0), it is estimated that the 
weakening of the observed metal lines in Group C giants corresponds to a change 
in equivalent width W given by AlogW=-—o-5. For the giants of Group A, 
we put Alog W< —o-15 which corresponds to the error of our estimation. 

The metal lines observed in low dispersion spectra of late type stars are strong 
lines for which log (W/A)~ —3 and may thus be expected to be formed on the 
damping part of the curve of growth. For such lines 


Ky 


We /Nf (1) 


where N is the number of absorbing atoms per gram of stellar atmosphere, f is 
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the oscillator strength which is constant for a given line and «, is the continuous 
mass absorption coefficient at the wave-length of the line. The damping in a 
strong line will be caused mainly by collisions with neutral hydrogen atoms and 
so the damping constant 6 will be proportional to the gas density p (cf. Unséld’s 
(2g) discussion on the metal line strengths in subdwarfs). At constant tem- 
perature, the equivalent width of the line is thus given by: 


log W =43(log p + log N — log «, + constant). (2) 


For stars of colour index C.I.~ + 1-3, the effective temperature will be taken as 
3800°. Wright’s investigation of late type giants suggests that the excitation 
temperature (7;) for neutral iron lines of low excitation potential should be at 
least 350° cooler than this giving 7;<3450°. In order to find the degree of 
ionization of the metals from the Saha equation, let us consider magnesium as 
representative of an abundant metal with an ionization potential typical of the 


iron group. For 7;<3450°, the fraction ¢ of Mg atoms which are ionized is 
given by: 


log (-*;) < —2:2—logp, (3) 


where p, is the electron pressure. Wright (/oc. cit.) gives p,=0-01 dynes/sq. cm 
for early K giants whilst Miss van Dijke (30) gives 0-003 dynes/sq.cm. It may 
be noted that Ueno (28) gives a much higher value (0-1 dynes/sq.cm) in a 
theoretical model. The metals must contribute most of the electrons and if (as 
we suspect) there is a metal deficiency in the globular cluster giants, this and the 
low surface gravity should mean that the electron pressure is lower in these stars 
than normal giants so that the metals should be largely ionized. As there is 
some doubt in this matter we shall consider two models: one with the metals 
largely neutral and one in which they are largely ionized. 

Let it be assumed that the continuous opacity is caused by negative hydrogen 
absorption. For @,= 5040/7, > 0-9, the mean mass absorption coefficient given 
by Chandrasekhar and Miinch (4) may be represented by 


log k = — 2°83 + 1-580, +logp,. (4) 

In an isothermal atmosphere (cf. (27, equation 14)) the density p is given by: 
log p = log g — log T,, — log & + constant. (5) 
Let ¢ and y be the ionized fractions of the metal and hydrogen respectively, then: 


log p, =log (y+ ¢/A) + logp + log T, + constant (6) 


and for small y: 
log y = — 13-546, + $ log T;—logp,—0-48. (7) 


At constant temperature, equations (4), (5), (6) and (7) yield the following 
relations for the density, electron pressure and continuous opacity :— 


p x git(y+ 4) ale cCKOp a ( +4)" wie (8) 
a A a e ¥ A g& . 


Case I: Metals largely neutral.—In this case (1—$)~1 and dxp,—. Sinceg 
is small we may put ¢=Cy where logl = (xy — Xm)6;- (9) 
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Here x, and x,, are the respective ionization potentials of hydrogen and a metal 
such as magnesium representative of the iron group. We then have: 


£\13 ; ¢\-18 
K, CK OP, K p(s + 4) and p cc g(x + 5) : (10) 


It may be noted that equation (10) and also (16) reduce to those obtained in a 
similar analysis by Schwarzschild, Spitzer and Wildt (27, eqn. 16) if the hydrogen 
ionization is ignored. Putting the values from equation (10) into equation (2) 
we have: 

log W= —4log (A(A + £)*) + constant (11) 
for the equivalent width of a line in the atmosphere of a star at constant g and 7. 
If the hydrogen/metal ratio changes from the solar value Ao to that for a globular 


cluster giant Ag, then the decrease in equivalent width of a strong metal line will 
be: 
A9g(Agt at 
Alog W=} log {SY eh, 12 
eW~t08 eae TOF ~ 
The theoretical models discussed in Section 2.2 suggest that an increase in A 
produces an increase in L/Lo and since 


mf) -m( eee ema) 


it will lead to a decrease in g. Our comparison is therefore between two stars 
having the same temperature and surface gravity; one is the globular cluster giant 
having mass 1:2Mo and a hydrogen/metal ratio Ag and the other is a more 
massive giant of luminosity class II with solar metal abundance Ag. The 
observed differences between the equivalent widths of cluster stars and standard 
stars of the same effective temperature and luminosity are given at the beginning 
of this section ; from this data and equations (g) and (12) we find that the hydrogen 
ionization is negligible and that for 
metals _ 1 


clusters of Group A, a > 3 solar abundance 


metals Par 
H 


solar abundance. 





clusters of Group C, 


‘These low metal abundances indicate that the electron pressures must be lower 
than in normal giants; consequently by equation (3), the metals are likely to be 
ionized and the initial assumption of case I is wrong. Furthermore it is likely 
that weak atomic metal lines contribute to the continuous opacity in these stars. 
We may determine the consequences of this by considering the continuous 
opacity to be wholly due to this cause, in which case: 
log k = log N + constant (14) 

and then from (1), (5) and (14) 

log W= —4 log N+ constant =} log A +constant. (15) 
In this case the strong metal lines actually increase in strength as the metal abun- 
dance falls—an effect which has been considered by Deutsch and others (6). 
The effect will be greatest for the clusters of Group A in which the metal content 
is relatively high, so it is likely that the metal abundances found for this group 
are overestimated. We shall therefore consider the metals to be largely ionized 
in both groups. 





Vol. 119 
Case II: Metals largely ionized.—In this case ¢ ~ 1 and (1 —¢) ocp, and hence: 
K, OK Op, © (y+A-)2gl2 and p x gl? (y+ A-1)-¥2, (16) 
We also have in this case: 
log N= log (1-4) —log A (17) 
log W= —} log (A*(y + A~*)) + constant (18) 


for a star of given g and 7. ‘The decrease in equivalent width in changing the 
abundance from Ag to Ag will be: 


; Ag? +A-1 
Alog W=}log {Sop o+ 3} (19) 


The hydrogen ionization y may be found from equations (4), (5), (6) and (7) to be: 


and hence 





log [y?(y + A-)] = — 27-086; + 5 log 7; + 1-586, — 3-79 —logg. (20) 


Since logg is known from (13), (19) and (20) may be solved for Ag using the 
observed values of AlogW. Assuming Ag =10', the hydrogen ionization is 
found to be negligible and we have for 


metals 


clusters of Group A, > —- solar abundance 


I 
4 
metals I 
~ —— solar abundance. 
H 100 


The second abundance has an uncertainty of at least a factor 4 because of the 
uncertainties in Alog W alone. As remarked above, these metal abundances are 
likely to be too high (A too small) because of the neglect of the contributions of 
the metal lines to the continuous opacity. These values may be compared with 
those found by Helfer, Wallerstein and Greenstein (9) for the giant stars in 
NGC 6341 (our Group C) and NGC 6205 (our Group A) using the curve of 
growth method. I am indebted to them for providing this information before 
publication. 


(Group A) NGC 6205 


clusters of Group C, 


metals 
H 


metals 
H 


They remark that these values could easily be in error by a factor of 4 either way 
but are unlikely to be in error by a factor as great as 10. Considering this, and 
the approximate natuce of the present analysis, the agreement is as satisfactory 
as could be expected. We conclude that the metal/hydrogen ratio is likely to 
lie in the range o-o1—0-001 the solar value for clusters of Group C while for 
Group A clusters it is relatively indeterminate but is likely to be of the order o-1 
of the solar value. ‘The two chief difficulties in this analysis are finding the source 
of the continuous opacity and fixing the excitation temperature. 

We have hitherto considered strong lines. A weak line formed on the 
Doppler part of the curve of growth will have an equivalent width 


I 
have — solar abundance. 
20 


(Group C) NGC 6341 have = solar abundance. 


We 2. (21) 


Ky 
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In this case, for no contribution to the electron pressure from hydrogen and with 
the metals largely ionized, the line weakening is given by: 

Alog W=log (+) (22) 
and a much larger change in line strength occurs for a given abundance change. 
It may be noted that Schwarzschild, Spitzer and Wildt (27) assume that the metal 
line strength is proportional to the partial pressure of atoms producing the line 
divided by the partial pressure of the particles producing the general opacity. 
It follows therefore that their arguments apply only to weak lines. 

It is known (14) that a (U—B)/(B—V) plot for the spectroscopically nearly 
normal giants in NGC 5272 shows an ultra-violet excess which has been inter- 
preted as being due to smaller blanketing in the ultra-violet due to a metal 
deficiency. This appears quite likely since although the strong lines are 
relatively unaffected, a metal abundance deficiency by a factor 10 will reduce 
the equivalent widths of the faint lines proportionally and a significant change 
in blanketing should occur. 

Finally consider the case of the subgiant branch of a globular cluster and 
take the colour index C.I.= +0°64 as representative of this group of stars in 
NGC 5272. Assuming 7T,=5495°, 7;=5040° and M,,=+1-2, we may 
solve equations (13), (19) and (20) for Alog W with different assumed metal/ 
hydrogen ratios. We then have: 

A-1=4x 1075 y=1'41 x 10% p=O75 

107% 1°95 x 107° 0°44 

5x 107% 2°04 x 10~* 0°32 
where p is the factor by which a strong metal line is weakened relative to a star of 
the same colour and luminosity but with solar metal abundance (A~'= 10-*). 
It will be noticed that the weakening becomes large when y > A~*; that is when 
the hydrogen supplies more electrons than the metals. It is thus seen that a 
metal deficiency which will produce little change in the spectra of the giant stars 
(as in group A clusters) may produce significant weakening of the metal lines in 
the subgiants. Since the subgiants contribute a considerable fraction of the 
integrated cluster light, we shall in the next section try to confirm this by con- 

sidering the integrated spectra of clusters. 


3. The interpretation of the integrated spectra 

3.1. The integrated spectrum of NGC 5272.—Let us consider the formation 
of an absorption line in the integrated spectrum of a cluster. The cluster may 
be considered to consist of n groups of stars so that for the ith group the equivalent 
width of a given line is the same (W,*) for all stars in the group. We shall in 
fact assume that all stars in the group have the same spectrum. Then if a, is the 
fraction of the total light of the cluster emitted by the ith group in photographic 
wave-lengths, the observed equivalent width of the line in the integrated spectrum 


will be 
> «We. (23) 


i=n 
At present this calculation can only be attempted directly for one cluster 
(NGC 5272) for which Sandage*(24) has found the luminosity function. For 
each 0-2 range of absolute photographic magnitude M,, to M,, +0-2 corrected 
according to Johnson and Sandage (14), the photographic luminosity function 
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$(M,,) is given by Sandage. The total contribution in photographic light from 
M,, to M,, +0°2 thus equals 


&(Myyy Myy 402) = [(M,,)10-°o0]) ¥ a (24) 


i=n 

where the summation of « on the r.h.s. of equation (22) is over the whole range of 
magnitudes of stars in the cluster. The plot of « against M,, (cf. Fig. 1 in (24)) 
shows a sharp peak around M,,=o which is produced by the stars of the hori- 
zontal branch. It is possible to assume that the contribution to « from the giant 
branch varies smoothly in this region and so the different contributions of the giant 
and horizontal branches may be separated in this range (—0-4< M,,< +08). 
The colour magnitude diagram (14) allows us to assign a colour to each group of 
stars within it and using Johnson and Morgan’s data (13), each group may thus 
be given a spectral type. As mentioned in Section 2.3 an approximate equivalent 
width (W,) may be assigned to each line for a given spectral type and therefore 
for each group from the data of Pannekoek (22), Wright (31) and Hiltner and 
Williams (10) if the stars have normal spectra. For abnormal spectra, a factor p 
must be applied which is appropriate to the given line. We can then find 
W,* =pW, and knowing «; from (24) find the equivalent width of the line in the 
integrated spectrum from (23). The equivalent width may then be converted 
back to a spectral type corresponding to the absolute intensity of the line for 
comparison with the observed spectrum. 


Tas.e III 


The contributions to the equivalent widths of hydrogen and metal lines in the integrated 
spectrum of NGC 5272 when spectra of stars are normal for colours 


Contribution to Contribution to Contribution to 
Group of Stars photographic equivalent equivalent width 
light of cluster width of H lines of metal lines 
per cent per cent per cent 

Giant branch B-V > +0°8 20°8 5°8 36°3 

Subgiant branch and main 

sequence for B~V < +0'8 50°2 29°! 52°0 

Horizontal branch 27°1 64°7 3°9 

Main sequence 

B-V>+08 1°9 0'4 7°8 

Table III] gives in condensed form the results of this calculation for NGC 5272, 
assuming that the spectra are normal, i.e. that they correspond to the colours. 
In the actual calculation the cluster was divided into twenty-one groups of stars 
but the condensed table with four groups is sufficient to show where the main 
contributions to the lines in the integrated spectra originate. Equivalent 
widths of different metal lines will, of course, show a different dependence on 
spectral type: the present results refer to the strong Fe 1 lines 14325 and A 4383 
which are representative of those used to obtain Radcliffe and McDonald (Fe) 
types. 

We consider two extreme cases. In the first the spectra correspond to the 
colours throughout the cluster and there is therefore no line weakening. In this 
case the (H) type for NGC 5272 is F7 and the (Fe) type is G2. In the second 
case, we assume that stars with (B— V)> +0°8 have normal spectra while those 
with (B— V)< +0°8 (i.e. the subgiants, horizontal branch and upper part of the 
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main sequence) have normal hydrogen lines but no metal lines at all. In this 
case the (Fe) type is F2 while the hydrogen type is of course unchanged (F7). We 
may compare these types with those observed by Morgan which are (H) = F8 and 
(Fe)=F1. It is seen that the observed spectrum agrees well with the case in 
which the metal lines in the subgiants are very weakened. This was the prediction 
for these subgiants (given at the end of the previous section) if the metal/hydrogen 
ratio was less than about o-1 of the solar value. We conclude therefore from our 
discussion of the integrated spectrum of a cluster of Group A that the metal/ 
hydrogen ratio is <o-1 of the solar value and not higher as would have been 
possible from a consideration of the giant stars alone. 

3.2. The integrated spectra of other clusters —Although NGC 5272 is the 
only cluster for which a full luminosity function is available, one can attempt some 
interpretation for other clusters using the data on the luminosity function of the 
giant and horizontal branches of four clusters given by Arp (2, Tables III and IV). 
It is assumed that the relative fractions of light coming from the giant, subgiant 
and main sequences are the same in all clusters and that only the contribution 
from the horizontal branch changes. The contributions to the total light from 
the horizontal branch of these clusters relative to that in NGC 5272 and from 
the rest of the cluster relative to that in NGC 5272 can be calculated and are 
given in Table IV. The cluster NGC 6356 of late integrated type is also 
considered. According to provisional unpublished data,* this cluster has only 
a weak stub on the red side of the normal position of the horizontal branch, 
ie. in the region (B—V)=+0-4. It has therefore been assigned a horizontal 
branch of the same strength as that in NGC 6205 (whose horizontal branch is also 
weak) at (B—V)= +0°4. 

Tas.e IV 
Calculated and observed integrated spectral types for five globular clusters 
NGC 5904 6205 7089 7078 
Spectral group (giant stars) A A B Cc 
Photographic contribution to 0°78 0°37 o"50 0°76 
horizontal branch relative 
to NGC 5272 
Photographic contribution to 
giant branch relative to 1°08 1°23 r18 1°09 1°23 
NGC 5272 
(H) type: 
calculated Fs-6 F8-—9 F8 F8 Gi 
observed F8(M) F8(M)_ F8(R) F6-8(MR) G2(M) 
(Fe) type: 
calculated Fs F6 Ag~Fo A4 Gs 
observed Fo(M) Fo(M)_ Fo(R) (0) (MR) Gs (M) 
Observed type M=McDonald. R= Radcliffe 


Spectral types have been calculated for the hydrogen and metal lines in five 
clusters by the method described in Section 3.1 under the following assumptions : 
(1) For NGC 6356, the spectra correspond to the colours of the stars, i.e. the 
abundances are sufficiently normal for no metal line weakening to occur. 
(2) For Group A clusters NGC 5904 and 6205, the spectra correspond to the 
colours of the stars for (B— V)> +0-8 but metal lines are absent in stars to the 
blue of this limit. 


*I am indebted to Dr A. R. Sandage for kindly telling me this in a private discussion. 
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(3) For the Group B cluster NGC 7089, the metal lines have 60 per cent of the 
normal strength for the star’s colour for (B— V)> +0°8 and are absent to the 
blue of this limit. 

(4) For the Group C cluster NGC 7078, the metal lines have 30 per cent of the 
normal strength for the star’s colour for (B— V)> +0°8 and are absent to the 
blue of this limit. 

Table IV gives the calculated and observed (H) and (Fe) types for these five 
clusters. The agreement between observation and calculation must be considered 
as good as could be expected considering the drastic simplifying assumptions 
which have been made. The greatest discrepancy is for the (Fe) types of 
NGC 6205 and 5904; it may be noted that these clusters belong to the group of 
intermediate (CH/Hy) type in which there appear to be systematic differences 
between Radcliffe and McDonald observations. 


Spectral group of cluster giants 
A B | a 


(CH/Hy) (Fe) \(CH/Hy) (Fe) |(CH/Hy) (Fe) 





= a) @) 
w Oo Ww 


Integrated type 


i] 











Fic. 5.—The integrated spectral types (CH|Hy) and (Fe) associated with the spectral groups 
of the cluster giants A, B and C. 


Fig. 5 shows the relations between the integrated (CH/Hy) and (Fe) types 
and the spectral group (A, B or C) of the cluster giant stars. The clusters 
containing giants of Group A are seen to have a wide range in the (Fe) type 
of their integrated spectra. This may be understood since we have seen that 
a wide range of metal abundance may produce practically normal spectra for the 
brighter giant stars. Such a range of metal abundance will however produce quite 
different metal line strengths in the subgiants which contribute at least half the 
photographic light of the cluster and so have a considerable influence on the 
integrated spectrum. Other differences in the integrated spectrum are likely to 
originate in the different relative populations of the horizontal branches of different 
clusters; it does not seem as though it is necessary to assume any large differences 
in the luminosity functions of the giant, subgiant and main sequences in different 
clusters apart from a scale factor. 

The present calculations contain many simplifications which later work 
should remove. To do this, however, more cluster luminosity functions together 
with photometric and spectrophotometric data on both clusters and normal disk 
stars will be needed before an accurate understanding of the composition of 
globular cluster stars is possible. It is not unlikely, for example, that the relative 
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abundances of different metals or their abundance relative to carbon or nitrogen 
is variable but investigations of this point are beyond the scope of this paper. 

We conclude, however, that the metal/hydrogen ratios relative to the Sun 
Ag "'/Ag™ derived here should be of the correct order of magnitude; they are 
summarized for different kinds of cluster below: 


Spectral group of Integrated type iB ae 
cluster giant stars (CH/Hy) Ag Ao 
A Gs-Go >1/10 
Go-Fs5 <1/10 
B F7-F2 > 1/100 
Cc Fs -F2 <1/100 


4. Correlations between spectral type and other parameters 

4.1. Spectral type and distance from the galactic plane.—Mayall (19) noticed 
originally that the clusters of later spectral type occurred near the galactic centre. 
Baade (3) has recently shown that Mayall’s spectral types show a correlation 
with the distance (Z) of the cluster from the galactic plane. He suggests that 


TABLE V 


The spectral type (CH/Hy) of integrated spectrum, distance from the galactic plane (Zp 
and distance from galactic centre projected onto the galactic plane (R,) in kpc. Assumed 
distance of Sun from galactic centre equals 8:2 kpc 


NGC (CH/Hy) |[2| R, NGC (CH/Hy) 
104 G 4°1 6°8 6333 F3-4 
362 12°6 8-4 6341 F2 
1261 23 19°9 6356 
1851 9°3 19°0 6388 
1904 6:2 18-2 6397 
2298 ‘ 33 6402 
2419 69 6440 
2808 ; 10°! 6441 

4147 10°9 6541 

4590 ; 8-7 6544 
5024 . 48 6584 
5139 , 6°4 6624 
5272 6°1 6626 
5286 8-4 6637 
5634 78 6638 
5694 

5824 36 

5904 

5927 

5986 

6093 

6139 

6171 

6205 

6218 

6229 

6254 

6266 

6273 

6284 

6293 

6304 
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clusters may be divided into two groups: those with integrated spectral type 
G3-5 which have a strong concentration to the plane and those of A5—Go which 
are a halo group which show very little concentration to the plane. 

Table V gives the integrated spectral type (CH/Hy) and the galactocentric 
distance projected onto the galactic plane (R,) for 63 clusters. The clusters are 
divided for convenience into three groups according to spectral type: G5—Go, 
Fg-F6 and Fs5-F2. As we have seen in Section 3, these groups are roughly in 
order of decreasing metal/hydrogen ratio. Fig. 6 shows the distribution of |Z] 
and R, for the three groups and shows the decrease in |Z| as the spectral type gets 
later. On the simplifying assumption that Z has a normal distribution about the 
galactic plane, we may find the variance }Z*/n—1 for each group. Snedecor’s 
variance ratio test then shows that the probability that any two of the observed 
variances could have come by chance from the same parent population is less than 
one in one hundred. We should consider whether the observed relation could be 
produced by a variation in the absolute magnitude of the RR Lyrae stars (upon 
which the distances are based) with integrated spectral type. ‘This is unlikely 
since Baade (/oc. cit.) has shown that a similar relation holds between spectral 
type and galactic latitude. Moreover, if the pulsation hypothesis of the RR Lyrae 
stars is correct, the assumed absolute magnitudes of these stars in clusters of 
earliest type are too faint and so the distances of these clusters are too small; a 
correction for this effect will therefore accentuate the spectral type-|Z| relation 
found above. 

It is seen in Fig. 6 that the ratio mean R,/mean |Z| increases as the integrated 
spectral type gets later which is an indication that the clusters of later type form a 
more flattened system than those of earlier type as distinct from a spherical system 
of smaller mean radius. This result is only tentative since there are likely to 
be many clusters (particularly of later type) which are unobserved because of 
obscuration. Schmidt (26), however, attempts to allow for this observational 
selection by effectively considering only high latitude clusters and also finds a 
systematic decrease in the flattening of the cluster system outwards. 

Whilst the largest difference in distribution undoubtedly exists between the 
groups of type G5—Go and those of F5—F2, we do not think that it can be assumed 
(as by Baade) that the intermediate group is merely a mixture of the two extreme 
groups produced by inaccurate spectral types. These inaccuracies must not be 
underestimated but the present spectra of both stars and integrated cluster light 
do show the presence of intermediate spectral types; so that some intermediate 
spatial distribution of these clusters is not unlikely. Observational selection 
undoubtedly influences the observed frequency distribution of clusters with 
spectral type. This may be roughly corrected by assuming (after Baade, loc. cit.) 
that the number of clusters with |b|<19° should be tripled. We then have: 


Range in spectral type 


Gs5-Go Fo9-F6 F5-F2 
Observed frequency 17 26 20 
Corrected frequency 45 56 32 











After correction, the intermediate group is still the largest which lends support 
to the hypothesis that a continuous range of spectral types exists and not two groups. 
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4.2. Spectral type and the RR Lyrae stars.—Arp (2) suggested that a correlation 
exists between the giant star spectra and the mean period of the RR Lyrae stars 
of type ab in a cluster: the giant stars of more normal spectral type (Group A) 
being associated with the variables of shorter mean period. Fig. 7 shows a plot 
of mean period against giant star spectral group (A, B or C) and also against 
integrated spectral type. The spread in periods among type ab variables in a 
cluster is of the order of 0-1 days or more so that only clusters containing ten or 
more variables to define a mean period have been used. The present meagre 
data give some support to Arp’s conclusion but clearly more data are required. 
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Fic. 6.—The distance from the galactic plane (|Z|) and the galactocentric distance projected 
on the plane (R,) in kpc as a function of the integrated spectral type (CH|Hy). 


























Iwanowska (12) finds a relation between the period of field RR Lyrae stars 
and their spectral type so a correlation is not unlikely between the periods of 
cluster RR Lyrae stars and the spectra of stars in the same cluster. The pul- 
sation hypothesis (15, 25) predicts that the absolute magnitude of these variables 
should be a function of both colour and period. Thus if the absolute magnitude 
of the horizontal branch is related to the spectral group of the cluster giants, the 
mean period of the RR Lyrae stars should be also if their mean colour is constant. 
Arp’s data (loc. cit.) shows that the gradient of population of non-variable stars 


39* 
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across the RR Lyrae gap is not the same for all clusters and so the mean colour of 
the variables is probably not constant*. This may partly explain the scatter in 
Fig. 7, but actual colours of RR Lyrae stars in clusters are needed to make this 
certain. 
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Fic. 7.—The spectral group of the cluster giants and the integrated spectral type (CH/Hy) 
plotted against the mean period of the RR Lyrae variables of type ab in the same cluster. 


Arp also noted that the number of RR Lyrae stars in a cluster depended upon 
the population gradient of non-variable stars across the variable star gap. We may 
express this gradient as m,—m,/n,+m,, when #, and n, are the number of non- 
variable stars on the blue and red side of the gap respectively (Arp, /oc. cit., Table 
IV). The gradient is given in Table VI together with N, the number of RR Lyrae 
stars in the cluster normalized to the same total stellar population as in NGC 5272 
using Christie’s integrated photographic magnitude as an approximate measure of 
the total population. Table VI also shows that the relative number of variables 
is large when the mean colour of the non-variable stars in the horizontal branch 


B-— V lies in or near the zone of instability (0-17 < B— V <0-39) and is few when 
(B—V) is well outside these limits. Neither (B-—V) nor the gradient 
Nn, — Nq/n, +N, however, show any correlation with the spectral group of the cluster 
giants. 

To summarize: the present classification of the spectra of cluster giants agrees 
with that of Deutsch. The classification of the integrated spectra in general 
confirms that of Mayall and Morgan at lower dispersion although small systematic 


*I am indebted to Dr A. R. Sandage for pointing this out to me. 
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differences are present. An analysis of the spectra shows that the metal/hydrogen 
ratio probably ranges from less than o-o1 to greater than o-1 of the solar value. 
Clusters with the lowest metal abundance show an extended spherical galactic 
distribution, and the greater the metal abundance, the more the clusters are con- 
centrated to the galactic plane. More photometric and spectrophotometric data 
of both cluster giants and subgiants are needed to find more accurate metal 
abundances. Also more data on cluster RR Lyrae stars are needed to confirm the 
suspected correlation between their mean period and cluster spectral type. 


TasB.e VI 
Density gradient (n,—n,/n,+n,) and mean colour (B—V) of the non-variable stars in the 
horizontal branch compared with normalized number of RR Lyrae stars (N) 
and spectral group of cluster giants 
m,—n ciiacegiiaa Spectral group 
ps thes} B-V N of cluster 
giants 


NGC 
my +N, 
6205 1°00 — 0°02 4°5 
7089 0°90 0°00 9°2 
6254 0°82 +0°09 ° 
6341 0°63 +0°08 20°6 
4147 0°58 +o'1! 160 
7078 0°47 +017 55°7 
5904 0°45 +0°16 113 
5272 0°07 +0°24 169 
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GLOBULAR CLUSTERS, III. AN ANALYSIS OF THE CLUSTER 
RADIAL VELOCITIES 


T. D. Kinman 
(Communicated by the Radcliffe Observer) 
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Summary 


The solar motion found from the radial velocities of 70 globular clusters 
is 167 +30 km/s with an apex which does not deviate from the direction of 
the Sun’s galactic rotation. A small and scarcely significant apparent 
contraction of the system is found which could be produced by systematic 
errors in the velocities found from low dispersion spectra. The radial 
velocities provide no conclusive evidence for differential motion in the cluster 
system either as a function of galactocentric distance or of integrated spectral 
type: large errors are inherent in these solutions, however, because of the 
small samples of clusters available, their unfavourable galactic distribution and 
large peculiar motions. 

It is shown that field RR Lyrae stars of the same period range as those 
occurring in globular clusters have the same solar motion and velocity 
dispersion as the clusters; it is suggested that these stars and the clusters have 
a common origin. 

A scatter diagram of the ratio of the radial velocity to the circular velocity 
against the angle subtended at the cluster by the Sun and the galactic centre 
confirms von Hoerner’s conclusion that the cluster orbits are mainly highly 
eccentric. A high eccentricity (e=0-8) is also found by considering the mean 
galactic rotation of the cluster system. 

A simple model is used to determine the distribution of mass in the Galaxy 
from the cluster motions. The results are consistent with Schmidt’s model 
except possibly at galactocentric distances beyond 12 kpc. 

The mean angular momentum per unit mass of the cluster system and of 
the rest of the Galaxy is shown to be the same. The bearing of this on the 
origin of the cluster system is briefly discussed. 





Introduction.—Recent analyses of the radial velocities of globular clusters 
have been based on Mayall’s (11) velocities of 50 clusters. After correction for 
local solar motion, Mayall found a solar motion with respect to the cluster system 
of 175 +25 km/s from the mean of four solutions. He considered that this result 
was probably influenced by the absence of velocities for southern clusters and 
that the probable value for the solar motion was 200+ 25 km/s. Since he found 
from other evidence (motions of the extragalactic nebulae and the Oort galactic 
rotation constants) that the Sun’s galactic rotation velocity was 280-300km/s, he 
suggested that the cluster system must have a rotational velocity of the order of 
80-100 km/s. However, he was unable to find any satisfactory evidence for a 
differential or general rotation of the clusters from the peculiar velocities of the 
clusters themselves. Later, Parenago (17) and also Perek (18) claimed to have 
found differential rotation in the cluster system but Schmidt (21), using the 
same observations, has been unable to find any evidence for this, 

(on the assumption of a solar galactic rotation of 216 km/s) he finds a mean 
rotational velocity of 80 km/s for the cluster system. 
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In the first paper in this series (8) radial velocities were given for 30 southern 
globular clusters, 10 of which had been observed previously. Radial velocities 
are thus available for 70 globular clusters (or roughly two-thirds of the total 
numbers of clusters known) and the effects of observational selection should 
be much reduced. It is therefore appropriate to attempt a fresh analysis to try 
to resolve the conflicting results found previously. In the first section we consider 
various solutions for the solar motion with respect to the cluster system as a 
whole, while in the second section the motion of different groups of clusters is 
discussed in an attempt to detect differential motion. In the third section the 
motions of the globular clusters are compared with those of some other high 
velocity galactic objects and the galactic orbits of the globular clusters are 
considered in a fourth section. Ina fifth section an attempt is made to determine 
the mass of the Galaxy from globular cluster motions while the brief concluding 
section deals with the origin of the system of globular clusters. 

1. The solar motion of the system of globular clusters—We may use the radial 
velocities of the clusters to determine the motion of the clusters as a whole 
relative to the Sun. It is more useful, however, to refer the cluster motion to 
a coordinate system which has the same motion as the nearby stars. The 
heliocentric cluster radial velocities were therefore corrected for the local solar 
motion relative to these stars. ‘The standard local solar motion of 20 km/s 
towards «= 18" and §= + 30° is of sufficient accuracy because of the dispersion 
and accuracy of the cluster velocities. This corrected velocity V, for a cluster 
with equatorial coordinates («,5) is related to the component cluster velocities 
X, Y and Z by 

X cosacos6+ Ysinacoss +Zsind=V, (1) 


where X, Y and Z are directed towards (0", 0°), (6%, 0°) and (+ 90°) respectively. 
1950°0 coordinates have been used throughout. A least squares solution for the 
70 clusters then gives the following solution : 

X= —73 + 33 km/s; 

Y= +93 +22 km/s; 

Z=—119+28 km/s. 


Unless otherwise stated, the errors quoted in this paper are r.m.s. errors. The 
total cluster motion (U) relative to the nearby stars is thus (X?+ Y?+Z*)** 
which equals 168+27 km/s. The apex of this cluster motion is 


a=(21527™)+(26™); S=(+45° 08’) + (6°08’). 


Converting to epoch 1900-0 and using the Lund galactic coordinate tables (13), 


which are of sufficient accuracy for this purpose, the galactic coordinates of the 
apex are: 


I= 5§7°84+7°5; b=—4°647°6. 


The apex of cluster motion therefore shows no deviation from the direction 
in the galactic plane at right angles to the galactic centre (/= 58°, b=0°) within 
the errors of the solution. No attempt has been made to weight the cluster 
velocities in this solution or later ones, since the errors in the component velocities 
are largely due to the actual dispersion in the cluster velocities rather than to 
observational errors. Mayall, in a similar solution for 50 clusters (11, p. 312, 
solution 6), found U=168+34km/s and J=51°-3+10°1, b= +8°-2+13°0 
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for the apex (Harvard = Lund pole); Mayall’s probable errors have been changed 
to r.m.s. errors for comparison purposes. It is seen that an increase in the number 
of clusters together with a more uniform distribution over the sky has produced no 
change in the value of the solar motion although the errors are slightly reduced. 
In all the following solutions we shall assume that the apex is in fact at right 
angles to the galactic centre and in the galactic plane (= 58°, b=0°). Let us 
also introduce a K term of motion of the clusters relative to the Sun. The equation 
of condition for a cluster with galactic coordinates (/, 6) is then 


UcosA+K=V, (2) 


where cos A =sin (/—328°)cosb. A least squares solution for the 70 clusters 
then gives: 


U=172+31 km/s; = —16+15 km/s. 


We find, in agreement with Mayall, that the K term relative to the Sun is of 
the same order as its error. It is therefore omitted in the next solution, leaving 
as the equation of condition 

UcosA=V, (3) 


for which the least squares solution is 
U=167 + 30 km/s. 
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Fic. 1.—The ordinate is the cluster radial velocity (corrected for local solar motion) in km|s and the 
abscissa is the cosine of the apical distance. New velocities are shown by filled circles and previously 
known velocities by open circles. Previously known velocities corrected by new observations are 
shown by half-filled circles. The solar motion solution U=167= + 30 km/s is shown by the line. 
The apical distances of clusters with unknown velocities are shown by the filled squares along the 


abscissa axis. 
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Following Mayall (/oc. cit., Fig. 3), this solution may be represented as a straight 
line in Fig. 1 which is a plot of the cluster velocity against the cosine of the apical 
distance (A). Clusters with previously known velocities are shown by open 
circles, clusters with new Radcliffe velocities are shown by filled circles and 
clusters with both previously known and Radcliffe velocities are shown by 
half-filled circles. The apical distances of the 35 clusters in Miss Sawyer’s 
catalogue (20) with unknown velocities are indicated along the abscissa axis: 
these are mainly faint clusters near the galactic centre and only a few have a 
sufficiently large apical distance to affect a solar motion solution. 








N 


Fic. 2.—The Sun (S) in the galactic plane has a circular velocity e, about the galactic centre (O). 
A cluster (C) has three component velocities: ©, about the galactic axis of rotation OZ, Rg along 
OC and @ at right angles to both @,. and Rg. 


The relative motion between the cluster system and the nearby stars may 
be expressed in terms of the galactic motions of the two systems. In Fig. 2, the 
Sun (S) is assumed to be in the galactic plane at a distance R, from the galactic 
centre (O). A globular cluster with galactic coordinates (/,5) is distant R from 
the Sun and R, from the centre of the Galaxy. R, is the distance R, projected 
onto the galactic plane. The circular velocity of the nearby stars is ©,. The 
velocity components of the cluster are ©, around the axis of rotation OZ of the 
Galaxy, R, radially outward from the centre of the Galaxy and a tangential 
component @ at right angles to these two. The observed radial velocity 
corrected for local solar motion is then given by 


a , R R—R, cosAcosb 
V.= (90,-0,.—)c 2 
f ( ® 0, By cos A+ Ry ( R, ) 





+ ze (®- F (Reosb— Rycos)) cosb) (4) 
g c 


where cos A =sinAcosb as before and A=1—1,, with 1, =328°. 
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It seems reasonable to assume that with a sufficiently large number of clusters, 
the tangential component ® will average zero and this term may therefore be 
omitted. We now see that equation (3) above implies that for all clusters: 

(2) }R,=0, i.e. that there is no expansion or contraction of the cluster 

system ; 

(6) ©.(Ro/R,) =constant, i.e. that the cluster system rotates as a solid body 

with angular velocity w,=0,/R,. 
Equation (3) may thus be written as 


V,= Ry(w,—,) cos A, (5) 


where w, is the angular velocity of the nearby stars. On these assumptions and 
taking R, = 8-2 kpc, the difference in angular velocity between the cluster system 
and the nearby stars is 20-4 + 3°7 km/s kpc. 

Let us now remove assumption (a) and assume that the cluster system 
expands so that R,=«R, where ¢ is constant. The equation of condition is then 


V..=(w,—,)R, cos A + «(R— Ry, cos A cos b) (6) 
and a least squares solution for 70 clusters gives (in km/s kpc) 
(w,—w,) =20°5 +36; «= —1-48+0°85. 


The velocity of contraction of the system is only slightly larger than its r.m.s, 
error. To investigate the reality of this result, we divide the clusters arbitrarily 
into two equal groups by the NGC numbers and again solve equation (6) to 
obtain (in km/s kpc) 
NGC 104-6284 
(w,—w,)=20°54+51; €= —1r2rt+1°08. 
NGC 6293-7099 
(w,—@,)=16-947°4; €= —2°93 $2°12. 
Both solutions indicate a contraction. The clusters NGC 2419 and NGC 7006 
are very distant from the galactic centre and they will have considerable weight 
in the solution for e. Their velocities are based on low dispersion spectra however 
and thus are of low weight; they may therefore be omitted and the solution for 
the remaining 68 clusters gives (in km/s kpc) 


(w,—w,)=20°0+ 38; e= —1°764 1°12, 


which again indicates a contraction. If these values of « are indeed due to a 
contraction then the cluster system will collapse in the order of 10° years. This 
implies a much more rapid rate of kinetic energy loss than is likely from Spitzer 
and Schwarzschild’s (22) relaxation times of such stars in their encounters with 
interstellar gas clouds. A more likely explanation is that it is caused by a 
systematic error of the order of 10 km/s such as may be present in the low 
dispersion spectra. The calcium H line which is blended with He has a large 
weight in these spectra. Its true wave-length is 0-35 A bluer for clusters of 
late integrated type than for those of early type (Paper I (8), Table VI), which 
corresponds to a velocity difference of 25 km/s. No allowance for this effect 
is mentioned in Mayall’s paper; if a mean wave-length was in fact used for all 
clusters, then the clusters of late type will show a spurious velocity of recession 
and those of early type will show a spurious velocity of approach. Since the 
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late type clusters occur predominantly near the galactic centre (Paper II (9)), 
this will produce a spurious contraction of the whole system. Ten clusters of 
early type observed by Mayall do show a mean velocity of approach of 4 + 8km/s 
compared with the Radcliffe observations (Paper I (8)), but since none of these 
Mayall velocities is wholly based on low dispersion spectra, the effect may be 
reduced in this case. 

It is thus concluded that the observed contraction of the cluster system could 
be produced by systematic errors in the velocities such as could be produced by 
errors in the standard wavelengths in the low dispersion spectra. Only the 


first term on the r.h.s. in equations (4) and (6) will therefore be retained in future 
solutions. 


2. Differential motion in the system of globular clusters.—In the previous section 
it was assumed that the cluster system rotates as a solid body; the linear 
rotational velocity ©, of a cluster is thus proportional to its distance (R,) from the 
galactic axis of rotation. In the simple case of circular orbits about a central 
point mass, this circular velocity should be inversely proportional to the square 
root of the radius of the orbit. We should therefore expect ©, to decrease with 
galactocentric distance. Since, however, the eccentricity and major axes of the 
cluster orbits are unknown and the central point mass approximation is only 
valid for the clusters which are most distant from the galactic centre, it is inadvisable 
(as Schmidt (21) points out) to try to fit the observed velocities to an analytical 
form. We therefore proceed empirically and divide the cluster system into two 
groups, each rotating as a solid body: 51 clusters for which the projected galacto- 
centric distance R,<gkpe and 19 clusters for which R,>g kpc, and solve 
equations (3) or (5) for both groups. The division at 9 kpe was chosen to give 
both solutions nearly equal weight; many more clusters were therefore needed 
in the first group because cos A is small for clusters near the galactic centre. 
The solution of equation (5), assuming R,=8-2kpc, gives the following 
differences in angular velocity between the two groups of clusters and the nearby 
stars : 


R,<9 kpe R.>9 kpe 
(w, — w,) = 18-2 + 5:0 km/s kpc (w, —w,) = 22°4 + 6-2 km/s kpe 


while the difference in linear rotational velocity U=(0,—0©,R,/R.) found from 
equation (3) is: 


R.<9 kpe R.>9 kpe 
U=149+ 41 km/s U=183 +51 km/s 


The difference in U between the two groups is thus 34 +65 km/s. Although the 
difference in rotational motion between the two groups is in the expected sense, 
i.e. the outer group have the larger solar motion and thus rotate more slowly, 
the errors of the solution are too large for the result to be statistically significant. 
There seems little point in attempting further solutions with a finer subdivision 
of R,, partly because of the uncertainties in the values of R, and partly because 
decreasing the number of clusters in each group would increase the already large 
errors in the solution for each group. We conclude that, although differential 
motion is very likely to exist, the errors produced by the peculiar motions of the 
clusters preclude its detection. 
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In Paper II (9) it was shown that the galactic distribution of the globular 
clusters depends on their integrated spectral type; those of earliest type have an 
extended spherical distribution and those of latest type have a distribution which 
is flattened to the galactic plane. In this case, we might expect on dynamical 
grounds that the clusters of later type should have a larger rotational velocity 
and hence smaller solar motion than those of earlier type. To test this, the clusters 
were divided into two groups according to their spectral type: Group I consisting 
of 17 clusters of spectral type Gs~Go and Group II consisting of 46 clusters of 
spectral type Fg-F2. The division was made to give the largest difference in 
spatial distribution between the two groups. Least squares solutions for U 
and (w,—w,) on the assumption R, = 8-2 kpc then give: 


Group I Group II 
Gs5-—Go Fg9-F2 
(w,—«,) g°8 + 10:0 km/s kpc 19°8 + 4:4 km/s kpc 
U 80 + 82 km/s 162 + 36 km/s 


The difference in U between the two groups = 82 +90 km/s. It is seen that the 
difference in the solar motions of the two groups is in the expected sense, i.e. the 
groups of later spectral type have the smaller solar motion and hence the larger 
galactic rotation, but again the errors of the solution are too large for the result 
to have statistical significance. The large errors in the solution for Group I are 
caused by the majority of the clusters of this group lying near the galactic centre 
and so having small values of cos A. Thus while differential motions of the kind 
indicated between clusters of different spectral type are very likely to exist, they 
are largely masked by the errors introduced by the peculiar motions of the clusters. 

3. A comparison of the solar motion of globular clusters with other galactic 
objects.—Large radial velocities are found for several other classes of galactic 
objects besides the globular clusters and together with the globular clusters these 
form the galactic halo. However it is necessary to restrict the discussion to such 
classes of objects which may be distinguished as belonging to the halo by some 
physical parameter (such as spectrum or light variation) other than high velocity. 
Clearly, the analysis of the velocities of a group of objects selected on the grounds 
of high velocity will yield a spurious solar motion. In this section we shall therefore 
deal only with the RR Lyrae stars and, rather briefly, with the subdwarfs. 

3.1. The RR Lyrae stars—Oort (14) analysed the radial velocities of 67 
RR Lyrae stars and found a solar motion of 130 + 21 km/s towards 1=53°+ 12°, 
b=+12°+12°. Since then velocities of more of these stars have become 
available so that a new analysis is possible. This was made for 129 stars whose 
velocities are either given in Wilson’s General Catalogue (27) or by Joy (7) 
and assuming an apex of /=58°,b=o0°. After correction for local solar motion, 
this gave a solar motion U of 134+20km/s. It was noticed by Struve (24) 
that the solar motion of the RR Lyrae stars increases with increasing period 
but he failed to discriminate between variables of the Bailey types ab and c. 
Such discrimination must be made, however, since the variables of Bailey type c 
of period 0-3-0°4 days are known to be associated in clusters with those of type ab 
of period o-5-o°6 days. In the present analysis, the 129 RR Lyrae stars were 
divided into five groups according to period and Bailey type and for which the 
following solar motions (U) were found: 
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Period range Mean U 
period 
of-25 — of-35 01-296 89 +29 km/s 
04-20 < P< 0-40 01-369 30+27 
01-40 < P<o04-50 01-460 105 +29 
01-50 < P<01-60 01-548 167 + 30 
01-60 < P<o0f-75 01-666 205 +45 
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Fic. 3.—(a) Solar motion (U) in km|s of field RR Lyrae stars as a function of their period in days. 
(6) r.m.s. random sight-line velocity (¢) in km/s after removal of solar motion component 
for field RR Lyrae stars as function of their period in days. 
The corresponding quantities for globular clusters are also shown. 


These results are also given in Fig. 3a, which shows the marked increase 
of the solar motion with increasing period. The mean periods of the ab type 
variables in different globular clusters lie roughly in the range 04-50-04-65 
(see Fig. 4a). It is seen that the solar motion of the RR Lyrae stars of this 
period range corresponds closely to that found for the globular clusters 
(167+30 km/s). Furthermore, the r.m.s. random sight-line velocity (c) after 
removal of the solar motion is also the same within the quoted errors (Fig. 35) 
for the field RR Lyrae stars of the period range 04-50-04-65 as for the globular 
clusters. 

There is therefore reason to believe that the RR Lyrae stars of a given period 
have essentially similar galactic motions whether they occur in globular clusters or 
in the general field. This is contrary to the findings of Oort (loc. cit.) who found 
a significant difference between the solar motions of the globular clusters and the 
RR Lyrae stars of all periods and therefore concluded that the field RR Lyrae 
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stars could not have escaped from the clusters. The present results however 
support a common origin for RR Lyrae stars of a given period and Bailey type 
in both clusters and the general field. Either the field variables result from the 
disruption of globular clusters or, more generally, their presumed origin from 
the interstellar medium occurred at the same epoch as for the clusters and with 
the same galactic distribution. 

These results suggest that the globular clusters whose RR Lyrae variables 
of type ab have the longest mean periods may be expected to have the greatest 
solar motion. They will therefore have a nearly spherical galactic distribution 
and a large mean distance (|Z]) from the plane; hence, according to Section 4.1 
of Paper II they will have the earliest spectral type. It was shown in Section 4.2 
of Paper II that the few available observations tend to confirm this. 

In the previous section, we were unable to demonstrate convincingly that a 
difference in solar motions does exist between clusters of different spectral types : 
this was due to the small number of clusters available and their unfavourable 
distribution in the sky for solar motion solutions. While it is unlikely that the 
number of spectroscopically observable clusters can be greatly increased in the 
near future, this is not the case with the RR Lyrae variables. Moreover, their 
distribution in the sky for a solar motion solution is relatively more favourable 
on account of their greater spatial density near the Sun. We must aim therefore 
to improve the solar motion of the field RR Lyrae stars of longest period and this 
should prove the best method of obtaining a lower limit to the Sun’s galactic 
rotation velocity. 

Fig. 4a shows the period—frequency function for the type ab RR Lyrae 
stars in globular clusters (taken from Sawyer’s catalogue (20)). The same 
function is shown in Fig. 45 for the type ab variables with galactic latitude 
b< +30° (taken from the General Catalogue of Variable Stars (4) but omitting 
those due to Gaposchkin (3) at the galactic centre). It is seen that variables with 
periods less than o%-45 are much more numerous in the latter group which lie 
near the galactic plane; this is in accord with the small solar motion of this 
group (cf. van den Bergh (25), Oosterhoff (16) and Kukarkin (1o)). 

It is possible that variables with periods ~ 04-45 and less may occur in globular 
clusters of late integrated spectral type (G3-5) which lie close to the galactic 
plane. Variables have been found in those clusters of this type which have been 
searched (NGC 104, 6356, 6712 and 6833 (20)) but as yet no periods have 
been published. Relatively few of such clusters are available for searching: 
no doubt some are not observed because of obscuration, It is also possible that 
such clusters were more numerous in the past but, lying near the plane, they 
have been broken up by the tidal action of interstellar clouds: a mechanism 
which has been discussed by Spitzer (23) to account for the scarcity of 
“galactic” clusters with ages > 10° years. 

3.2. The subdwarfs.—Fricke (2) has analysed the radial velocities of 
78 subdwarfs and has found a solar motion of 148 + 28 km/s towards /=75°+8° 
and b= —2°+9°. Apart from an unexplained deviation of the apex in longitude, 
this is similar to the solar motion found for all the RR Lyrae stars taken together. 
It seems likely that, as in the case of the RR Lyrae stars, we are dealing with a 
mixture of different kinds of stars; indeed it is reported (§) that Greenstein 
has found that the space velocity dispersion in the subdwarfs depends on their 
metal line weakening: the stars with the weakest lines having the greatest velocity 
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dispersion and hence presumably the greatest |Z|. This is in qualitative agreement 
with what was found for the globular clusters in Paper II: the clusters with the 
weakest metal lines (earliest integrated spectral type) have the greatest |Z]. 
By analogy with the RR Lyrae stars we might expect that globular clusters 
containing subdwarfs with a given line weakening (or metal/hydrogen ratio) 
will have the same galactic motions and hence presumably the same origin as 
subdwarfs of the same physical characteristics in the general field. The same 
should also hold for all the constituent stars of globular clusters although, apart 
from the RR Lyrae stars and subdwarfs, it may not be so easy to identify the 
field stars with the same physical characteristics as their counterparts within 
clusters. In some cases it is possible that groups of the field stars having the 
same origin as globular clusters may retain a loose dynamical association by which 
they could be identified: they would thus exist as a very loose cluster. It is 
interesting to note in this connection that Eggen* has recently found that RR Lyrae 
itself is a member of a group of stars with a common motion. 











0.4 0-6 
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Fic. 4.—Frequency of occurrence of RR Lyrae variables of Bailey type ab as a function of period 
(a) in globular clusters 
(b) in the general field for |b| < 30°. 


4. The galactic orbits of the globular clusters.—The eccentricities of the cluster 
orbits have been discussed by von Hoerner (26) and Perek (19) on the basis 
of Mayall’s velocities. Contrary to an earlier discussion by Edmondson (1), 
using only 26 clusters, they both conclude that the cluster orbits are highly 
eccentric. Von Hoerner considers the ratio ug between the radial velocity of 
a cluster (u) (corrected for the Sun’s galactic rotation and for local solar motion) 
to the circular velocity at the cluster. He assumes a point mass for the Galaxy 
at the galactic centre such that the circular velocity at the Sun is 260 km/s and 


* I am grateful to Dr Eggen for giving me this information in advance of publication. 
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R, is taken as g kpc; in this case u=uV(R,)/780, where R, is the distance 
from the galactic centre to the cluster. If # is the angle subtended at the cluster 
by the Sun and the galactic centre (see Fig. 2), then in order that a globular 
cluster should not escape from the Galaxy: 

for circular orbits, uy <sin#; 

for straight line orbits through the galactic centre, uy < +/2.cos#; 

for elliptical orbits, uy < 4/2, 
since the escape velocity on this model equals 4/2 times the circular velocity. 
Von Hoerner excludes clusters nearer to the centre than 6kpe so that the 
point-mass model of the Galaxy may be approximately valid; he then finds, 
with only two exceptions, that u%»<4+/2.cos# so that the orbits must be highly 


eccentric. 
Apart from the approximate nature of the point-mass assumption, the distances 


of the globular clusters are poorly determined (Paper I (8), Table VIII): 
probable errors of the order of 30-50 per cent being common. These errors 
will therefore introduce considerable uncertainty in the calculated values of #3, 
particularly for clusters near the galactic centre, while, since u, depends upon 
V R,, this will also be affected but toa lesser degree. Von Hoerner’s calculations 
have been repeated assuming R,=8-2 kpc and 0,=216 km/s, so that 


uVR, . 


ae 619 ° 


clusters with R, <6 kpc have again been omitted. 
The results are shown in Fig. 5 in which up is the ordinate and # is the 


abscissa; the scatter of points again favours the highly eccentric as opposed to 


circular orbits. The horizontal lines indicate probable errors in # for selected 
clusters determined from the probable errors in the distances alone (Paper I (8), 
Table VIII), while the vertical lines indicate probable errors in u) determined 
from the probable errors in the distances and radial velocities alone. 

The only cluster which appears to be necessarily escaping from the Galaxy 
is NGC 5694; the distance of this cluster is very uncertain however, and its 
velocity is based only upon low dispersion spectra so that this result must be 
treated with reserve. NGC 3201, whose velocity and distance are relatively 
well-determined, deserves special mention. It is situated 11° from the solar 
antapex of galactic motion and has the large velocity of recession (corrected for 
local solar motion) of 481 km/s. It therefore appears that for any reasonable 
value of the circular velocity of the nearby stars this cluster must have a retrograde 
orlit. 

The accuracy of this analysis is clearly not high but does show that the values 
©,=216 km/s and R,=8-2 kpc are at least compatible with few if any of the 
clusters having velocities in excess of the velocity of escape and the majority 
having highly eccentric orbits. One would not expect to find clusters with 
velocities in excess of the escape velocity. As remarked earlier, there is probably 
little energy exchange between the clusters and the rest of the Galaxy; hence 
if, as is thought, the cluster ages ~ > 10° years, any cluster whose original velocity 
was excessive would have escaped long ago. 

Schmidt (21) has pointed out that many highly elongated orbits would be 
expected a priori since the dispersion in the cluster velocities is of the same order 
as the circular velocities. Let us take a very simplified model in which a cluster 


rT) 
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describes an elliptic orbit of semi-major axis a and eccentricity e around a point 
mass M. If p is its constant angular momentum per unit mass, then we have: 


. neon os PES 
e= =) where A= ou’ 


Here R is the time averaged mean radial distance of the cluster from the mass M, 
and taking the mean value of R, for all the clusters may be put as 11-8 kpc. 
Since the orbits may be inclined at random to the galactic plane, we put p equal 
to +/2.0,.R,, where ©, is the linear rotational velocity about the galactic axis 
of rotation. We have found that 


(6, 0, R) = 167 km/s 
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Fic. 5.—The ordinate (ug) is the ratio of the cluster radial velocity after removal of the 
component due to the Sun’s galactic rotation, divided by the circular velocity at the cluster. The 
abscissa (8) is the angle subtended at the cluster by the Sun and the galactic centre. 

The solid horizontal line at uo= 4/2 is the limit for elliptical orbits. 
The dotted line (ug=sin 8) is the limit for circular orbits. 
The dashed line (uy== +/2 . cos 3) is the limit for straight-line orbits. 
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and so assuming ©,=216 km/s, Ry=8-2 kpc and R,=11°5 kpc (which is the 
mean value weighted according to cos A), we have 0,=69 km/s. The quantity 
GM is equal to ©,2R, and hence we find that the eccentricity equals 0-8, i.e. the 
orbits are highly eccentric. 

It is interesting to note that Newkirk (12) also finds highly eccentric orbits 
for the field RR Lyrae stars: out of 50 stars discussed on the basis of proper 
motions and radial velocities, only 12 were found to have eccentricities less than 
0-6. This lends further support to a common origin for the field RR Lyrae 
stars and the globular clusters. 

5. The mass of the Galaxy from globular cluster motions.—Consider a cluster 
moving in a Keplerian orbit of eccentricity e and semi-major axis a about a point 
mass galaxy of mass M. The perimeter of this orbit is 


while the period is 27a**#/V(GM). The time averaged mean galactocentric 
distance of a cluster (R,) equals a(1 + }e*), so that the mean square velocity of 
the cluster in its orbit (v*) is 

GM nl —a- oA oe y 

z (+44) (1 é* ae Pv taal? Ee 
The coefficient containing e is unity for a circular orbit, 0-88 for e=o-8 and 0-61 
for a straight line orbit; following the conclusions of the last section we shall 
take the coefficient 0-88. 

Since only the radial velocity and not the space motion of the clusters are 
known, the mean square velocity has been deduced on two assumptions. If v 
is the residual radial velocity of a cluster after removing local solar motion and 
a solar galactic motion of 216 km/s, we assume either : 

(1) the peculiar motions of the cluster are random so that 


E 


n 


vt = 


for a group of a clusters; 

(2) the peculiar motions are along the line through the cluster passing through 
the galactic centre which makes an angle # with the line joining the cluster 
to the Sun. In this case, weighting according to cos? and omitting 
clusters with cos? < 0-3, 


2 saaty ore 


> cos# 





vo? = 


for a group of n clusters. 

In order to consider approximately homogeneous groups, the clusters were 
divided into 4 groups according to spectral type: Gs5—G3, G2—Go, Fg-F6 and 
Fs5-F2. The results are shown in Fig. 6, in which the mass calculated on the 
two assumptions is plotted against R,. The open circles are for assumption (1) 
(random motion) and the filled circles are for assumption (2) (straight line motion). 
The increase in M as the mean galactocentric distance of the group increases is 


40* 
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likely to be due to the greater part of the galactic mass encompassed by the clusters 
with greater R,. For comparison, the solid curve gives the mass of the Galaxy 
within a distance R, of the centre computed from Schmidt’s model (21, Table 5). 
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Fic. 6.—The galactic Mass M(R,) within a galactocentric distance Rg (in units of 104 gm) 
as a function of Ry. 

The smooth curve from Schmidt's model 

The filled circles are masses deduced assuming that the clusters have straight line orbits through 
the galactic centre. 

The open circles are found on the iption of randomly oriented orbits. 





Within the range of Schmidt’s model (0<R,<12-3 kpc) the agreement is 
satisfactory considering the approximate nature of the present analysis and 
suggests that straight line orbits are a better approximation than random motions. 
The same solar circular velocity (216 km/s) is however used in both Schmidt’s 
model and the present analysis and so the two results are not completely 
independent. The presence of a substantial galactic density for R,>12 kpc 
suggested by the group of twenty clusters with spectral type F5—F2 on either 
assumption is surprising. Whilst it is difficult to assess just how trustworthy 
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this result may be, in view of the very simple model which has been used and also 
the uncertainties in the cluster distances and the solar circular velocity, it is 
suggested that it should be treated with reserve. 

6. The origin of the globular clusters.—We can find the mean angular momentum 
per unit mass of the cluster system (f,) about the axis of galactic rotation since 
this equals the product of the mean linear rotational velocity (@,) and the mean 
galactocentric distance projected onto the galactic plane (R,) of the clusters. 
Using values of ©, and R, which are given in Section 4, we have 


b.=794 km/s kpe. 

We may also find the corresponding quantity for the rest of the Galaxy (p,) 
using Schmidt’s model (2x). The angular momentum of an annular zone 
the Galaxy of mass m(r), circular velocity @(r) and radius r is r.m(r).O(r), 
so that the mean angular momentum per unit mass is 





by= 2° ae = 800 km/s kpe. 
The agreement between p, and , is very striking; it is better than could have 
been expected from the limited accuracy of the data available although it may 
be noted that the same solar circular velocity (216 km/s) has been used to find 
both , and #, and that an error in this velocity will produce an error in both p, 
and p, in the same sense. 

Thus, although the galactic distribution of the cluster system differs markedly 
from that of the disk population (which largely determines f,), the mean angular 
momentum per unit mass of the two systems is the same and this must also have 
been true of the media from which these systems were formed. It seems likely 
therefore that the clusters were formed from representative samples of the same 
medium as the rest of the Galaxy. This lends some support to the proposal 
(cf. Oort (15)) that the Galaxy originated in a large turbulent gas cloud. This 
cloud contracted as the turbulent velocities decayed and, since angular momentum 
was conserved, formed the present rotating disk of gas and stars. The globular 
clusters are considered to have been formed from the gas when its distribution 
was still nearly spherical and since little energy exchange can have taken place 
between the clusters and the rest of the Galaxy since then (see Section 1), their 
present distribution and motion must reflect that of the original gaseous cloud 
at the time of their formation. It has been shown (Paper II) that the clusters 
with the most extended and spherical galactic distribution have the earliest 
spectral type and hence a lower metal/hydrogen ratio than clusters with a less 
extended and flatter distribution. One possibility is that the clusters with low 
metal content were formed first when the primordial gas was most extended and 
the clusters with higher metal content later when the cloud had contracted and 
flattened somewhat. In this case, the metal content of the Galaxy must have 
increased with time—possibly by a mechanism such as that considered by Hoyle 
(6). Alternatively, the clusters may have been formed nearly simultaneously 
but the metals (perhaps in the form of dust) may have contracted to the plane 
more rapidly than the gas. This latter alternative might be rejected if it could 
be shown conclusively that the interaction between gas and dust were sufficient 
for no separation to be possible in the time required for the gas to contract to the 
plane. 
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To summarize : 

The cluster system shows a solar motion of 167 + 30 km/s with no deviation 
of the apex. Differential motion in the cluster system is expected from differences 
in galactic distribution of clusters of different spectral types; it has not been 
possible to confirm this from the radial velocities because the small samples of 
clusters available, their unfavourable galactic distribution and large peculiar 
motions produce large errors in the solutions. 

Field RR Lyrae stars of the same period range as those occurring in globular 
clusters have the same solar motion and velocity dispersion as the clusters; 
it is suggested that these stars and the clusters have a common origin. 

It is confirmed that the clusters probably have highly eccentric orbits as is 
expected from the low galactic rotation of the cluster system. An attempt to 
find the distribution of galactic mass from the peculiar motions of the clusters 
gives results which are, in general, consistent with Schmidt’s model within their 
limited accuracy. 

The mean angular momentum per unit mass of the cluster system is the 
same as the mean for the rest of the Galaxy. ‘This is in accord with the hypothesis 
that the clusters were formed in a spherical gas cloud which was in the process 
of contracting to form the galactic disk. 

Unfortunately, it is unlikely that the number of cluster radial velocities can 
be greatly increased in the near future and hence any improvement made in the 
present solar motion solution. As was pointed out in Section 3, however, the 
number of radial velocities of RR Lyrae stars of long period (>0%-6) may be 
increased (particularly if new stars are discovered) and these should give a solar 
motion which is a better approximation to the Sun’s galactic circular velocity. 
A radial velocity programme which includes such stars in the Southern Hemisphere 
is at present being undertaken at the Radcliffe Observatory. 
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Note added in proof.—Two RR Lyrae stars have been identified (A. D. 
‘Thackeray, A. J. Wesselink and M. W. Feast (M.N., in press)) as belonging to 
NGC 104 (47 Tuc) which has an integrated spectral type G3; one of type ab has 
a period 0-74 days and the other of type c a period of 0-37 days. Although this is 
a small sample, the discovery of such long periods in a cluster of late integrated 
type weakens the correlation between period and integrated spectral type discussed 
here and in Paper II. 





SOME FORMULAE CONNECTING HARVARD x AND y 
EQUATORIAL COORDINATES AND STANDARD COORDINATES 
FOR THE MAGELLANIC CLOUDS 


A. F. Wesselink 
(Communicated by the Radcliffe Observer) 


(Received 1959 March 2) 


Summary 


The paper deals with the essentially practical problem of identification of 
stars in the Magellanic Clouds. Formulae are given which convert Harvard 
x and y equatorial coordinates and standard coordinates into each other. 
Their use for the identification of objects and for the ruling of prints is 
described. In order to simplify the vast amount of future work that is likely 
to be concerned with the Magellanic Clouds, the advantages of the systems of 
standard coordinates in preference to the Harvard x, y systems are pointed out. 





Notation and definition 


x, y: the rectangular coordinates introduced by the Harvard Observatory 
in H.A., 60, IV. Each Magellanic Cloud has its own system. In the first 
instance a system is defined by a negative of the corresponding Magellanic Cloud 
on which the square net of a reticule has been photographed. The original 
plates were taken with the 24 inch Bruce Telescope at the former Harvard Station 
at Arequipa, Peru, and are of size 14x 17 inches. The reproductions published 
in H.A., 60, IV are about half this size. 

a, 5: equatorial coordinates referred to 1875-0. The epoch 18750 was 
chosen as this is the epoch of the CPD and the AGC. In the derivation of the 
formulae AGC positions have been used. For future application and use, the 
CPD positions seem to be the most useful. 

a», 59: the equatorial coordinates of the centres of the Harvard plates as used 
in the present investigation. The coordinates differ slightly from the centres 
published in H.A., 60 and are rounded off values of the latter. We have (again 
for 1875-0) 


Small Magellanic Large Magellanic 
Cloud Cloud 
t= Oh 5o™Mco8-o ay = 5226™008-0 
89= —73°15'00" 59 = — 69°00'00” 
Aa =a— a%. 

Standard coordinates.—The standard coordinates as used here comply with 
the definition in general use. The origin of the coordinates are the positions a», 
5) as given above. € and 7 are expressed in seconds of arc. The 7-axis is 
directed towards the pole of 1875-0. 

The problem of identification.—There are essentially two methods of identify- 
ing celestial objects. In the one method, charts are used showing the object 
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and the stars in its immediate neighbourhood. This method is usually quite 
satisfactory but when many objects, widely spread over a large area of the sky, 
have to be identified and when moreover the background is crowded, many 
separate charts are required. The other method is followed by the Harvard 
Observatory in the Magellanic Clouds. It defines the relevant object by its 
rectangular coordinates x and y (see above). ‘There is little doubt that the 
method works quite satisfactorily both in assigning x and y to a new object, and in 
finding the star when x and y are given, when the original plates, with the 
imprinted mesh system, are at hand. However, practice shows that the investi- 
gator who has no access to these original plates is confronted with almost 
insuperable difficulties in the same problem. For this problem, then, the investi- 
gator could until now make use of: (1) the two published reproductions of both 
clouds as given in H.A., 60; (2) Tables I and III, also of H.A., 60; and (3) the 
HDE which gives for a large number of CPD stars Harvard x and y in the LMC. 

Though the reproductions mentioned above show the reticule lines clearly, 
they are useful only for the identification of the brighter objects in the least 
crowded areas. 

The use of Tables I and III of H.A., 60, for the identification on photographs 
obtained with large telescopes by interpolation, fails on account of the small 
number of these standard stars available on the same photograph. 

The HDE can be used in identifying LMC objects by interpolation between 
the CPD stars with known coordinates x and y. But this cannot be done in the 
SMC and in those regions of the LMC which are outside the HDE. 

The formulae presented in the present paper in Tables I and II permit 
(among other things) the calculation of Harvard x and y from equatorial coordi- 
nates and vice versa for any object for which the one set of datais known. Since 
reproductions of photographs obtained with even the larger telescopes usually 
contain a sufficient number of CPD stars, it is now easy to rule these prints in x 
and y. In this way it has been found possible to identify objects given by their 
x and y in any part of the Clouds. 

Derivation of formulae.—For 28 AGC stars near the Small Magellanic Cloud, 
Harvard coordinates x and y are given in Table I, p. 89, of H.A., 60. The table 
contains, though in an implicit manner, the most accurate definition of the x, y 
system of the SMC available. 

Standard coordinates were calculated for the 28 stars mentioned above, using 
the 1875-0 positions as given by the AGC, 

Then least square solutions were made of the following form :— 


af +byn+e,=y—7 
a,x + byy +eg=x-—€ 


axthy+q=y—y. 


The results are presented as formulae (5), (6), (7) and (8) in Table I of the 
present paper. As a further check formulae (7) and (8) were derived from (5) 
and (6) by transformation. 

In the course of this work it was found that the data of H.A. 60 (Table I, p. 89) 
were discordant for the stars AGC 193 and 1449. They were excluded from the 
solutions. In a similar manner, the data given in Table III (p. 98) for 21 AGC 
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stars near the Large Magellanic Cloud were treated. The results of the com- 
putations will be found as formulae (5), (6), (7) and (8) in Table II of the present 
paper. As in the case of the SMC the results for the LMC were further checked 
by transformation. 

The discordant star AGC 6746 was excluded from the solutions. 

The residuals in the various solutions for both the SMC and the LMC are of 
the same order of magnitude. The mean error of a single equation of condition 
was found +2” generally. 

Formulae (1)-(4) were adapted from the formulae as given by Smart (1). 
Formulae (9)-(12) were derived from the formulae (1)--(8) by substitution and 
rearrangement. 

The formulae as presented in Tables I and II were all checked on several stars 
having known equatorial, standard and x,y coordinates. The standard co- 
ordinates of these latter stars were computed in a manner independent of formulae 
(1) and (2). 

It may be found useful to note that formula (1), (2), (9g) and (10), (3) and 
(4); and (11) and (12) have the same denominator. 

The formulae transform to an accuracy of about one second of arc. 

Preference of standard coordinates over Harvard x and y.—The question 
arises what the uses of these coordinate systems for future work could be. To 
dispense with rectangular systems altogether, and to confine oneself to the sole 


Tas_e I 


Small Magellanic Cloud 


10° sin Ax 
1°39721 cos An+ 4°64244 tan (— 8) 





_ 957571 cos Aa— 288196 tan (— 8) 
1°39721 cos An+ 464244 tan (— 8) 
g 
59445 +0°957579 
197513 —0°288207 
59445 +0°957577 
x= +12740 + 1:00694£ + 0°001087 





tan Aa= 





tan (—8)= 


cos An 





y= + 10423 —0°00135£ + 1006777 
&= — 12642 +.0°99312x —o'00108y 
2 = — 10370 + 0'00133x +.0°99332y 
_ 1007116 sin (Ax +4™175-2) + 58834 tan (—8) 
bk 1°39721 cos Aa +4-64244 tan (—8) 
pe 978618 cos (An +o0™-198-0) — 241759 tan (—8) 
¢ 1°39721 cos Aa + 464244 tan (—§) 








— 12642 +.0°99312x — 0:00 108y 
49515 +O°OO127xX+O°90SIISY 


tan Aa= 





200502 — 0'00038x — 0°28626y 
49515 +o°oor27x +0051 IS 


tan (—§)= cos Aa 
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Taste II 
Large Magellanic Cloud 
10® sin Aa 
1°73742 cos An +4°52612 tan (—8) 


_ 933580 cos Aa ~ 358368 tan (—§) (aL) 
7 1°73742 cos Aa + 4°52612 tan (—8) 


é 
tan Aan= L 
73919 +0°93358y Gh) 
192565 — 0°358377 
cos Aa 
73919 +0°93358y a) 
x= +12140+ 1°00642£ +.0°011047 (sL) 


y= +10610—o0°01111€ + 1006317 (6L) 


f= 





(1L) 








tan (—§)= 





£= — 11945 + 0°99349x — OOI0gTY (7L) 
n= — 10676 + o'o1096x + 0°99362y (8L) 


_ 1006910 sin (Aa +7™08*-9) + 50991 tan (— 8) (oL) 
oe 1°73742 cos Aa + 4°52612 tan (—8) : 





, . 957969 cos (An+ 2398-5) — 312607 tan (—8) 


L 
1°73742 cos Aa+4°52612 tan (—8) =: 





tan Aa = ogee Oa cae (111) 
63952 +0°01023x +0'92762y 





196391 —0°00393x —0°35608y cos Aa (12L) 
63952 +0°01023x +.0°92762y 


tan (—8)= 





use of an equatorial system, might seem attractive were it not that, due to the 
curvature of the mesh lines of that system, it is rather inconvenient to rule prints 
accurately. Nor is it, for the same reason, so convenient and accurate to make 
measurements of position as a rectangular system. If then, for these reasons, it 
is advantageous to continue using a rectangular system, the author favours the 
use of the systems of standard coordinates in preference to the x, y systems 
because of the greater simplicity in the relations with «, 5 of the €, 9 (see Tables I 
and II of the present paper). Since transformation between «, 8 and the rect- 
angular coordinates and vice versa may become a frequent occupation, this — 
choice in favour of the standard coordinates may merit serious consideration. 

It is a pleasure to thank Mr D. Lodge for his efficient help in calculating 
standard coordinates for a number of stars. I am indebted to Mr J. v. B. Lourens 
of the Royal Observatory, Cape, for advice on the approximate methods (2) in 
working out standard coordinates. 


Radcliffe Observatory, 
Pretoria: 


1959 February 24. 
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A SURVEY OF RADIO SOURCES AT A FREQUENCY OF 159 Mc/s* 
D. O. Edge, 7. R. Shakeshaft, W. B. McAdam, }. E. Baldwin and S. Archer 


Summary 


The Cambridge four-element interferometer has been used at a frequency 
of 159 Mc/s to determine positions and flux densities for 471 radio sources 
lying between declinations —22° and +71°. Information concerning the 
angular diameters of some of the brighter sources has also been obtained 
and the majority of these have diameters less than 6’. 

Most sources have an isotropic distribution but there is a concentration 
of the intense ones towards the galactic plane. ‘There is also evidence for an 
excess of sources in the region of the belt of background radiation crossing 
the galactic plane at /=o°. 

If the number-flux density distribution is compared with that expected 
from a uniform distribution of sources in space, a deficit of the more intense 
sources is found. 





THE POSITIONS, FLUX DENSITIES AND ANGULAR 
DIAMETERS OF 64 RADIO SOURCES OBSERVED AT A 
FREQUENCY OF 178 Mce/s* 


B. Elsmore, M. Ryle and Patricia R. R. Leslie 


Summary 


Measurements have been made at a frequency of 178 Mc/s of the position, 
flux densities and angular diameters of 64 radio sources between declinations 
+65° and —06°, both to provide a system of reference sources and to aid in a 
further search for optical identifications. 

The observations were made with a transit interferometer in which one 
aerial had dimensions 1450 ft x 65 ft; this element was separated by 2570 ft 
from a second smaller aerial mounted on a North-South railway track 1000 ft 
in length. In the present measurements, observations of each source were 
made with eight different positions of the movable aerial and the declination 
was derived from the relative phases of the interference patterns recorded 
with the different positions. The right ascension was determined in the usual 
way from the mean of the observed phases. 

The average probable errors of the 64 positions are + 18-5 in R.A. and 
+15 of arc in declination. The optical identification of six of the sources 
has been confirmed by these observations. 


* The full text of these papers is published in Memoirs of the R.A.S., 68, Part II, 1959- 
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PE3IOME JIOKJIAIOB, B TIEPEBOJE HA PYCCKHM A3bIK 


IBYOKEHHA B COJIHLIE HA YPOBHE ®OTOC®@EPHI. 
VI. BPAITEHHE COJIHLIA HW ®OTOC@EPHAA LIMPKYJIALMA 


I. X. Tlaacxemm 


Ilo cnhekrporpaMMaM, MOJIY4CHHBIM 21-23 ampena 1953 r., ObLIM M3MepeHbI B 332 TOUKAaX 
BHMMOl NOBepxHOCTH Tina 

K 
3afaty u K ceBepy. Pemenua mo cnocoOy HaMMeHBIIIMX KBafipaTOoB TaKHx 
ypaBHeHH [JIA pa3JIMGHbIX WIMPOT yKaSbIBaloT Ha aCHMMCTPHIO ee 
SKBaTOPH@IbHOM 30He HW Ha HasIMUMe MCDHIMOHAIBHbIX TeucHHil. QBWKCHHA YKa3bIBaloT 
Ha CYLeCTBOBaHHe B (poTocibepe UMPKYJIANMH, aHalOTHuHOl WMPKYIAUMM, MMCIOUIclicaA B BepxHell 
atmoctbepe Gemam. Ora rumoresa paccmaTpMBaeTcaA B cBeTe uccie¢gzoBaHua Jlaiirxunna, u 
mpeasiaraetcaA HaOmofaTeIbHbIi MeTO Ce MpOBepKH. 


MHTEHCHBHOCTb COJIHEGHOrO H3JIVYEHHA Y 114 
®, Caitiedu u P. M. Tydu 


HaGmoyenme M3IyYeHHA WeHTpa couHeuHOrOo AMCKa BOIM3H A=II,Iow (901 cm~') myTem 
HeMoOcpeACTBeHHOTO CpaBHeHHA C M3JIyGeHHeM YepHOrO Tena mpH 1300 °K Gpino 
BecHoH 1958 r. MinTeHCHBHOCTh COuIHeUHOrO H3TyueHHA Sblla UpHOspKeHHO peyuMpoBaHa K 
H3JIV4UCHHIO YEPHOTO Tesla IPH MOMOUIH Bpalljaioulleroca CeKTOpa, H OTHOMICHHe MHTCHCHBHOCTeH 
ObLIO 3aTeM H3MEpeHO IPH MOMOLIM KasMOpoBaHHOrO sIeKTpHUecKoro ocnaGuTeA, T.c, CHCTeEMa 
6bima He3aBHCHMa OT MHeiHOCTH mpHeMHMKa. EyMHCTBeHHbI¢ BCIOMOTaTesIbHbI¢ M3MepeHHA, 
KOTOpbie TpeGoBasIMCh, 3TO aGcomoTHOM OTpaxKaTenHOH cnocoOHocTM OfHOrO sepKasa 
H H3Me€HCHHA OTpaKaTesHOM CnOcoOHOCTH C H3M¢CHCHHeM yTJla OTpaoxKeHHA Apyroro sepKasia. 

Ilupuna HICH 0,031 MIM 2,5 cm~! Gblia JocTaTOuHO Maza, UTOOLI BITYCTHTS THIS aTMOC- 
dbepHoe “ oKHO”, QJIA KOTOPOrO HOrsolleHHe ONMChIBaeTCA 3aKOHOM JlamGepTa, 4TO NosBONAeT 
MPOBeCTH SKCTPanONAUMIO K MHTCHCHBHOCTH BHe SemJIM. 

CpeqHee B3BelI¢HHOe M3 NATH H3MepeHHi UcHTpaNbHO MHTCHCMBHOCTH COCTAaBJIACT 


IgG1, 10 1) =2,408 + 0,016 x 10" spr/cm* crepagman . cek 
(B ICM MHTepBaya (JIMH BOJH). SOTO SHAaYCHMe SKBMBaJIeCHTHO cH TpasbHo TemMepaType 
Toft 1,10 4) = 5036 + 30 °K. 


MATHUTHOE HM THMHAMUMUECKOE ITIOJIA BHE MIPOTOSBESIIbI 
JI. Mecmea 


Coxnmalouasica MmpoTosBesyza, opMHpyloulaica B MAaFHHTHOM OGsaKe, 
MarHHTHbI MOTOK, €CJIM TOJIbKO MJIOTHOCTS TWia3Mbl H¢ CTaHOBHTCA CJIMIIKOM HM3KOH. B stolt 
cTaTbe HCCyIeslyeTCA CTpyKTypa BHeLIHero NOJIA HM CTO BUIMAHMe Ha TeueHMe BeulecTaa. I[Ipeqno- 
JlaraeTCa, 4YTO WeHTpoGerkKHad CHUIa BCe BPeMA NpeHeOperxKHMO Maia, 2 3Be31a B OO Make HeNOABwKHA. 
Tlone B 3Be3e HM Ha GecKOHeUHOCTM CUMTaeTCA ——* NOTOK——H30T¢pMMU4CCKMM, 

Bo3MOoKHbIe¢ paBHOBECHbi¢ NOJIA BKINOUAIOT B CeOA: 

1. TlonomfambHoe Mose, CO BCCMH CHJIOBbIMH JIMHHAMH, YXORAWIMMM B GecKOHeETHOCTS, HO C 
HeOOJIbLIMM TOPOHaJIbHbIM KOMIOHCHTOM, JOCT@TOUHO CHJIBHbIM, YTOObI COXpaHHTS 
ioe cia6biM HECMOTPA Ha BTeKaHHe BellecTBa. DTo none cabo Ha SKBaTOpe, HO CHAOBBIC 
JIMHMM MIpeTepiieBaloT TaM pe3KH¢ H3BMBbI, Nopowkyad CcaMO*OKYCcMpyIOUMiicA SKBaTOpHaNbHbIt 
aspaAr. 
‘ : IlonomwyasbHoe Mose C 4YaCTLIO THHHM KOHCUHOM MpOTAYKCHHOCTH, OTICJICHHLIX KOJIBIIOM 
HeiiTpasIbHbIx ToueK Tuna X oT GecKOHe4HBIX JINHMM NOA, KOTOpoe MMeeT CTpyKTypy, Nog_oGHy10 
NoOsIKO, OITMCaHHOMy B I. I. 

3. Ckpyuennoe 6e3 cKpyamBalomlero ycHMA um mouTa GeccumoBoe mone, orpaHiueHHoe 
KOHe4HOM NoBepxHoctsi S, u nourn GesBuxpeBoe none BHe S. 

Ecnm none o6maka MepBoHayasibHO COsepxKHT TOpOMMabHbIli MOTOK, NpOHNSbIBaloumi 
MOJO TaJIbHbIC Ss 10. SOOT ee B NPOTHBHOM CJlygae THITbI I Mm 2 GynyT 
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3AMEUAHHE O MAFHHUTHOM TOPMOXKEHHMH BPAIIAIOWEHMCA 3BE3IJIbI 
JI, Mecmea 


IIlpom3BosbHOe MOIOMAabHOe MarHHTHOe Moe, coeqHHmIoUlee GOBICTpO Bpalllaloilyioca 
3Be3ay C MejIeHHO Bpatlaloulelica OGonouKOM, GymeT Tak 3aKpyueHoO, 4uTo GSyyeT MepeHOCHTE 
yrioBoll MOMeHT 4epes IPpOMExKYTOUHY1O OONacTb C HU3KOM NoTHOCTEIO. [lone 6ygeT cosqaBaTS 
paBHbie H TIPOTHBOMOJIOZKHO HalipaBJICHHbi¢e Mapbl CHJI, IIPWJIOXKCHHbIC COOTBCTCTBCHHO K 3Be3]1¢ 
vu K OG6ono¥Ke. B mpomexkyTouHoi oGnacTH Mose, BepoATHO, He GygeT co3qaBaTb GoubUIMX CHI 
oxaTua u 6GygeT upHOnmKeHHo GeccunoBbim. B ormmune oT adcosiomuo GSeccHoBoro Nona, 
paccMOTpeHHOro B mpefbrayuielt paGoTe aBTopa, KOTOpoe He MOKET BOSHHKAaTb B pe3ylIbTaTe 
necopMallMM NOOK AaIbHOrO MONA, B JaHHOM CJIyyuae He Gy eT BOSHMKAaTb TOPOMAaJILHOrO MOTOKa, 
3aMbIKaloulero MOJOMMaIbHy!O MeTINO, a BIeKTpHuecKHi TOK He GyfeT TeUb B TOM 2Ke CaMOM 
HallpaBJICHHH BO BCeCxX TOUKaX TIeTJIM. 


3BESIIHbIE IPYIIIbI. IV. fPYOMMA 3BE3], C BOJIBUIMMH CKOPOCTAMH 
I'PYMBPHIJDK 1830 H EE CBA3bC IWIAPOBbIMH CKOIMJIEHHAMH 


O. JIow. Seeen u A. P. Candudo 


Mmeroulmeca aHHble 0 COOCTBeHHBIX JBHKCHHAX HM JIY4YeBbIX CKOPOCTAX OBLIM HCIOJIb3OBaHbI 
JIA WOKasaTesIbCTBa CYICCTBOBaHHA JBMwKylleliica rpymnbi cyOKapmHKoB (rpynma 3Be3q pym6- 
pup 1830), BKuOuaoUlel sBesqu RR |.yr. Ha ocnHopanuMu cyulecrByloluleii saBMCHMOCTH 
M@Ky HaGJOJCHHbIM YJIbTpa (HONeETOBLIM H3ObITKOM H CMCIIICHHeEM H¥DKe HOPMAJILHOH riaBHoi 
MocueOBaTesIbHOCTH CyOKapsIHKH B rpynme [pymOpuypK 1830 GbLIM OTOXKTeCCTBIICHbI KaK 3Be3bI 
riaBHO MOCIeOBaTesIbHOCTH WiapOBbIX CKOMIeHHH. OTo OTOMKReCTBMeHHe IPHBOAMT K 3Haue- 
HHIO MonyIa m—M=14™,2 ana maposoro cKkommeHHa M13, B pesymbTaTe “ero OKa3bIBaeTCA 
4uTO (JIA JByx MepeMeHHprx Tuna RR Lyr My~+o™,s, a ana HamOonee APKMX 3Be3q, CKOILIe- 
Hua, My=—2™,3. Tlokasprpaetca, uro cama 3Be3qa RRLyr, Ana KoTOpoH Ha OCHOBe 
rpymmoporo Mapasiiakca MbI uMeem My~-+0™,8, yyospneTBopseT COOTHOIIeHHIO nepHoj- 
aMILIMTya iid MepeMeHHbIX B cKOMIeHHH M3 Hu TO OHA HMMeeT LBeCTOBOM H30bITOK 
(moKpacHeHHe), paBHbIii OM,o5 mo OTHOMeHMIO K 9THM 3Be3qamM. IIpHpaBHMBaa CBeETHMOCTb 
RR Lyr K cpeqHelt cBeTHMOCcTH NepeMeHHBIX B M3, MbI ouyuaeM (JIA MOYIA CKOIVICHHA 
sHaueHne m—M=15™,o. Mobi He Mpom3BemM ompezeneHHe MOfyJIA OOBIUHBIM JIOrHUCCKHM 
cmoco6omM IyTeM COBMeUIe¢HHA TiaBHbIX mMociteqoBaTembHocteli B M3 u rpymme IpymOpnmpx 
1830, Tak KaK MOKa3aTeIM IBeTa 3Be3 TiaBHOM MocneyoOBaTesbHOCTM B M3 MoryT cojepxKaTS 
CHCTeMaTHYeCcKy10 OWHOKy. 

Tak Kak NIPHCyTCTBMe MepeMeHHbIX Tuna RR Lyr B 3Be3{HbIX rpyllax MOMKET WaBaTb e{MHCT- 
BeCHHY!O BO3MO*KHOCTS TOUHOTO ONpeesICHHA CBETHMOCTei STHX 3Be3, IpesCTaBIAeTCA BaxKHbIM 
MIPOH3BeCTH CHCTeMaTH4eCKHMe NMOMCKH TaKHx rpymm. Mbi mpypenm B Tammie MMeIolmecs 
JJaHHbie OTHOCHTEeJIBHO (1) BCeX 3Be3jl, OTHCCEHHbIX K CyOKapJIMKaM Ha OCHOBaHHM OcOoGeHHOCTeli B 
MX CII¢KTpax, (2) 3Be3q, KOTOPbIe, CyqA NO MX HaOMIOJCHHBIM yIbTpadHoueTOBbIM M3ObITKaM, 
110 BCeli BCPOATHOCTH ABJIAIOTCA KpaiiHuMu CyOKapsMKaMH, XOTA H CUHTAIOTCA OObeKTaMM rlaBHoi 
Mocuse{OBaTeNbHOCTH, M (3) 3Be3, KOTOpbie MoryT GnITs mH60 ruraHTamu, 2HG6o cyOruraHTaMu, 
mH60 3Be3aMH TOPHSOHTA@JIbHOi BeTBH wWapoBbix cKomeHuii. IIpocrpaHcTBeHHbIe CKOPOCTH 
3Be3, STHX TpexX KaTeropHit ObLIM BbIUMCIICHbI C MapaliakcaMH, ONpej[eJICHHbIMM IIpH NOMOLIM 
(My, B—V)-quarpamMp! Ju1a cKoneHua M13, npoKasMOpoBaHHoii ucnombsya rpynny pym6puwpxK 
1830. TIlokasaHo, uro c Gombmoi cTeneHbIO BePOATHOCTH HeKOTOpOe UYMCIIO STHX 3Be3q 
o6naylaeT OOLIMM JBYOKeHHeM B IIpOCTpaHCcTBe H YTO C IPHBeUeHHeM JONONHMTeIBHOTO MaTepvasia 
yaaetca OOHapyxKHTb pay rpynn, aHanormuHEix rpymme TpymOpuypx 1830. Tak xe TOUHO 
MOKa3aHO, 4TO MPOCTpaHCTBeHHbI¢e CKOPOCTH “YeTbIPex NepeMeHHBbIX THMa RR Lyr C OTHOCHTEeJIBHO 
XOPOIIO OMpeesieHHbIMH COGCTBCHHbIMH JBHOKCHHAMHM COBMagaloT C JBH)KCHHAMH HCKOTOpPbIX 
OKpecTHBIX CyOKapmuKoB, uTo nano My~+o™,6 ana sTux “eTbIpex MepeMeHHBIX. 


O CYIJECTBOBAHHUW CYBKAPJIMVKOB HA JIMAFPAMME (Mg,;, log T;) 
A. P. Candudoe u O. JIoe. Dezen 


PaccmaTpuBaeTca BIIMAHHe OcIaOneHHA cbpayHTodeposbix HHH Ha HaOMOMaembie MOKasaTe1H 
peta 3Be3n. Yurenbi 3ajlepoKKa HempepbIBHOrO H3JIYYeHHA JIMHHAMH HM OOYCNOBIeHHO? JIMHHAMH 
H3MeHeHHe TeMMepaTypHoro pactipesenenua B comHeuHol armocdepe. Ecim 6nr pce MCTOUHMKH, 
cosfaIoumme JIMHHM, OLLI yaanenbi c Comma, TO H3MeHeHMe HaGsHOmaeMbIx NoKasaTene WBeTa, 
cormacHo TeopuHm, cocrapiano 6h1 A(U—B)=o0%,32 un A(B—V)=0™,17. Oru m3aMcHeHHA 
onmpenenmot “ mMHM0 TermMuHoro sexta ” (“ blanketing line”) Ha qHarpamme (U—B, B—V), 
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M3 KOTOpod nompaBKa K HaGmOqaemomMy NoKasaTemo ypeta B—V, ob yeon tena 
JIMHHAMM, MO>KET GbITS HaiifleHa 43 HaGOAaemoro yrbTpaduoseroBoro us6nITKa 6(U—B). 

IlpuBeyeHbI sBe3qHbI¢ BeEIMUHHLI M Beta cyOKapmMKoB, Nexauyx Gonee uem wa 0,25 
HIDKe HOPMaJIbHOH riaBHo nocnegopatenbHocTw. TlompaBKu 3a UBeT, BbIBCCHHBIC B32 
aeMbIX YJIbTpa*bMOeTOBbIX M3ObITKOB, NPHaHbI K HaGmMOTaeMbIM MoKasaTenam upeta B—V. 
OTH NOMpaBKH ABMAIOTCA JOCTaTOUHbIMH, 4TOObI NepesqBHHYTS yYMepeHHEie cyOKapmHKM Ha 
rilaBHyl0 nocueqoBatembHocTs Tuay. Jienaerca BbiBog, uTo SGombumMEcTBO CyOKapJIMKOB Haxoj- 
HTCA HWDKe rlaBHOM MOcue_OBaTebHOCTH Ha JMarpamme (My, B—V) Tomko Bchea_CTBHe BRMAHMA 
cna6bix cbpayHrodeposBbix HHH Ha ux HaGmogaeMbIe BETA. 

Uniexbi (BwxKylleiica rpynne pym6pyypx 1830 He NOHOCTBIO BOCCTaHaBIMBaloTCca B riaBHolt 
NOCI€OBaTeIbHOCTH HM, BO3MO)KHO, YTO 9TH 3Be3bI WelCTBHTeMBHO NeKaT HIDKe rnaBHoll 
TOCJI€¢OBaTeIbHOCTH Ha AMarpammMe (Magol, log Te). Tlpeqpaputenbubie 3aKOUCHHA M3 STHX 
jaHHBIxX cuegyioume: (1) HempospauHocTs cyOKa ckopee, no-BHMMOMy, o6ycnOBNeHa 
cBo602HO-cBO6O2HBIMH nepexogzamu H mu He, em NOTIOUCHHeM TADKeSIbIMK 
auieMeHTaMH; (2) oTHOmIeHHe H/He ana 3Be3y rpymm: [pym6pug* 1830 mMomeT HeMHOrO OTIH- 
YaTLCA OT AaHaJIOTHUHOPO OTHOMICHHA JIA sBe3q cKONIeHHA Duan. 

Tennuunasn (blanketing) Mofens aeT COOTHOMIeCHHe MOKy PpasHOCThIO BextHUHH 
m060H 3Be3qbI OT riaBHOM nocneyoBaTembHocTH [vag u HaOmo7aemMEiM YO-melatanen Beta. 
STO COOTHOUIeHHe NPHMeCHACTCA K 3BE31aM B PpasJIMUHbIX raylaKTHYCCKHX CKONMJICHMAX H K 3Be3 aM 
o6ulero NONA, OTKyTa MOUWy4sIOTCA pesyNbTaTbI: (I) METO NOMTOHKH riaBHbIX NocHesOBaTeL- 
HocTeili JIA ONpeyeueHHA MOJYIA paccTosHHi ramakTMueckmx CKonIeHH ommOoueH, ecm He 
BHECeHbI NompaBKM 3a pasymuHbIe YOD-nsGnirKu; (2) HaGsnOTaeMbIii pas6poc B O10 B 
Yo- u3ObITKAX 3Be3f, NONA, KOoTOpbie O6LIUHO HasbiBaloTca “ KapMKaMM rnaBHOl nocnreqOB- 
aTeJIBHOCTH ”’, NOKa3bIBaeT, “TO NMOIHaA HMCTHHHAA WMpHHa TiaBHO nmocieqoOBaTesBHOCTH, 
onpeylenaeMad STHMM 3Be31aMH, CocTaByaeT NO Kpaiineli Mepe O™,60. Tak Kak 3Be3qBI NOAA C 
CHJIBHbIMM JIMHMAMM, CyOKapJIMKM HM 3Be3{bI pasIHUHBbIX PavlaKTHYCCKMX CKOMJIeHHH 
KOHTHHYYM Ha Juarpamme (B—V), . TO MpecTaBAeTCA BepOATHBIM, 4TO (a) cyGKapuMKH 
o6pasyi0oT KOHTHHYyyM Ha Juarpamme (My, V), o6benmnancs c Gonee HOPMaJILHBIMM SBesamMn 
C CHJIBHBIMM JIMHBAMH, u (6) cytecrBy1oT beta A MO@K]Iy COjlepyKaHHeM TSDKCIIbIX SJIEMCHTOB B 
OTICNBHbIX 3Be37aX NOJA ¥ B SBe3Tax pasMUHBIX ralaKTHYeCKHX CKONeHHH. 


PATIMOHU3JIYYEHHE HOPMAJIBHBIX TAJIAKTHK. I. HABJIIONEHHA M31 HW M33 
HA UACTOTE 158 MW 237 MIL. 


P. X. Bpayn u C. Xasapd 


Bpuio MpoBezeHO M3y4eHHe paqMOMSJIVYeHHA APKMX HOPMAJILHLIX TawakTHK Ha YaCcTOTEe 
158 Mru. B gannoii crarse coo6imaeTca o pesymbTaTax, orHOCAUIMKXcA K M31 m M33. Maurer- 
PAIbHble PayMOBeIMUMHD! (mr) STHX WBYX TalakTHK Obl M3MepeHbI HM CpaBHeHbI Cc mx doTor- 
pacbHuecKMMM BesIHUMHAaMH (mp); HalifleHO, YTO pasHOCTS (mr—mp) NOYTH OFMHaKOBa Ain OGenx 
rayiaKkTHK. 

PacrspeyeneHue MHTeHCHBHOCTH NonepeK M31 yKasbipaeT Ha HalIMUMe AByX TNaBHbIX KOMIO- 
He€HTOB: JMCKa HM KOpoHbI. Kopona o6ycnosnmpaet 90% oGuiero usmyueHMA, HMeeT OTHOIMCHHe 
BHJIMMBIX Oceli, paBHoe NPHMepHO 0,6; ee GombWad OCb MPMMepHO B TPH pasa AmMHHee Gonnmolt 
OCH BHJMMOi TyMaHHOCTH. PesymbTaTbI NOKasbIBaIOT, YTO M33 TaloKe HMCeT KOPOHY; asIHTI- 
THUHOCT, MSIYYeHHe HM pasMepbl KOPOHbI NO OTHOMEHHIO K BHIMMOM TyMaHHOCTH BecbMa GnMSKH 
K 3Ha4eCHHAM, HalijeHHEIM Wa M31. 

CpapHenne M31 u M33 c nameit Tanaxrukoli ykasbipacT Ha TO, 4TO H3NyYeHHe KOPOHBI 
TanaKTHKH CylUleCTBeCHHO CHJIbHee, 4M H3y4eHHe KOPOH YKa3aHHbIxX ABYyX TasiakTHK. 


KOCMOJIOrMA CTAIMOHAPHOTO COCTOAHMA, PACCMOTPEHHAA 
COrJIACHO OBIEH TEOPHH OTHOCHTEJIBHOCTU 


B. [[siieudcon 


Mutepnpetalima KOCMOJIOrHH cTaliMoHapHOrO Salen meanness "aren aa 
OTHOCHTCJIBHOCTH, Ha¥aTad Max-Kpu, B ganHol paGoTre mpofomxena 
BHelliHeH HM BHYTPpeHHeli MCTPHK UIA CTaTHUeCKOrTO HaOO*KeHMA chepueciar commie 
CKOHIC¢HTPHPOBaHHOM HIM paspeoxxKeHHoll Macchi Ha “ cralMOHapHy10 ” 

Haiigzeno, YTO NOCTOMHHaA TIONH@A TWJIOTHOCTE p,» M OTpPHlaTebHOS JABNEHHe po (ok 
nl0-BHMMOMy, NpeyctaBsiaror cool xapaKTepHCTHKM €CTeCTBEHHOTO COCTOABHA 
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ycroitumBocTw cTaljMoHapHolt cpenbr. Mccnegopanue moKa3biBaeT, YTO 39TO COCTOAHMe MODKET 
CYUICCTBOBaTh MOK Ay ralakTHKaMH B CTalMOHapHoH Moyen Bcenenuolt. Paccmotpenue MMpOBbIXx 
JimMHMa CpeqbI B OKPeCTHOCTM CTaTHUeCCKOrTO BOSMYIICHHA MPHBOAMT K BbIBOLY, 4TO NoROGHOE 
BosmyuleHHe O6nagaeT cnocoGHocTBIO 2HG60 MpOTHBOeHCTBOBaTh KOCMHYCCKOMY OTTAJIKHBaHHO, 
mm60 YCHJIMBaTb ero B 3€BMCHMOCTH OT TOFO NpeCTaBJIAeT JIM OHO KOHICHTpallMio WIM paspexKeHHe 
MaTepHH MO OTHOIICHHIO K ecTecTBeHHOMYy cocTosHHIO. IIpeymaraeTca MexaHH3M 3apooxKqeHHA 
MaTePHH, IPH KOTOPOM CKOPOCTb 3apO?KCHHA CHJIBHO 3aBMCHT OT JIOKaJIbHOrO rpaBHTaljMOHHOrO 
BIIMAHHA; OHA OKaSbIBaeTCA MOJIO*KMTebHOM TaM, re cCpesa paclIMpAeTcA, HM OTPHUaTeNbHOlt 
(aHHMTWIAWMA) TaM, rye OHA CKMMACTCA. 

OnpeyeneHb!I mpejesibI pasMepoB HM IWIOTHOCTH JIA CTaTHYeCKHX BOSMYLIICHHH ecTecTBeHHOrTO 
coctosHua. TIlokasaHo, 4TO B BOSMYIIeHHOM cpeye 3HaK rpaBHTallMOHHOrO NOJA 3aBMCMT OT 
JIOK@JIBHOrO COCTOMHMA CpebI HM B Cuy4ae pasperxKeHHA TWIOTHOCTh MHEPIMH, B HbIOTOHOBCKOM 
CMBICII€, OKa3bIBaeTCA OTPHaTebHO. YKasaHbI TeopeTHueCcKHe YCIOBHA, HEOOXOZMMBIC WIA 
OTHOCHTeJIBHOM yCTOHMYMBOCTH CCTECTBEHHOrO COCTOAHHA. 
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ycTOMUMBOCTH CTal[MOHapHoli cpeybi. MccneyqopaHue moKa3bIBaeT, YTO ITO COCTOAHMe MOKET 
CYUI€CTBOBaTb MCK Ly raslaKTHKaMH B CTalHOHapHOH Moje1H BcenenHoi. PaccmorpeHHe MHpoBbIx 
JIMHHH CpesbI B OKPeCTHOCTH CTaTHYeCKOrTO BOSMYLIICHHA MIPHBOMT K BbLIBOY, 4YTO NMOWOGHOE 
BOSMyleHHe OOsIafaeT CnOCcOOHOCThIO J1H6O NPOTHBOJeCHCTBOBaTh KOCMHYCCKOMY OTTAJIKHBaHHIO, 
nm60 YCHJIMBaTb ero B 34BHCHMOCTH OT TOPO Mpe/{CTaBJIACT JIM OHO KOHUCHTpallHio HIM paspexKeHHe 
MaTePHHM TIO OTHOUWICHHIO K ecTeCTBeHHOMY cocToAHHW. IIpemnaraetcA MexaHH3M 3apoKeHHA 
MaTepHH, IPH KOTOPOM CKOPOCTh 3apOK[CHHA CHJIBHO 3aBHMCHT OT JIOK@JIBHOrO rpaBHTallMOHHOro 
BIIMAHHA; OHA OKa3bIBaeTCA MOJIO*KHTCIbHOH TaM, rae Cpeya pacluMpxAeTcH, HM OTpHUaTesbHOit 
(aHHHMPHJIAWMA) Tam, re OHA CKMMAaCTCA, 

OnpejleseHb! mpesesibI pasMepOB HM MWIOTHOCTH JUIA CTaTH4eCCKHX BOSMYLICHHH ecTecTBeHHOrO 
cocronHus. IlokasaHo, 4uTO B BO3MYLIICHHOM cpejle 3HaK PpaBHTal[MOHHOrO MOJIA 3aBHCHT OT 
JOKaJIbHOrO COCTOHHHA CpejibI HK B Cly4ae pasperxKeHHA IWIOTHOCTb HHEPUHH, B HbIOTOHOBCKOM 
CMbIC/I€, OKaSbIBaeTCAH OTPHaTesbHOH. YkKasaHbl TeopeTH4ueCKHe YCNOBHA, HeOOXOAMMBIe JIA 
OTHOCHTCJIBHOH YCTOHUYMBOCTH CCTCCTBCHHOrO COCTOAHHA, 
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